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Summary 

A collection of analytical studies is presented related to unconstrained and constrained 
aircraft energy-state modeling and to spacecraft motion under continuous thrust. With 
regard to aircraft unconstrained energy-state modeling, the physical origin of the singular 
perturbation parameter that accounts for the observed two-time-scale behavior of aircraft 
during energy climbs is identified and explained. With regard to the constrained energy- 
state modeling, optimal control problems are studied involving active state-variable 
inequality constraints. Departing from the practical deficiencies of the control programs for 
such problems that result from the traditional formulations, a complete reformulation is 
proposed for these problems which, in contrast to the old formulation, will presumably 
lead to practically useful controllers that can track an inequality constraint boundary 
asymptotically, and even in the presence of two-sided perturbations about it. Finally, with 
regard to spacecraft motion under continuous thrust, a thrust program is proposed for 
which the equations of two-dimensional motion of a space vehicle in orbit, viewed as a 
point mass, afford an exact analytic solution. The thrust program arises under the 
assumption of tangential thrust from the costate system corresponding to minimum-fuel, 
power-limited, coplanar transfers between two arbitrary conics. The trajectory equation 
describing the above exact analytic solution is identical in form with the trajectory equation 
corresponding to Keplerian motion (motion with zero thrust). This solution can be used to 
satisfy boundary conditions corresponding to arbitrary coplanar transfer and escape 
problems. The thrust program can be used not only with power-limited propulsion 
systems, but also with any propulsion system capable of generating continuous thrust of 
controllable magnitude, and, for propulsion types and classes of transfers for which it is 
sufficiently optimal the results of this report suggest a method of maneuvering during 
planetocentric or heliocentric orbital operations, requiring a minimum amount of 
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computation, and thus uniquely suitable for real-time feedback guidance implementations. 
The results pertaining to the thrust program and to the exact analytic solution of the 
equations of motion are summarized in Appendix H, and then generalized to a much wider 
class of thrust programs, given the name the "Keplerian class", in Appendix I, supplied at 
the very end of this report. It should be emphasized that the Appendices (with the exception 
of A and C) are bn integral part of the research reported in the fourth chapter of this report. 
Specifically, any reader who knows anything about Keplerian motion (spacecraft motion in 
the absence of thrust) can get a preliminary idea about the primary contribution of Chapter 
four by examining just two Tables, namely, Table H.l (page 135) of Appendix H, and 
Table 1.1 (page 139) of Appendix I. 
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Chapter I 


Introduction 

In spite of all the recent technological breakthroughs regarding the speed of 
numerical computing, and the latest advances in new computational techniques, analytical 
methods still continue to play an essential role in problems related to aircraft and spacecraft 
trajectory optimization and optimal guidance, and to flight mechanics in general. One 
reason for this is that there is always much physical insight to be gained from any analytical 
procedure that extracts as much information as possible from the differential equations 
governing a problem, even without attempting to solve them. For example, in systems 
exhibiting multiple-time-scale behavior, the fundamental physical process that gives rise to 
such behavior can sometimes be uncovered and understood by an analysis of the governing 
differential equations. A second reason is certainly that all real physical problems are first 
reduced to mathematical formulations or models before any attempt is made toward their 
solution. Formulating a real physical problem in precise mathematical language is 
equivalent to asking precise questions, and one cannot expect to obtain the right answers 
unless one asks the right questions. Analytical techniques can sometimes be useful and 
result in further progress just by helping one ask the right questions. Finally, a third reason 
for the importance of such methods is that the real-time solution of the common two-point 
boundary value problems encountered in the field of aircraft and spacecraft trajectory 
optimization still remains a dream. Accordingly, the numerical computation of optimal 
trajectories and optimal feedback guidance laws has still to be aided and supplemented 
considerably by analytical methods. 

The major objective of this report is the application of such methods to the areas of 
aircraft and spacecraft trajectory optimization related to energy-state modeling and Non- 
Keplerian orbital motion. 

A systematic approach is presented in the second chapter for identifying the 
perturbation parameter in singular perturbation analyses of aircraft optimal trajectories and 
guidance. The approach is based on a nondimensionalization of the equations of motion. It 
can be used to evaluate the appropriateness of forced singular perturbation formulations that 
were used in the past for transport and fighter aircraft. It can also be used to assess the 



applicability of energy-state approximations and singular perturbation analyses for 
airbreathing, transatmospheric vehicles with hypersonic cruise and orbital capabilities. In 
particular, the family of problems related to aircraft energy climbs is considered. For 
energy climbs constrained to a vertical plane it is shown that the singular perturbation 
parameter can logically be taken as the maximum allowable longitudinal load factor of the 
vehicle. Two-time-scale behavior is suggested when this load factor is sufficiently less than 
one. 

In the third chapter the focus is on optimal control problems involving active state- 
variable inequality constraints. Such problems arise naturally in the context of aircraft 
trajectory optimization whenever functions of the state variables (such as dynamic pressure, 
aerodynamic heating rate, etc.) are constrained along a trajectory. Traditional analytical 
techniques that work fairly well in the absence of such constraints break down when such 
constraints are present, and therefore new methods are needed for an understanding of 
these problems. Departing from this premise, a transformation technique is introduced in 
the third chapter that splits the class of all piecewise continuous (in time) controllers that 
track a given hypersurface in the state space of a dynamical system into two disjoint 
classes. The first class contains all controllers that track the hypersurface in finite time. The 
second class contains all controllers that track the hypersurface asymptotically. Four 
theorems are presented that describe the two classes. The results are applied to the study of 
optimal control problems involving active state- variable inequality constraints. The 
controllers obtained from the traditional formulation of such problems are typically finite- 
time and one-sided, that is, they break down when a disturbance throws the system toward 
the prohibited side of a state-constraint boundary. These features tend to make such 
controllers quite unattractive from a practical point of view. This report proposes a 
reformulation of such problems in which the optimization is carried out only with respect to 
asymptotic controllers. The reformulated problem leads to controllers that are 
approximately optimal, asymptotic, but still one-sided. However, if the state constraint is 
regarded as a soft constraint, then one can show that there may exist controllers that are 
asymptotic, two-sided, and result in the same optimal value of the performance index 
corresponding to the original problem, that is, they are practically optimal, but at the 
expense of violating the state constraint. 

The topic in the fourth chapter is orbital motion under continuous thrust. A 
continuous thrust program is presented for which the equations of two-dimensional motion 
of a space vehicle in orbit, viewed as a point mass, afford an exact analytic solution. The 
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thrust program is proportional to an arbitrary throttling parameter, and arises from the 
optimization problem corresponding to minimum-fuel, power-limited, coplanar transfers 
between two arbitrary conics. It is shown that, for this problem, the assumption that the 
thrust is tangent to the flight path results in the complete analytic solution of the system of 
state-costate equations governing the optimal trajectories. This approximation of tangential 
thrust is made only in the costate equations, affecting only optimality, resulting in the 
elimination of the costates, and giving rise to the thrust program, for which the state 
equations can be solved analytically with no further approximations. The most striking 
aspect of the motion suggested by this solution is the fact that it is described by a trajectory 
equation identical in form to the one corresponding to Keplerian motion (motion with zero 
thrust)! The difference is that, in the latter the angular momentum is a constant, while in the 
former it is a linear function of time, with slope equal to the throttling parameter. This 
similarity allows one to speak of the two motions by using more or less the same 
mathematical vocabulary. The transverse (horizontal) component of the thrust program is 
inversely proportional to the radial distance from the inverse-square attractive center, while 
the radial (vertical) component is such that the thrust is in the flight-path direction. The 
constant of proportionality (throttling parameter) appearing in the thrust program can be 
appropriately selected to satisfy the boundary conditions corresponding to arbitrary 
coplanar transfer and escape problems. To document these facts, and demonstrate the 
existence of such particular solutions, several examples are given of such maneuvers. 
Questions concerning the optimality of such trajectories are also dealt with, and hints are 
provided, suggesting that there should be at least some classes of transfers for which the 
resulting trajectories are sufficiently optimal, both for power-limited and constant ejection 
velocity types of propulsion. The results pertaining to the thrust program and to the exact 
analytic solution of the equations of motion are generalized to a much wider class of thrust 
programs, given the name the "Keplerian class", in Appendix I. 
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Chapter II 


Aircraft Energy-State Modeling and Singular Perturbations 

2.1 Introduction 

The methods of matched asymptotic analysis in singular perturbation theory are 
based on the presence of small parameters in the differential equations of motion, that give 
rise to multiple time-scale behavior. It has been noted by several authors 1 - 2 that, in spite of a 
wide number of papers attesting to the applicability of singular perturbation methods to 
optimization problems in aircraft flight mechanics, few have been successful in first casting 
the equations of motion in a singular perturbation form. Exceptions are Refs. 1-4. Two 
methods of analysis for time-scale separability are proposed in Ref. 1. Both of these 
methods are based on an estimation of the state variables’ relative speeds. In Ref. 2 a 
rescaling to nondimensional variables is recommended. However, it is noted that the proper 
scaling transformation is not obvious, even if the time-scale separation of the variables is 
well-known from analysis or experience. Both of these papers (and in particular Ref. 1) 
provide extensive references to earlier studies which employ so-called forced singular 
perturbation formulations, in which the perturbation parameter e, nominally equal to one, is 
artificially introduced as a bookkeeping parameter, in a formal expansion of the solution 
about e = 0. In particular, there exists a large number of publications on the optimization of 
aircraft energy climbs (see for example Refs. 5-8), none of which identify an appropriate 
perturbation parameter in terms of the relevant problem parameters. This disparity is 
particularly disturbing, especially considering the number of years that have passed since 
such analytical techniques were first introduced in the flight mechanics literature. Note in 
contrast, that there have been applications of singular perturbation theory to other areas of 
flight mechanics, in which the perturbation parameter was clearly identified in terms of the 
relevant problem parameters. For instance, with regard to aero-assisted orbit transfer, the 
perturbation parameter happens to be just the ratio of the atmospheric scale height to the 
minimum trajectory radius (see for example Refs. 9, 10). In any singular perturbation 
analysis, it is advantageous to identify the perturbation parameter £ in terms of the relevant 
problem parameters (which in general include the boundary conditions), so that the 
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physical process that gives rise to the two-time-scale behavior is clearly understood. Then, 
the range of parameter values for which the perturbation analysis is possible can be easily 
identified. In fact, knowledge of time-scale separability present in the system dynamics, 
and success in exploiting this characteristic to obtain approximate solutions, cannot in itself 
be a rigorous justification for artificially introducing £. 

In this chapter an attempt is made to partially rectify the situation described above by 
presenting a systematic (albeit still ad hoc) approach for identifying the singular 
perturbation parameter e through nondimensionalization of the problem variables. Attention 
is focused on nonlinear optimization problems in flight mechanics, though most of the 
considerations that are presented apply in other fields as well. The main motivation for 
collecting and stating these considerations is to define the thought process by which it is 
possible to arrive at a suitable scaling of the aircraft energy climb problem. Of particular 
interest, from the point of view of future potential applications, is an assessment of energy- 
state approximations and singular perturbation analyses for airbreathing, transatmospheric 
vehicles with hypersonic cruise and orbital capabilities. 


2.2 Subsonic-Supersonic Regimes, Flat Earth Approximation 


Consider atmospheric flight of a conventional aircraft, viewed as a point mass, in a 
vertical plane over a flat Earth. The equations governing such flight can be reduced to a 
three-state model in mass specific energy E, flight path angle y, and altitude h. The vehicle 
mass, m, is assumed to be constant. The equations are: 


dE V(T-D) 
dt m 


( 2 . 1 ) 


dy _( L ( gcosy ^j 
dt ImVj l V J 


( 2 . 2 ) 


— = Vsiny (2.3) 

where L, D and g denote the lift, the drag and the (constant) gravitational acceleration. It is 
assumed that the atmosphere is stationary, and that the thrust, T, is directed along the flight 
path. The specific energy (mechanical energy per unit mass of the vehicle) E and the speed 
V are related by: 
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(2.4) 



and E rather than V has been employed as a state variable. 

For a singular perturbation analysis, Eqs. (2.2) and (2.3) are commonly written as: 


M L ]-I 

f gcosy') 

(2.5) 

dt KmVj 1 

k V J 

dh __ . 

— = V sin Y 
dt 


(2.6) 


where 8 is artificially introduced, and its nominal value is said to be equal to 1.0. The main 
purpose of the present chapter is to avoid an artificial introduction of e at the outset, thus, 
Eqs. (2.2) and (2.3) will be retained, while Eqs. (2.5) and (2.6) will be used only as a 
guide for the natural introduction of e. 


2.2.1 Nondimensional Form 

The first step in seeking a natural introduction of the perturbation parameter e is to put 
Eqs. (2.1) through (2.3) in nondimensional form. To this end one may start by defining the 
set S: 

S = {t 0 ,E 0 ,h 0 ,V 0 ,T 0 ,D 0 ,L 0 } (2.7) 

The elements of the set S are at this point arbitrary positive quantities, and the only 
restriction that one imposes upon them is that: 

^ has dimensions of time 
E 0 has dimensions of energy per unit mass 
ho has dimensions of length 
V 0 has dimensions of speed 
T 0 , D 0 , and L 0 have dimensions of force 
Using the elements of S to define the nondimensional quantities: 



( 2 . 8 ) 
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(2.9) 



L = 


_L 

L 0 


Eqs. (2.1) through (2.3) can be put into the following nondimensional form: 


«v(tt 0 -ddJ^O 

dt l,E 0 mJ 


( 2 . 10 ) 


dy 

dt 



Eplp ) 

mV 0 J 





( 2 . 11 ) 


dh 

dt 


ltT inY 


( 2 . 12 ) 


The goal is now to put Eqs. (2.10) through (2.12) in the traditional singular perturbation 
form. Multiplying both sides of Eqs. (2.1 1) and (2.12) by (1^/ V 0 ^ ) results in: 




Vot-oJ 


dt 



(2.13) 


vVo. 


dh 

dt 


= Vsin y 


(2.14) 


Comparing the set of Eqs. (2.10), (2.13), and (2.14), with the set of Eqs. (2.1), (2.5) and 
(2.6), it is evident that one can make the two sets similar by imposing the following four 
conditions on the elements of the set S: 


T 0 ” T> 0 

(2.15) 

T. t.V. 

U U U — J 

E 0 m 

(2.16) 

L„h„ 

-“-f = 1 
mV 0 2 

(2.17) 

V„ ! 

(2.18) 
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If one defines £ as: 


V 0 1 0 

then, Eqs. (2.10), (2.13) and (2.14) assume the form: 

f = V ( T -D) 


(2.19) 


( 2 . 20 ) 


£ dx JL-cosy^j (2.21) 

dt 1 V j 


e— = Vsiny (2.22) 

dt 

To summarize, it was shown in the present section that it is possible to introduce a 
parameter £ naturally into the equations of motion (Eqs. (2.1) - (2.3)) by first introducing a 
set of arbitrary positive quantities S (see Eq. (2.7)) to scale the variables of interest, and 
then by imposing four conditions (Eqs. (2.15) - (2.18)) on these quantities so that the 
resulting nondimensional equations assume the traditional singular perturbation form (Eqs. 
(2.20) - (2.22)). Note that only one of the arbitrary quantities in S is uniquely determined at 
this point. Combining Eqs. (2. 17) and (2. 18) it follows that Lo is given by: 

L 0 = mg (2.23) 


2.2.2 Specifying a Particular Nondimensional Form 

As shown in the previous subsection, only four conditions are imposed on the seven 
elements of set S in transforming the equations of motion to the traditional singular 
perturbation format. This means that one can specify three of the elements of S to fit one's 
convenience and then determine the remaining four using Eqs. (2.15)-(2.18). The first 
conclusion therefore is that in general the value of £ is quite arbitrary. For example, by 
choosing h^,, V 0 , and ^ in two different ways £ can be made arbitrarily small or large. The 
separability of the time scales on the other hand is a property of the system and not of the 
particular nondimensional form of the equations of motion that is chosen. One therefore 
should expect that if the system does indeed possess the property of time scale separability. 
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it will exhibit it no matter what the actual value of e is. This is precisely the reason for the 
success of so many singular perturbation treatments of the past in which £ was introduced 
artificially and its nominal value was said to be "fixed" at one. 

Although there is no unique way of specifying a particular nondimensional form of 
the equations of motion, one can argue that there is at least one choice for the elements of 
set S that results in additional physical insight. First, in order to maintain the relationship in 
Eq. (2.4) in the transformed variables, a fifth condition is introduced: 

E 0 =gh 0 (2.24) 


which together with Eq. (2. 18) implies: 

V 2 

E = — + h (2.25) 

2 

Using Eqs. (2.16) and (2.24) in Eq. (2. 19), it follows that £ can be written as: 

£ = — (2.26) 
mg 

Now, only two among the seven elements of set S need to be specified. Then, the five 
conditions, Eqs. (2.15) - (2.18) and (2.24), uniquely determine the remaining elements. 

Eq. (2.26) implies that £ depends only on T 0 and is independent of the value of the 
remaining elements of S. The question therefore arises as to whether there is a particular 
choice of T 0 for which the resulting value of £ can be used as a strict criterion for the 
applicability of a singular perturbation analysis to Eqs. (2.20)-(2.22). The answer to this 
question is negative in general, because, in a given time interval, it is both the relative 
magnitudes of the three quantities: 

dE dy dh 
dt ’ dt ’ dt 


and the boundary conditions of interest that determine the validity of a singular 
perturbation analysis. Specifically, for an aircraft to exhibit the well-known two-time-scale 
behavior in a given time interval, it is necessary that for some choice of control the 
relations: 


d§ , , dy 
dt " dt 


(2.27) 
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(2.28) 


dE ,, dh 
¥ " dt 

be valid in that interval, and that the required change in E be sufficiently large to permit the 
boundary layer responses in h and y to reach their equilibrium values. Hence, it will be 
assumed that thJ boundary conditions are such that the resulting change in E is sufficiently 
large. Then, the conditions in Eqs. (2.27) and (2.28) further require that the net change in 
E during the boundary layer response is sufficiently small to permit approximating E as a 
constant (to zero order in e) in the boundary layer analysis. In addition, the aim here will be 
only to identify whether this two-time-scale property is a consequence of the inherent 
dynamics of the aircraft, and not whether it is a consequence of using a high gain control 
solution for L. Therefore, it will be assumed that the control L resulting from the boundary- 

layer analysis is of order one in Eq. (2.21). 

Under the above assumptions, there is a choice for T 0 for which the value of e can be 

used as a criterion for the existence of time intervals in which two-time-scale behavior is 
exhibited. If the choice of T 0 is such that dE/dt in Eq. (2.20) is at most of the same order of 
magnitude as £dy/dt and edh/dt in Eqs. (2.21) and (2.22), then, a value of £ sufficiently 
less than one indicates the possible existence of such intervals. By suitably choosing V 0 one 
can restrict V to be of order one. Then, by selecting the flight condition where the 
difference between thrust and drag, T-D, reaches a maximum, and by choosing T 0 to be 

equal to this difference: 


T„=(T-D).„ (2.29) 

one can guarantee that dE/dt is of order one, and both dy/dt and dh/dt are of order l/£. For 
this choice of T 0 , £ is given by: 


(T-D), 
1 mg 


(2.30) 


and is equal to the maximum longitudinal load factor of the vehicle. 

Eq. (2.30) actually represents an upper bound for £ (i.e. £<£,) since it is obtained by 
selecting the flight condition where the difference between thrust and drag, T-D, reaches a 
maximum The logical choice for V 0 is the speed at this flight condition. One can also adopt 
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the less conservative viewpoint of evaluating e along the energy climb path that represents 
the reduced solution (the solution obtained for e = 0). The value of e as a function of E can 
then be used as a measure to distinguish energy levels at which a singular perturbation 
analysis may be appropriate from other energy levels at which it may not be valid. 

It is important to emphasize here that having a high-gain control solution for lift, L, 
does not change the above conclusion, since high lift results in further time-scale separation 
between the flight path angle dynamics and the energy dynamics. It is precisely for this 
reason that the aim was only to identify whether the two-time-scale property is a 
consequence of the inherent dynamics, and not whether it is a consequence of using a high 
gain control solution for lift. It was done with the hope that this would lead to a conclusion 
independent of the performance index. One might also wish to exclude the situation 
wherein the open loop dynamics are not two-time-scale, but the closed-loop dynamics are 
two-time-scale. This would be the nonlinear counterpart to the so-called "cheap control 
problem" in linear quadratic optimization". 

Much can be anticipated from Eq. (2.30) for conventional aircraft without exact 
numerical evaluation. The quantity (T-D) max divided by mg is approximately equal to siny^ 
where y max is the maximum climb angle that can be maintained at a given energy level 
without loss of airspeed. It follows therefore that e,<l for all such aircraft types. For 
transport aircraft siny max is approximately 0.1, while for many fighter aircraft siny max is 
approximately 0.8 or less. This suggests that the forced singular perturbation analysis used 
in past studies of optimal aircraft trajectories is valid for most conventional subsonic and 
supersonic aircraft. 

A second upper bound for e can also be evaluated in terms of the quantities (T/mg)^ 
and (L/D)^ for a given aircraft. Since L is approximately equal to mg along the reduced 
solution corresponding to an energy climb path, it follows that: 

£<£2 (2.31) 

where is defined as: 


e 2 = 




0-'D). 


(2.32) 


Estimates of are given in Table 2. 1 . 
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Table 2.1 Estimation of e 2 based on Eg. (2.32) 


Parameter 

m. 

Transports 

Fighters 

(T/mg)^ 

0.25 

0.90 

(l/dl 

13 - 15 

4-7 

^2 

0.17-0.18 

0.65 - 0.76 


2.3 Hypersonic Regime 

Consider the flight of a hypersonic and possibly transatmospheric vehicle, viewed as 
a point mass, in a vertical plane over a spherical, non-rotating Earth, of gravitational 
strength p. The equations governing such flight can be reduced to a four-state model in E, 
m, y and radial distance from the center of the Earth, r. The equations are: 

dE = v ( t 1 t ~ D) (2.33) 

dt m 


f-'fcv-’O 


(2.34) 


dy _ f L ^ _ ( jxcosy ^ f Vcosy 
"dt" ~ V mV J~\ Vr 2 J + l r 


— = Vsiny 
dt 


(2.36) 


where T is now the maximum available thrust at a given speed and altitude. The control 
variables are L and T|, where 0 < r\ < 1 is introduced as a nondimensional throttling 

variable. E and V are now related by: 


E 



(2.37) 


Eqs. (2.35) and (2.36) now assume the form: 


e 


dy 

dt 



(2.38) 
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dr . 
e — = Vsinv 
dt 1 


(2.39) 


and again in order to avoid the artificial introduction of e, Eqs. (2.38) and (2.39) will serve 
only as a guide for the natural introduction of e. 

2.3.1 Nondimensional Form 


In order to put Eqs. (2.33) through (2.36) in nondimensional form one can now 
define the set of arbitrary positive quantities: 


and impose the restrictions that: 
tg has dimensions of time 
E 0 has dimensions of energy per unit mass 
nio has dimensions of mass 
r 0 has dimensions of length 
V 0 has dimensions of speed 
f 0 has dimensions of mass per unit time 
T 0 , D 0 , and L 0 have dimensions of force 
Using the elements of S to define the nondimensional quantities: 

+ t c E m r 

t = - ; E = — ; m = — ; r = - ; ' 


V = ^- 
V„ 


(2.40) 


(2.41) 


f=! 

fft 


T = — 


(2.42) 


Eqs. (2.33) through (2.36) can be put into the following nondimensional form: 

dE_ V(r 1 TT 0 -DD 0 ) r to v o >> 
dt m U 0 m 0 , 


(2.43) 


ik kv,,) 


(2.44) 
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dy 

'dt 



V 

A 


L t 1 
^0 l 0 

m 0 v o j 


\ 

cosy 

Ht 0 1 

+ 

r Vcosy^| 

(\ t ^ 
v 0 l 0 

l Vr 2 J 

V r 2 
v o l o7 


l r J 

r 0 J 


(2.45) 


£ 

dt 


Vo 
V r c ) 


Vsiny 


(2.46) 


The goal is again to put Eqs. (2.43) through (2.46) in the traditional singular perturbation 
form, thus, multiplying both sides of Eqs. (2.45) and (2.46) by (r 0 / V 0 ^ ) results in: 


V oW 


dy 

dt 


mV 


V l r ^ 

*0 

m 0 v oy 


^cosy^ 


Vr 2 


A 


V 0 2 r 0 y 


''^Vcosy^ 


r 


(2.47) 


r _£o_|£ = vsiny (2.48) 

Ivjdt 

Comparing the set of Eqs. (2.43), (2.44), (2.47) and (2.48) with the set (2.33), (2.34), 
(2.38) and (2.39) results in the following /ive conditions on the elements of set S: 


H 

O 

II 

O 

o 

(2.49) 

T t V 
1 o L o v 0 _ J 

E 0 m 0 

(2.50) 

f 0 *0 _ j 

(2.51) 

m 0 


L ° r o =1 
m 0 V 2 

(2.52) 

f =1 

(2.53) 

V 0 2 r 0 



By defining £ as: 



Eqs. (2.43), (2.44), (2.47) and (2.48) assume the traditional singular perturbation form: 


dE _ V(r)T-D) 
dt m 

(2.55) 


(2.56) 

dyf L 1 fcosy^j fVcosy^ 

(2.57) 

£ dt [mvj [vr 2 jA r J 

dr .. . 

£— = Vsiny 
dt 1 

(2.58) 


2.3.2 Specifying a 


Particular Nondimensional Form 


For the hypersonic case, only five conditions on the nine elements of set S are needed 
in order to put the equations of motion in the traditional singular perturbation form. Thus, 
one can specify four of the elements of S to fit one's convenience and then determine the 
remaining five using Eqs. (2.49)-(2.53). 

Again, in order to maintain the relationship in Eq. (2.37) in the transformed variables, 
a sixth condition is introduced: 



which together with Eq. (2.53) gives: 



2 r 


(2.59) 


(2.60) 


If one thinks of r 0 as a radial distance, then Eq. (2.53) restricts V 0 to be the circular orbital 
speed at r 0 . Similarly, Eq. (2.52) restricts L 0 to be the centripetal force that a point mass 
would experience in a circular orbit at r 0 . Using Eqs. (2.50), (2.53), and (2.59) one 
obtains: 

T r 2 

e = -2-2- (2.61) 

pm 0 
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Hence, by picking three among the nine elements of set S arbitrarily, the six conditions 
Eqs. (2.49)-(2.53), and (2.59) uniquely determine the remaining elements. 

The question arises again as to whether there is a particular choice for these three 
elements for which the resulting value of e can be used as a criterion for the applicability of 
a singular perturbation analysis to Eqs. (2.55)-(2.58). The right-hand-sides of Eqs. (2.55) 
and (2.58) can bi made of the same order of magnitude by choosing the ratio T 0 /rr^ as: 


T 0 f TlT-D ^ 

m 0 V m / max 


(2.62) 


Choosing r 0 as the sea level radius r SL , r and V are of order one. Also, for these choices of 
T 0 /mo and r 0 , dE/dt is of order one, and both dy/dt and dr/dt are of order 1/e. By choosing 
f 0 as the value of f at the flight condition where the ratio (r|T- D)/m is a maximum, dm/dt 
can also be made of order one. With the above choices of T 0 /mo, r 0 , and f 0 , e becomes: 


e = 


(x 1 \ 

£sl 

V J 


( riT-D 
V m 


(2.63) 


The right-hand-side of Eq. (2.63) is the maximum longitudinal load factor of the 
vehicle in units of sea-level g's, and actually represents an upper bound for e since it is 
obtained by selecting the flight condition where (qT- D)/m reaches a maximum. One can 
again adopt the less conservative viewpoint of evaluating e along the energy climb path that 
results from the reduced solution (the solution obtained with £ = 0). The value of £ as a 
function of E can then be used as a measure to distinguish energy levels where a singular 
perturbation analysis may be appropriate, from other levels where it may not be valid. 

A hypersonic flight vehicle employing an airbreathing propulsion system and sized 
for acceleration to orbital velocity necessarily employs a multimode propulsion system. An 
example might include turbojet, ramjet, scramjet, and rocket modes. Each mode of 
propulsion when considered in its operating Mach regime can be characterized by a 
corresponding £. Available models of this vehicle type exhibit large values of excess thrust 
at low hypersonic Mach numbers. In fact, Eq. (2.63) will produce an £ that is greater than 
one over such flight phases. Experience with hypersonic vehicle dynamics reported in Ref. 
12 indeed suggests that the assumed time scale separation is not valid in these phases. 
However, as will be seen in the next section, over most of the trajectory corresponding to 
the reduced solution Eq. (2.63) results in an £ that is less than one, just as in the flat Earth, 
subsonic-supersonic case. 
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2.4 Numerical Validation 


It was shown in the preceding sections that for aircraft energy climbs that take place 
in a vertical plane, the singular perturbation parameter e can logically be identified as the 
maximum longitudinal load factor of the vehicle, measured in units of sea-level g's. In 
order to further explore the implications of this result, numerical evaluations of e will be 
presented in this section for several types of vehicles. As stated earlier, it will be assumed 
that the required change in specific energy is sufficiently large to permit the boundary layer 
responses in altitude and flight path angle to reach their equilibrium values. 

For a given aircraft, it may be sufficient to evaluate a single upper bound for e, valid 
for the entire envelope, in order to provide a hint regarding possible two-time-scale 
behavior. If the resulting value of this upper bound is sufficiently less than one, and if the 
boundary conditions are appropriate, then two-time-scale behavior is implied for any 
energy climb that the aircraft is allowed to perform. If however the resulting value of the 
upper bound turns out to be greater than one, then no conclusion can be drawn. The way to 
proceed in the latter case would be to evaluate a less conservative (smaller) upper bound for 
£ and apply the same reasoning. Unfortunately, the less conservative the upper bound, the 
more computation one has to perform in order to evaluate it. In particular, e in Eq. (2.30) 
or Eq. (2.63) can be evaluated as a function of energy E, using all the assumptions made in 
the evaluation of reduced solutions in aircraft energy climbs (Y=0, L=mg etc.). By 

evaluating in this sense, and at each energy level the absolute maximum value of the 
longitudinal load factor one obtains a curve C on the e-E plane. The interesting properties 
of this curve are that for a given aircraft it need only be constructed once and that it lies 
above any other curve that may be evaluated similarly, but along the actual reduced solution 
corresponding to the specific problem of interest. In other words, points on curve C 
represent upper bounds for £ at the corresponding energy levels. The portions therefore of 
curve C where £ is sufficiently less than one immediately show the energy levels where 
two-time-scale behavior (boundary-layer transitions along constant E) can be expected. If 
there are any portions of curve C where £ is greater than one, then no conclusions can be 
drawn as to the possible two-time-scale behavior at the corresponding energy levels. In the 
latter case one has again to evaluate a still less conservative upper bound for £ at these 
energy levels. Such less and less conservative upper bounds would of course eventually 
lead to the maximum value of the longitudinal load factor evaluated as a function of E along 
the reduced solution corresponding to a specific problem. 
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If one is interested in the possible two-time-scale behavior of a vehicle along a 
particular trajectory (corresponding to a specific problem) then, the least conservative upper 
bound for e would be the maximum longitudinal load factor encountered along that (exact) 
trajectory. Such a test for time-scale separability would require computation of the (exact) 
trajectory first and would not be very helpful. Hence, it makes much more sense to start 
with a more conservative (greater) upper bound and to proceed with less and less 
conservative upper bounds. 

Numerical evaluations of £ are presented in Figs. 2.1 through 2.8 for four types of 
vehicles. For each type there is a plot showing the variation of the maximum longitudinal 
load factor of the vehicle with energy E, and one or more plots showing the variation of the 
longitudinal load factor with E along the reduced solution corresponding to a specific 
optimization problem. The odd-numbered figures are the energy-climb figures for the four 
vehicles. They distort the fact that the altitude profiles contain jumps, because they only 
needed to be constructed roughly, since they only served for the evaluation of £ as a 
function of the energy E, given in the even-numbered figures. 

Figs. 2.1 and 2.2 show the results for an F-8 fighter 13 . The two optimization 
problems considered for this case were minimum time to a specified energy and minimum 
time to a specified downrange position. The reduced solutions corresponding to these 
problems are obtained by maximizing (with respect to V) at each energy level the quantities 
(T-D)V for the former and [(T-D)V]/(V 0 -V) for the latter. In this case V 0 is the maximum 
possible cruising speed of the aircraft and D is calculated at L=mg. Fig. 2.1 shows the 
actual paths in the envelope corresponding to these reduced solutions and to the maximum 
longitudinal load factor of F-8. Fig. 2.2 shows the results for E evaluated along these 
climb paths. Since the maximum longitudinal load factor of F-8 remains below 1.0 in Fig. 
2.2, it is reasonable to assume that for any optimization problem, if the required energy 
gain is sufficient, the transitions to the reduced solution will take place at nearly constant E, 
exhibiting the well known boundary layer structure. 

Figs. 2.3 and 2.4 show similar results for an F-15 fighter 14 . Again, the problems of 
minim um time to a specified energy and minimum time to a specified downrange position 
were considered. A maximum dynamic pressure constraint of 1500 lbf per square feet is 
imposed on the climb paths in each case. Due to the large thrust to weight ratio of F-15, the 
£ levels in Fig. 2.4 are much higher than those of F-8 (i.e. in comparison with Fig. 2.2). In 
particular, there is a small region at low energy where £ exceeds one, implying that time- 
scale separation at these energy levels may not be appropriate for these problems. 
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Figs. 2.5 and 2.6 show the results for a conventional transport 15 . In this case, 
however, the two optimization problems considered were minimum fuel to a specified 
energy and minimum fuel to a specified downrange position. The reduced solutions 
corresponding to these problems are obtained by maximizing (with respect to V and T|) at 
each energy level the quantities [(T-D)V]/f for the former and [(T-D)V]/(fY 0 -f 0 V) for the 
latter. Here V 0 represents the most fuel efficient cruising speed of the aircraft and f 0 is the 
fuel consumption rate at this cruising flight condition 7 . The magnitude of e in Fig. 2.6 
remains small in comparison to that in Figs. 2.2 and 2.4. Two-time-scale behavior for the 
entire transport aircraft envelope is therefore suggested. 

Finally, Figs. 2.7 and 2.8 show the results for a hypersonic vehicle model, used by 
NASA and termed the "Langley Accelerator” 16 . The only optimization problem considered 
in this case was minimum fuel to a specified energy, the reduced solution corresponding to 
which is obtained by maximizing the quantity [(TiT-D)V]/(mf) at each energy level (mass is 
not constant in this case). A maximum dynamic pressure constraint of 2000 lbf per square 
feet is imposed on the climb paths for this case. This particular vehicle model employs a 
multimode propulsion system, sized for acceleration to orbital velocity, and consisting of 
turbojet, ramjet, scramjet, and rocket modes. Optimal switching between propulsion modes 
was not included in the generation of Figs. 2.7 and 2.8. Instead, within allowable Mach 
regimes the cycle that maximizes [0iT-D)V]/(mf) and 0iT-D)/m was selected. The points of 
cycle transitions are shown in the figures. In Fig. 2.8, e is plotted against the speeds at 
which the constant energy contours intersect the zero altitude axis. The reason for this is 
that E is negative in this case so that its size is no longer intuitively obvious. Thus, sea- 
level speed was used as the abscissa, because at hypersonic speeds practically all the 
energy of the vehicle corresponds to kinetic energy. The calculated value of e will likely be 
reduced if a practical method for cycle transition is employed. Note that as the energy 
increases the boundary of the envelope is approached and e goes to zero. This is a basic 
characteristic of all aircraft (see also Figs. 2.2, 2.4, and 2.6), suggesting that transitions to 
the reduced solution can be treated as boundary layers with relatively greater success at 
higher energy levels. The physical explanation for this rests bn understanding the behavior 
of the difference between thrust and drag. At low energy levels both the speed and altitude 
are low, resulting in high thrust and low drag, so that the difference between thrust and 
drag is high. This large amount of excess power can be used to effect a change in energy 
even during a transition. However, at higher energy levels, either speed or altitude or both 
are high, and the corresponding difference between thrust and drag is low. Thus, 
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transitions to the reduced solution at higher energy levels occur for the most part by 
interchanging kinetic for potential energy (or vice versa), with the total energy remaining 
more nearly constant. The expression found for e (see Eqs. (2.30) and (2.63)) captures and 
quantifies this effect in a straightforward manner. 


2.5 Conclusions 

A systematic procedure has been introduced to identify a singular perturbation 
parameter in the differential equations of motion for both the conventional (subsonic- 
supersonic, flat Earth) and the transatmospheric (hypersonic, spherical Earth) flight 
regimes. The procedure uses a set of arbitrary scaling constants to nondimensionalize all 
the variables of interest. Nondimensionalization alone is not sufficient to clearly identify if a 
system will exhibit two-time-scale behavior. However, there is a useful choice of the 
scaling constants that leads to the conclusion that two-time-scale behavior can be expected 
when the maximum longitudinal load factor is sufficiently less than one. The important 
point here is that this statement is valid independent of the performance index that is being 
optimized. 

This explains the past successes in singular perturbation treatments of aircraft energy 
climbs, despite an inability to explicitly identify a perturbation parameter. These 
observations also apply to the family of future hypersonic vehicles. If such a vehicle is 
employed as a passenger transport, its acceleration will necessarily be constrained in the 
interest of human comfort. To constrain the maximum longitudinal load factor of such a 
vehicle to be sufficiently less than one would imply two-time-scale behavior for any type of 
energy climb that such a vehicle would be allowed to perform. 
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Chapter III 


State-Constrained Energy-State Modeling 

3.1 Introduction 

State-variable inequality constraints are commonly encountered in the study of 
dynamic systems. The study of rigid body aircraft dynamics and control is certainly no 
exception. For instance, a maximum allowable value of dynamic pressure is usually 
prescribed for aircraft with supersonic capability. This limit is required to ensure that the 
vehicle's structural integrity is maintained. Given a typical state-space description of the 
vehicle dynamics, this limit constitutes an inequality constraint on the vehicle state. Such 
dynamic pressure bounds are commonly encountered during fuel-optimal climb for 
supersonic transports 19 , for rocket powered launch vehicles such as the U.S. space 
shuttle 20 , and for single-stage-to-orbit air-breathing launch vehicles 21 . 

State-variable inequality constraints have been studied extensively by researchers in 
the field of optimal control. First-order necessary conditions for optimality when general 
functions of the state are constrained have been obtained 22 - 24 . However, the direct 
construction of solutions via this set of conditions proves difficult. Moreover, the 
controllers derived from such traditional formulations of the problem suffer from serious 
practical flaws. They typically tend to track the hypersurface representing a state-constraint 
boundary in a finite time, which makes the traditional asymptotic boundary layer theory 
non-applicable (see below), and they break down whenever a disturbance causes the 
system to violate an active state-constraint. Accordingly, most practitioners seeking an 
open-loop control solution rely on direct approaches to optimization that employ penalty 
functions for satisfaction of state-variable inequality constraints 25 . As a rule however, 
algorithms employing such methods are computationally intense and slow to converge. 
Consequendy, they are not well-suited for real-time implementation. 

From a singular perturbations point of view 26 - 27 , in the absence of a state-variable 
inequality constraint (i.e. when no constraint is active), the initial boundary-layer solution 
for the class of systems being considered is an infinite-time process. A solution is sought 
which asymptotically approaches the reduced solution (solution with e=0). However, when 
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a state-constraint is active in the reduced solution, the boundary-layer problem can be of 
finite-time in the stretched time variable 28 29 . Thus, traditional techniques concerning the 
asymptotic stability of the boundary-layer system are not applicable, and cannot be used to 
construct an approximate boundary-layer solution. The presence of an active state-variable 
inequality constraint also introduces the possibility of discontinuous costate variables at the 
juncture betweeli constrained and unconstrained arcs. A Valentine transformation can be 
used to convert the constrained problem to an equivalent unconstrained problem of 
increased dimension 30 . Smoothness is regained in the process, but at the expense of 
introducing a singular arc along the state-constraint boundary 30 31 and to little or no 
advantage when seeking a solution for real-time implementation. 

As an alternative, this chapter proposes a complete reformulation of optimal control 
problems involving active state-variable inequality constraints. Since in practice it is always 
the asymptotic controllers that have the most desirable properties, maybe the optimization in 
such problems should be carried out not over the class of all controllers, but only over the 
class of asymptotic controllers that track a given active state-constraint boundary. It is 
shown in Section 3.2 that a transformation technique can be used to isolate and describe 
completely this class of asymptotic controllers. If a minimum over the class of all 
controllers exists, then the reformulated problem is guaranteed to have an infimum. The 
results in Section 3.3 suggest however that this infimum for the reformulated problem 
corresponds to a finite-time controller and is not achieved over the class of all asymptotic 
controllers. The situation is somewhat reminiscent of H-infinity control theory for linear 
systems in which one seeks a strictly proper, stabilizing controller to minimize the H- 
infinity norm of a closed-loop transfer function. The minimum of this norm over all strictly 
proper, stabilizing controllers does not exist. Its infimum, however, does exist and 
corresponds to a proper controller. Thus, just as in H-infinity theory, the question arises 
naturally in the present case as to how one can find an asymptotic controller that somehow 
approximates this infimum. Although there are no general answers yet, a procedure is 
presented in Section 3.3 that does supply one with insight at least for a simple example. 
The procedure leads to a controller that is approximately optimal, asymptotic, but still one- 
sided. 

Finally, if the state-constraint is regarded as a soft constraint, then an example shows 
that there may exist controllers that are asymptotic , two-sided, and result in exactly the 
same optimal value of the performance index corresponding to the original problem, that is, 
they are practically optimal, but at the expense of violating the state-constraint. Such 
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controllers, however, do not correspond to stationary solutions of the optimization 
problem. Thus, at the present, a systematic procedure for finding them does not exist and 
remains a topic for future research. Preliminary results on the topic covered in this chapter 
can be found in Ref. 32. 


3.2 Construction of Arbitrary Nonlinear Feedback Control Laws for 
a Dynamical System, that Track a Given Hypersurface 

Consider the dynamical system: 

x(t 0 ) = x 0 (3.1) 

y(to) = yo (3.2) 

where x, f are vectors of the same dimension, and y, g and u are scalars. It will be of 
interest to describe the set C of all piecewise continuous (in time) control laws u(x(t),y(t),t) 
that track a given hypersurface in the state space of the above system, given by the scalar 
equation: 

S(x,y) = 0 (3.3) 

that is, if u=u(x,y,t) is a specific control law belonging to C, then the system of Eqs. (3. 1), 
(3.2), driven by u(x,y,t) for t > ^ (and assuming that S(x 0 ,y 0 ) is not zero) will eventually 
reach the hypersurface given in Eq. (3.3) and stay on it thereafter. A control law can drive 
the system onto the hypersurface either in finite time or asymptotically. Accordingly, the set 
C is the union of two disjoint sets F and A. The set F contains all control laws that track the 
hypersurface in finite time. The set A contains all control laws that track the hypersurface 
asymptotically . The purpose in this section is to give a complete description of the sets F 
and A. 

Let Z be the set of all piecewise differentiable, scalar functions of the real variable a, 
defined and invertible for all a in [0,1], and satisfying the boundary conditions: 

z (0) = 0 ; z(l) = — S(x 0 ,y 0 ) (3.4) 


— = f(x,y,u) ; 
^ = g(x,y,u) ; 
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Assume that the range of S is contained in the range of z for all z in Z. Then, consider a as 
an independent variable and make the transformation from y to a: 

z(a) + S(x,y) = 0 (3.5) 

Differentiating Eq. (3.5) with respect to time and using Eqs. (3.1) and (3.2) results in: 


dz da | 
da dt 


dS cl x dS / \ a 

— f(x,y,u) + — g(x,y,u) = 0 


(3.6) 


Now, let da/dt play the role of a new control, P, by defining: 



(3.7) 


then, Eq. (3.6) becomes: 


dz D 9S ./ v 3S , \ n 

— p +— f(x,y,u) + — g(x,y,u) = 0 
da dx dy 


(3.8) 


Assume now that the hypersurface given by Eq. (3.3) is first order in u, that is, the total 
time derivative of S(x,y) is explicitly dependent on u. If this is not the case, the results of 
this section can be generalized to the case where the hypersurface given by Eq. (3.3) is 
higher order in u. Also, assume that Eq. (3.5) is invertible in y, and that Eq. (3.8) is 
invertible in u. These assumptions result in the two equations: 

y = h(x,z) (3.9) 


u = 


k| x,z,~— P 
da 


(3.10) 


for y and u respectively. The system of Eqs. (3.1) and (3.2) has now been transformed 
through Eqs. (3.5)-(3.10) to the equivalent system: 


dx 

dT 


f x, h(x, z(a)), kf x,z(a),^P 


(3.11) 
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(3.12) 



with initial conditions x(t 0 )=x 0 and a(to)=l. 

3.2.1 Finite-Time Tracking 

Theorem 3.1. For any given finite time t f and any function z(a) in Z, there exists 
a control law u=u(x,y,t) in F that drives the system of Eqs. (3.1) and (3.2) from its initial 
state at t=t„ onto the hypersurface given by Eq. (3.3) at t=tf and keeps it on the hypersurface 
thereafter. 

Proof 3.1. Let t f be a finite time and z(a) be a function in Z. One can use z(a) to 

obtain the equivalent transformed system of Eqs. (3.11), (3.12). Then, one can use as the 
control p the function: 


P = - 


t f 



forto<t<t f (3.13) 


p = o 


leading to the time variation for a: 


a = 


tft 
If — 


for tf < t (3.14) 


for t,, < t < t f (3.15) 


a = 0 for tf < t (3.16) 

Thus, a is driven from 1 to 0, in tf, < t < tf and stays at zero for t f < t. Accordingly, 
due to the boundary conditions (Eq. (3.4)) on the function z(a), the original system of 
Eqs. (3.1), (3.2), is driven from its initial state at t=tf, onto the hypersurface given by Eq. 
(3.3) in ^ < t < t f and stays on the hypersurface for tf < t. The feedback controller u(x,y,t) 
that will perform this task for the system of Eqs. (3.1), (3.2) can be found from Eq. 
(3.10). P is given by Eqs. (3.13), (3.14), z is equal to -S(x,y) from Eq. (3.5), and since 
z(a) is invertible on [0,1], dz/da can be expressed as a function of z and therefore of 
-S(x,y) (Q.E.D.). The procedure is best illustrated through an example. 
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Example 1 

Consider the system: 


dx 2 

* =y - u ; 


dy 

j = u ; 

dt 


o 

X 

II 

o 

(3.17) 

*-»■ 

o 

II 

o 

(3.18) 


and assume that it is desired to track the straight line: 

y - 1 = 0 

in a specified final time tf. As the function z(a) one can select: 

z ( a )= fe) n(1+c0 


(3.19) 


(3.20) 


Note that z(0)=0, and z(l)=l-y 0 as required by Eq. (3.4). From Eq. (3.20) one is led to the 
transformation: 


\ 2° y + <x ) + y-1 = ® 


(3.21) 


Differentiating Eq. (3.21) with respect to time, one obtains: 

(!z£lYJL ) +u=0 

l In 2 Al + aJ 

resulting in the transformed system: 


(3.22) 


^ = 1 -fiz2i ) to(1 +a) _[riz 5 i' r_e_ 

dt l In 2 ) \ ln2 \l + a 


(3.23) 


— = 8 (3.24) 

dt 

with initial conditions x(t 0 )=x 0 and a(to>=l. Now, if one uses the control (5 given in Eq. 
(3.13), then, from Eq. (3.22) one obtains: 
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(3.25) 


Izlo 

ln2(t f -t 0 )' 


U = — — ; . exp< -In 2 


i-y 

l ~y 0 . 


This open-loop control can be shown to drive the system of Eqs. (3.17), (3.18) from its 
original state at t=t„ onto the line y=l at t=t f . To obtain a closed-loop (feedback) controller 
one simply replaces y 0 by y and ^ by t in Eq. (3.25) to obtain: 


i-y 

ln4(t f -t) 


(3.26) 


which again drives the system to y=l at t=t f . 

Theorem 3.2. Let u=u(x,y,t) be any control law in F that drives the system of 
Eqs. (3.1), (3.2) from its initial state at t=to onto the hypersurface Eq. (3.3) at t=t f and 
keeps it on the hypersurface thereafter. Then, there exists a function z(a) in Z, such that, 
when the system of Eqs. (3.1), (3.2) is transformed using z(a) to the system given by 
Eqs. (3.1 1), (3.12), the control p for the transformed system corresponding to u(x,y,t) is 
given by Eqs. (3.13), (3.14). 

Proof 3.2. The function z(a) for 0 < a < 1 is found by solving the system of 
differential equations: 


dx 

da 


=-( t f- t o) f ( x »y^ u ( x .y.O) 


(3.27) 



(3.28) 


subject to the boundary conditions: 


x(a = l) = x 0 ; z(a = 1) = -S(x 0 ,y 0 ) (3.29) 

where y in Eqs. (3.27), (3.28) is a function of x and z through Eq. (3.9) and t in u(x,y,t) is 
a function of a through Eq. (3.15), that is, 

t = t f -a(t f -t 0 ) (3.30) 
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Note that the function z(a) found from the solution of Eqs. (3.27) through (3.29) 
satisfies Eq. (3.5) for all t > ^ (a and t are related by Eq. (3.15)). Thus, it also satisfies the 
boundary condition z(0)=0, since at t = t f , a and S(x,y) are both equal to zero (Q.E.D.). 

Theorems 3.1 and 3.2 indicate that once the control p for the transformed system of 
Eqs. (3.11), (3.12) is fixed as in Eqs. (3.13), (3.14), there is a complete correspondence 
between functions in Z and piecewise continuous control laws u(x,y,t) in F that track the 
hypersurface given by Eq. (3.3) in a finite time q. That is, for every element of Z there 
exists an element of F and, more importantly, for every element of F there exists an element 
of Z. The correspondence therefore established between F and Z through the selection of (3 

as in Eqs. (3.13), (3.14) is onto. 

3.2.2 Asymptotic Tracking 

Theorem 3.3. For any function z(a) in Z, there exists a control law u=u(x,y,t) in 
A that drives the system of Eqs. (3.1), (3.2) from its initial state at t=t 0 onto the 

hypersurface given by Eq. (3.3) asymptotically. 

Proof 3.3. Let z(a) be a function in Z. One can use z(a) to obtain the equivalent 

transformed system of Eqs. (3.1 1), (3.12). Then, using the control P given by: 


P = -a (3.31) 

leads to the exponential time variation for a: 

a = e" (3.32) 

Thus, cc is driven exponentially from 1 to 0. Accordingly, due to the boundary 
conditions (Eq. (3.4)) on the function z(a), the original system of Eqs. (3.1), (3.2), is 
driven asymptotically from its initial state at t^ onto the hypersurface given by Eq. (3.3). 
The feedback controller u(x,y,t) that will perform this task for the system of Eqs. (3.1), 
(3.2) can be found from Eq. (3.10). P is given by Eq. (3.31), z;s equal to -S(x,y) from 
Eq. (3.5) and since z(a) is invertible on [0,1], dz/da can be expressed as a function of z 
and therefore of -S(x,y) (Q.E.D.). Again, the procedure can best be illustrated through an 
example. 
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Example 2 

Consider again the system: 


dx , 

¥ = y ‘ U ; 


dy 

dt 


= u ; 


o 

X 

II 

o 

(3.33) 

© 

II 

'C 

o 

(3.34) 

the straight line: 


) 

(3.35) 


As the function z(a) one can again select: 


z(a)= (i^) ln(1+a) 


(3.36) 


satisfying z(0)=0, and z(l)=l-y 0 as required by Eq. (3.4). From Eq. (3.36) one is led to 
the transformation: 


(^2 L ) ln(1+ “ ) + y “ 1 = 0 


(3.37) 


Differentiating Eq. (3.37) with respect to time one obtains: 


Izlo 

In 2 


1 + a 


+ u = 0 


(3.38) 


resulting in the transformed system: 



flzlo] 

In 2 J 


ln(l + a)- 



(3.39) 



(3.40) 
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with initial conditions x(t 0 )=x 0 and a(t 0 )=l. If now one uses the control p given in Eq. 
(3.31), then, from Eq. (3.38) one obtains: 


u = 


e k 0-y)_i 

ke k(, - y) 


(3.41) 


where k is defined as: 


k = 


In 2 

i-y 0 


(3.42) 


This open loop control can be shown to drive the system of Eqs. (3.33), (3.34) from its 
original state at t=to asymptotically onto the line y=l. To obtain a closed-loop (feedback) 
controller one simply replaces y 0 by y in Eq. (3.41), to obtain: 


u = 


i-y 

In 4 


(3.43) 


which again drives the system of Eqs. (3.33), (3.34) asymptotically toward y=l. Note that 
as expected, in the asymptotic case u does not depend explicitly on t. 

Theorem 3.4. Let u=u(x,y,t) be any control law that drives the system of Eqs. 
(3.1), (3.2) from its initial state at t=to onto the hypersurface given by Eq. (3.3) 
asymptotically. Then, there exists a function z(a) in Z, such that, when the system of Eqs. 
(3.1), (3.2) is transformed using z(a) to the system given by Eqs. (3.11), (3.12), the 
control P for the transformed system corresponding to u(x,y,t) is given by Eq. (3.31). 

Proof 3.4. The function z(a) for 0 < a < 1 is found by solving the system of 
differential equations: 

a^ = -f(x, y ,u(x, y ,t)) (3.44) 

a ^] = | f ( x , y , u(x , y , t )) + | g ( x , y , u (x, y .,)) (3.45) 

subject to the boundary conditions: 

x(a = 1) = x 0 ; Z (a = l) = -S(x 0 ,y 0 ) (3.46) 
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where y in Eqs. (3.44), (3.45) is a function of x and z through Eq. (3.9), and t in u(x,y,t) 
is a function of a through Eq. (3.32), that is. 


t = In 



(3.47) 


Note that the function z(a) found from the solution of Eqs. (3.44) through (3.46) 
satisfies Eq. (3.5) for all t > tg (oc and t are related by Eq. (3.32)). Thus, it also satisfies the 
boundary condition z(0)=0, since as t approaches infinity a and S(x,y) both tend to zero 
(Q.E.D.). 

Theorems 3.3 and 3.4 indicate that once the control P for the transformed system of 
Eqs. (3.11), (3.12) is fixed as in Eq. (3.31), there is a complete correspondence between 
functions in Z and piecewise continuous control laws u(x,y,t) in A that track the 
hypersurface given by Eq. (3.3) asymptotically. That is, for every element of Z there exists 
an element of A and, more importantly, for every element of A there exists an element of Z. 
The correspondence therefore established between A and Z through the selection of P as in 
Eq. (3.31) is onto. 


3.3 Significance for Optimal Control Problems Involving Active 
State- Variable Inequality Constraints 

The ideas presented in Section 3.2 can be appropriately utilized in the study of 
optimal control problems involving active state-variable inequality constraints. A common 
feature of such problems is that when a portion of the optimal trajectory rides the 
hypersurface representing a state constraint boundary, the optimal transition to this 
hypersurface from an initial point that does not lie on it occurs in finite time. Consequently, 
the corresponding optimal feedback controllers for such problems, that can be obtained by 
well-known analytical or numerical methods, suffer from two flaws that tend to eliminate 
almost completely their practical usefulness. First, such feedback controllers are finite-time, 
meaning that, for two-time-scale systems traditional asymptotic boundary layer theory is 
not applicable. Second, if a disturbance throws the system instantaneously toward the 
prohibited side of the hypersurface, the feedback scheme breaks down and there is no 
"optimal" way of returning to the hypersurface 28 . 

Asymptotic controllers on the other hand, capable of tracking from both sides a 
hypersurface representing a state constraint boundary, presumably won't suffer from either 
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of the above two flaws. This observation, and the ideas presented in Section 3.2 suggest 
that once one knows that a portion of the optimal trajectory for a problem rides such a 
hypersurface, one can change one's point of view and try to optimize the system over all 
asymptotic controllers capable of tracking that hypersurface from both sides. 

Consider therefore that one wishes to minimize the performance index: 

>f 

J = 0(x(t f ),y(t f )) + J L(x,y,u)dt (3.48) 

l 0 

subject to the dynamical equations: 

dx 

— = f(x,y,u) (3.49) 

at 

^7 = g(x,y,u) (3.50) 

dt 

the boundary conditions: 

x(t 0 ) = x 0 ; y(t 0 ) = y 0 ; t 0 fixed ; t f free (3.51) 

and the state-variable inequality constraint: 

S(x,y) < 0 (3.52) 

As in Section 3.2, x, f are vectors of the same dimension, and L, S, y, g and u are 
scalars. Again, it will be assumed that: 

S(x 0 ,y 0 )<0 (3.53) 

and that the optimal trajectory reaches the hypersurface: 

S(x,y) = 0 (3.54) 

at a finite time t=t, and stays on it for t > t,. Thus, in order to avoid the problems with the 
finite-time controllers mentioned above, one would now like to optimize J over the class of 
all asymptotic controllers, which has already been denoted by A, capable of tracking the 
hypersurface given by Eq. (3.54) from both sides. 
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Based on the above observations, a new way of looking into the problem can now be 
formulated in two steps: In the first step (Section 3.3.1) the optimization is carried out over 
all asymptotic controllers, but the state constraint is retained. This however would lead at 
best to an asymptotic controller that is approximately optimal, does not violate the state 
constraint, but is one-sided, that is, it breaks down on the prohibited side of the constraint 
boundary. Since this one-sidedness would tend to eliminate the practical usefulness of such 
a controller, in the second step (Section 3.3.2) the state constraint is discarded, and by 
means of a simple example it is shown that it is possible to find an asymptotic, two-sided 
controller that although (unfortunately) violates the state constraint, it results in exactly the 
same optimal value of the performance index, corresponding to the optimal, finite-time, 
one-sided controller for the original problem (that does not violate the constraint). 

3.3.1 First Reformulation: Optimization over all 
Asymptotic, One-Sided Controllers 

Let z be an arbitrary function in the set of functions Z defined in Section 3.2. The 
transformation reads: 


z(oc) + S(x,y) = 0 (3.55) 

with the function z(a) subject to the boundary conditions: 

z (0) = 0 ; z(l) = -S(x 0 ,y 0 ) (3.56) 

This leads as in Section 3.2 (see Eqs. (3.5) through (3.12)) to: 

~P + 3 ~f( x .y» u ) + T-g(x,y, u ) = 0 (3.57) 

da dx dy 

and to the transformed system: 


y = h(x,z) (3.58) 

u = k|x,z,-^pj (3.59) 
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^ = f^x,h(x,z(a)),k|x,z(a),^pjj 


(3.60) 


— = p (3.61) 

dt 

with initial conditions x(t 0 )=x 0 and a(t 0 )=l. The performance index to be minimized 
assumes the form: 

t f 

J = o|x(t f ), h(x(t f ), z(a(t f )) j j + j L^x, h(x,z), k^x, z >^ P j j dt (3.62) 

‘o 

As seen in Section 3.2 (Theorems 3.3 and 3.4), once the control P for the 
transformed system of Eqs. (3.60), (3.61) is fixed at P=-a, there is a complete 
correspondence between functions z(a) in Z and piecewise continuous control laws 
u(x,y,t) in A that track the hypersurface given by Eq. (3.3) asymptotically. That is, for 
every element of Z there exists an element of A and, more importantly, for every element of 
A there exists an element of Z. Therefore, to find the "best" asymptotic controller, it is 
natural to fix the control P at: 

p = -a (3.63) 

and to then try to determine the "best" function z(a) in Z that minimizes J. This leads 
directly to the off-line optimization problem: 

Minimize: J = 0(x(0),h(x(0),0)) + J*^— ^L(x,h(x,z),k(x,z,— av))da (3.64) 

o 


subject to: 



(3.65) 


- 7 — = v (3.66) 

da 
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-z(a) < 0 


(3.67) 


and to the boundary conditions x(l)=x 0 , z(0)=0, and z(l)=-S(x 0 ,y 0 ). Note that the state 
constraint is still being retained through Eq. (3.67). This point is important, because, if Eq. 
(3.67) is not imposed, then the problem presumably will not have a solution. One can see 
from Eq. (3.64) that for every function z(a) in Z that satisfies Eq. (3.67) there corresponds 
a number J,(z). Thus, one can define the set of real numbers: 

J = { J,(z): z is in Z, and inequality (3.67) is satisfied} (3.68) 

and state the off-line optimization problem, given by Eqs. (3.64) through (3.67), by asking 
for the minimum of J. Although it is not known at the present if and exactly when such a 
minimum exists, it is possible to state the following theorem that provides one with an 
important partial answer: 

Theorem 3.5. If the original problem given by Eqs. (3.48) through (3.52), 
involving an active state variable inequality constraint, has a minimum, then the off-line 
optimization problem given by Eqs. (3.64) through (3.67) has an infimum. 

Proof 4.5. Assume that the problem given by Eqs. (3.48) through (3.52) has a 
minimum, which will be denoted by J min . This immediately implies that the set J is bounded 
below by J min , that is, J min < J,(z) for all z in Z. Since any set of real numbers that is 
bounded below has an infimum, J has an infimum (Q.E.D.). 

Example 3 

To illustrate the above idea one can apply it to the following problem: 

Minimize: J = J(l-y + u 2 )dt (3.69) 

o 

subject to: 

^ = u ; y(0) = 0 (3.70) 

dt 

y-l<0 (3.71) 

The solution to this problem for 0 < t < 2 can be shown to be 28 : 


35 



(3.72) 


r t 

y = i-t ; u hi 

7 4 2 

leading to the finite-time controller: 

u = (1 - y) 1/2 (3.73) 

Although this controller is optimal, it is clearly useless from a practical point of view since 
it breaks down when y exceeds 1. It cannot be used to track the line y=l in the presence of 
two-sided perturbations about y=l. At t = 2, y reaches the value y = 1. For t > 2, y and u 
stay constant at 1 and 0 respectively and there is no further contribution to the performance 
index J. In order to optimize J over all asymptotic controllers that track y = 1 and satisfy 
Eq. (3.71), one can now use the transformation: 

z(tx) + y-l = 0 (3.74) 



(3.75) 

(3.76) 

(3.77) 

(3.78) 

(3.79) 

(3.80) 
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leads to the off-line optimization problem for the function z(a): 


Minimize: J, 


subject to: 


= v ; z(0) = 0 ; z(l)=l 


-z(a) < 0 

The Hamiltonian associated with this problem is: 


H = — i- av 2 + A,v — riz 
a 


where T| is a Lagrange multiplier. This results in the optimality condition for v: 


and in the costate equation for X: 


dX 1 

— = — + T1 
da a 


The solution for X is: 


X = In — f-qa 
a 


(3.87) 


where A is an integration constant. Combining Eqs. (3.85) and (3.87) with Eq. (3.82) one 
obtains the differential equation for z: 


^ = .m In |A|_n 


da V2a J a 2 
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which has the general solution: 


z 


4 





a 


J 



(3.89) 


B being a seconjd integration constant. A quick inspection now reveals that there are no 
values of A, B, T|, for which the function z(a) given by Eq. (3.89) can satisfy both of the 
boundary conditions z(0)=0, z(l)=l. Therefore, the off-line optimization problem posed by 
Eqs. (3.81) through (3.83) has no minimum. As guaranteed however by Theorem 4.5, this 
problem does have an infimum, since it is bounded below by J min , where J min is the 
minimum value of the performance index in Eq. (3.69), that is. 



(3.90) 


In order to see what kind of controller u the function z(ot) found in Eq. (3.89) implies, one 
can use Eqs. (3.76), (3.80), and (3.88) to find: 

u = -i(. + ln|A|)-^ (3.91) 

For T|=0, comparing with Eq. (3.72), it is seen that Eq. (3.91) implies the optimal, finite- 
time controller found before. Indeed, Eq. (3.91) can be shown to lead to Eq. (3.72) if the 
boundary conditions y(0)=0 and y(2)=l are imposed and the equation dy/dt— u is 
integrated. There is no contradiction however with either Theorems 3.3 or 3.4, since there 
is no function z in Z corresponding to the finite-time controller u given in Eq. (3.91). 

The above result implies that the minimum value of J, 1^,,= 4/3, found in Eq. (3.90) 
is not only a lower bound of J, but the actual infimum itself. This situation is somewhat 
reminiscent of H-infinity control theory for linear systems in which one seeks a strictly 
proper, stabilizing controller to minimize the H-infinity norm of a closed-loop transfer 
function. The minimum of this norm over all strictly proper, stabilizing controllers does not 
exist. Its infimum does exist, however, and corresponds to a proper controller. Thus, just 
as in H-infinity theory, the question arises naturally in the present case as to how one can 
find an asymptotic controller that somehow approximates this infimum. Although there are 
no general answers yet, a procedure will now be presented that does supply some insight, 
at least for the above example. 
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First, one can show using Eq. (3.88), that, for the off-line optimization problem 
posed by Eqs. (3.81) through (3.83), the following integral Cn=0 when the trajectory does 
not ride the constraint): 



(3.92) 


representing the "total energy" stored in the signal v(a), diverges. This suggests that if one 
imposes the isoperimetric constraint: 


J v 2 da = k (3.93) 

0 

on the off-line optimization problem posed by Eqs. (3.81) through (3.83), where k is a 
given, finite, strictly positive number, one may have a hope of finding a function z(a) in Z, 
that is, one that does satisfy the boundary conditions z(0)=0, z(l)=l, and the state 
constraint, Eq. (3.83). The Hamiltonian associated with this new problem, posed by Eqs. 
(3.81) through (3.83), and Eq. (3.93) will read: 

H = — + (|i + a)v 2 + Xv — t|z (3.94) 


where p is a constant Lagrange multiplier. For each value of k there corresponds a specific 
value of p and vice versa. The corresponding optimality condition for v now becomes: 


X_ 

2(p + oc) 


(3.95) 


The costate equation for X remains unchanged as in Eq. (3.86) and leads to the same 
solution for X given in Eq. (3.87), while the differential equation for z now reads: 


dz _ 1 ^ 

A 

r|a 

da 2(p + a) 

a 

2(p + a) 


(3.96) 


With T|=0 Eq. (3.96) leads to the solution: 
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(3.97) 




a 

f i , 

A 

J(n+9 

0 





satisfying the boundary condition z(0)=0. The boundary condition z(l)=l is satisfied by 
the choice of A (as a function of |i.) that guarantees that: 
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In 




(3.98) 


Finally, the particular value of |t is evaluated from Eq. (3.93) once a value for k has been 
specified. The state constraint, Eq. (3.83), will still have to be satisfied, but forgetting 
about that for a moment, it is possible to show that, for nonzero ji, the improper integrals 
on the right-hand-sides of Eqs. (3.97) and (3.98) do exist. Here, a small indication will be 
supplied to convince the reader that this is indeed so: Consider the value of k corresponding 
to the value p=l. Then, the integral on the right-hand-side of Eq. (3.98) can be evaluated 
easily 33 and, Eq. (3.98) can be solved for A to yield: 


ln|A| = - 


(24 + 7t 2 ) 
121n2 


(3.99) 


The corresponding asymptotic controller u can be found from Eqs. (3.76), (3.78), (3.80), 
(3.96), (3.99), and the fact that p=l. The result is: 


e~‘ f 24 + 7t 2 ^ 

2(l + e-)[ 121n2 \ 


(3.100) 


It can be shown that this controller drives y asymptotically from y=0 at t=0 to y=l as t 
approaches infinity. This fact, however, is already guaranteed by Theorem 3.3. To 
construct an asymptotic feedback controller, one will have to use Eq. (3.100) to integrate 
the equation dy/dt=u. This will lead to an expression for y as a function, of t. Elimination of 
t between this expression and Eq. (3.100) will lead to the desired asymptotic feedback 
controller. Or, alternatively, one can use the transformation defined by Eq. (3.74) with 
z(a) supplied by Eq. (3.97) and express a as a function of y. Then, using Eqs. (3.96) and 
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(3.76) with P=-a should result in the same asymptotic feedback controller as before. 
Although the procedure cannot be carried out analytically for either case, it can be carried 
out numerically. 

It will now be shown that, if one does not impose the state constraint, Eq. (3.67), on 
the off-line optimization problem, Eqs. (3.64) through (3.67), then this problem may not 
have a minimum. Indeed, consider again Example 3, and the associated off-line 
optimization problem without the state constraint, that is, Eqs. (3.81) and (3.82), but not 
Eq. (3.83). Consider now the control: 



if 0 < t < c (3.101) 


u (t) = -(b- l)e“ (, ' c) if c < t (3.102) 


where b is a positive parameter that exceeds 1, and c=b 2 . This control does violate the state 
constraint, Eq. (3.71), by driving y (see Eq. (3.70)) from y=0 at t=0 to y=l as t 
approaches infinity according to: 



if 0 < t < c (3.103) 


y(t) = l + (b-l)e (t_c) 

if c < t 

(3.104) 

The function z(a) corresponding to Eqs. (3.101) through (3.104) can be found from Eqs. 
(3.32), (3.47), and (3.74) as: 

i \ , In a 

z(a) = 1 + 

b 

if 0<ln(l/a)<c 

(3.105) 

N 

II 

cr 

1 

n 

R 

if c<ln(l/a) 

(3.106) 


Substituting from Eqs. (3.101) through (3.104) into Eq. (3.69), one can evaluate the 
corresponding value of the performance index as a function of b, which turns out to be: 
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From Eq. (3.107) it is seen that as b tends to infinity, J tends to minus infinity, 
meaning that the problem imposed by Eqs. (3.81) and (3.82) does not have a minimum. 
Note, however, that there is a value of b, for which the corresponding value of J is the 
same as the optimal value of J (see Eq. (3.90)) for the problem posed by Eqs. (3.69) 
through (3.71), with the state constraint included! This value of b can be found by solving 
the cubic polynomial equation: 


2 



4 

3 


The result is 33 : 


b = 



(3.108) 


(3.109) 


which is approximately equal to: 


b=l. 30497 (3.110) 

The maximum value of y is always equal to b, meaning that, if one regards the constraint 
given by Eq. (3.71) as a soft constraint, then one can recover the optimal value of J with 
an asymptotic, two-sided controller, at the expense of violating the constraint by no more 
than 0.30497. 

3.3.2 Second Reformulation: (Optimal) Asymptotic, Two-Sided 
Controllers that Violate the State Constraint 

The above results hint toward the formulation of the following second off-line 
optimization problem: One again uses the transformation defined by Eqs. (3.55), (3.56), 
(3.61), (3.63) and by: 


a = e~ ( ’~' o) ; dcc = -adt (3.111) 

Then, two functions 4* and M are selected, and the following companion problem is cast: 
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Minimize: 


(3.112) 


I = 'P(x(0),h(x(0),0)) + J M(x,z, v)da 

o 


subject to the isoperimetric constraint (the original performance index now plays the role of 
an isoperimetric constraint): 

i 

0(x(0),h(x(0),0)) + J^jL(x,h(x,z),k(x,z,-av))da = J mn (3.113) 
0 


and the dynamics: 


dx 

da 



x,h(x,z),k(x,z,-av 


)) 


(3.114) 


dz _ 
da 


(3.115) 


x(l) = x 0 ; z(0) = 0 ; z(l) = -S(x 0 ,y 0 ) (3.116) 

Note that there is no state-variable inequality constraint anymore! Specifically, with 
this second reformulation one hopes to find controllers that although will violate the given 
state constraint, will result in the same optimal value of the original performance index, 
corresponding to the problem posed by Eqs. (3.48) through (3.52), which has been 
denoted on the right-hand-side of Eq. (3. 1 13) by J mjn . 


3.4 Conclusions 

The class of all piecewise continuous (in time) controllers that track a given 
hypersurface in the state-space of a dynamical system can be split into two disjoint classes. 
The first class contains all controllers that track the hypersurface in finite-time. The second 
class contains all controllers that track the hypersurface asymptotically. This splitting of the 
two classes can be used to reformulate optimal control problems involving active state- 
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variable inequality constraints, so that, in the reformulated problem the optimization is 
carried out only over the class of asymptotic controllers. 

The original problem leads to optimal controllers that are finite-time and one-sided. 
The reformulated problem leads to controllers that are approximately optimal, asymptotic, 
but still one-sided. However, if the state-constraint is regarded as a soft constraint, then 
one can find controllers that are asymptotic, two-sided, and result in the optimal value of 
the performance index corresponding to the original problem. Thus, they are practically 
optimal, but at the expense of violating the state-constraint. From a singular perturbations 
point of view this suggests that, such controllers can be used in a boundary-layer system, 
to track the reduced solution corresponding to a specific problem, when this reduced 
solution happens to ride a state-constraint boundary. However, such controllers do not 
correspond to stationary solutions of the optimization problem, so at the present, a 
systematic procedure for finding them does not exist and remains a topic for future 
research. 
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Chapter IV* 


Spacecraft Motion under Continuous Thrust 

4.1 Introduction and Motivation 

Exact analytic solutions are scarce in the field of Space Mechanics, but not as scarce 
as in other fields. Until most recently 34 - 36 there have been only three cases in this field in 
which the equations governing the two-dimensional translational motion of a point mass, 
under the influence of inverse-square gravitational forces, and in the presence or absence of 
continuous thrust, have been analytically integrated or reduced to quadratures. The first 
case 37 pertains to Isaac Newton's celebrated solution of the two-body problem and hardly 
needs any comment. The corresponding motion has been called Keplerian, as a tribute to 
Johannes Kepler and his laws of planetary motion. The second case 38 , which is not so 
well-known, corresponds to Euler's reduction to quadratures of the problem of a point 
mass moving on a given plane, and in the vicinity of two fixed, inverse-square centers of 
gravitational attraction. Contrary to the first case, the second case has to date largely 
remained just a mathematical contribution, without any practical applications resulting from 
it. Both of the above cases relate to the natural motion of a point mass in the absence of 
thrust. By contrast, the third case relates to the two-body, planar motion of a point mass, 
under the action of constant radial thrust. The reduction to quadratures corresponding to the 
third case was first given by Tsien 39 , and a more detailed treatment was supplied by 
Battin 34 . It must be noted here that one can find, in general, many ad hoc thrust programs 
for which the equations of motion can be reduced to quadratures. However, the practical 
usefulness of the corresponding solutions usually tends to be minimal due to serious flaws, 
such as the inability to satisfy given boundary conditions, excessively large thrust levels, 
noneconomical fuel consumption, etc. In fact, the third case happens to be more or less just 
such a case, and as far as the author knows, it is the only such case that has received some 
attention in the literature 34 . The thrust program (constant radial thrust) used by Tsien is 
quite ad hoc, and accordingly there is little that can be done with it during actual orbital 
operations. For example, the angular momentum of a vehicle about a planet cannot be 
changed using radial thrust! The only problem that can be treated using Tsien's solution 

* This Chapter and the Appendices have been sole-authored by Nikos Markopoulos. 


45 



turns out to be the problem of escape from a circular orbit. Due to the radial character of the 
thrust such escape usually turns out to be markedly non-optimal 39 . Tsien, being fully aware 
of these facts, treated also in the same paper the problem of escape from a circular orbit, 
using constant transverse (horizontal) thrust, and concluded by observing that the latter was 
much more economical than the former. The equations of motion however cannot be 
reduced to quadratures for the latter type of thrust. 

Tsien’s paper turned out to be just one among a series of papers, that appeared during 
the same period, in which the authors focused primarily on the problem of optimal escape 
from a circular orbit, and compromised in obtaining only approximate analytic solutions, 
but at the significant advantage of using more reasonable (and useful) thrust programs. The 
papers by Benney 40 , who assumed tangential thrust, by Lawden 41 - 42 , who determined the 
optimum thrust direction for minimizing expenditure of rocket propellant, and by Long 43 , 
who studied the possibility of escaping along hyperbolic orbits, are among the notable ones 
in this series. A recent review written by Lawden 35 supplies an extensive reference list on 
this subject. 

Lawden 41 - 42 found that, for all practical purposes, optimum thrust is in the flight-path 
direction, as assumed by Benney 40 . Such thrust will usually be referred to as "tangential" in 
this chapter. However, as in the solution found by Tsien 39 , corresponding to constant 
transverse thrust, the detailed solutions for the flight path and mass loss obtained by 
Benney 40 and Lawden 4142 apply only to cases of large or small thrust. Moreover, in the case 
of small thrust, the solutions obtained are valid only in the initial and intermediate portions 
of the escape trajectory. As noted by Lawden 42 , the assumptions required to obtain such 
solutions are invalid in the final portion of the trajectory, as escape speed is approached and 
the instantaneous or osculating ellipse can no longer be considered close to a circle. More 
recently, Boltz 36 , rather than considering a constant value of the tangential thrust 
acceleration, assumed that the ratio of the thrust to vehicle weight in orbit is fixed. His 
approximate analytic solution obtained with this constraint is valid for any constant value of 
this ratio and describes the motion along the full extent of the escape trajectory. 

This chapter shows how to use the tangential thrust assumption, that apparently 
played a key role in many studies in the past, to propose a continuous thrust program for 
which the equations of two-dimensional motion of a space vehicle in orbit afford an exact 
analytic solution. To come up with this thrust program, the following important aspect of 
optimal control theory is utilized: Optimal control formulations are extremely useful not 
only because they result in control programs that are optimal, but also because they result in 
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control programs period! For a given problem, a control program representing the absolute 
optimum can usually be obtained only by solving a two-point boundary value problem by 
numerical methods. Sometimes however, it turns out to be advantageous to sacrifice 
(hopefully a small) part of optimality, by means of a reasonable approximation, for the sake 
of obtaining a control program that is simple enough to afford an analytical description, yet 
at the same time sophisticated enough to have the ability to steer the system in its state- 
space in a satisfactory manner. An approximation altering the state equations results in state 
trajectories that are approximate. An approximation on the other hand altering the costate 
equations affects only optimality ! It results in state trajectories that are exact , but nearly- 
optimal at best! 

The optimization problem that will be used to apply the above idea is one that has an 
analytically most convenient performance index. The problem of minimum-fuel, power- 
limited transfers between two coplanar (Keplerian) conics comes with a quadratic 
performance index in the control (thrust acceleration) and has proven in the past to be 
amenable to analytical treatment. For transfers between arbitrary elliptical orbits this 
problem has been traditionally attacked using the so-called averaging method 44 - 5 ', which is 
based on the assumption, that, when the duration of transfer is long enough the orbital 
elements are slowly changing, so that, over any given revolution around a planet their 
changes can be computed using rates of change averaged over the mean anomaly. This 
averaging solution provides many useful insights when both the initial and final conics are 
ellipses of low eccentricity, but, as its underlying assumption implies, its accuracy 
decreases as the eccentricities of the elliptical orbits increase (and the motion around the 
planet becomes highly nonuniform with time), and it breaks down when the vehicle does 
not revolve around the planet, that is, when either the initial, or final conics, or both are 
hyperbolic. Also, the averaging eliminates any information about the dependence of the 
solution on the true anomaly. A different approach is clearly needed if one is to improve on 
the accuracy of the averaging method, and to also account for the cases in which either one 
or both of the initial/final conics are hyperbolic. Numerical optimizations 44 - 45 ' 47 49 - 5 ' for such 
problems tend to suggest that there exists a large subclass of power-limited transfers 
between two coplanar elliptic orbits, corresponding mainly to cases where the changes in 
orientation and eccentricity are small, for which, as the transfer duration increases, the 
thrust levels decrease, and the direction of the optimal thrust acceleration tends to coincide 
more and more with the direction of the tangent to the optimal flight path throughout the 
transfer. This last observation once again suggests that tangential thrust is the rule rather 
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than the exception in spaceflight. Although there exist cases in which the thrust is markedly 
nontangential, it seems that such cases arise primarily from a requirement to satisfy given 
boundary conditions in a prescribed time duration, and only secondarily from optimality. 

With this last remark, one can embark on the analysis. The solution of any problem 
starts before its formulation, so, first, in Section 4.2, the relevant equations of motion are 
written in a new state-space format, for which the right-hand-sides of the equations assume 
quite a simple form, free of highly nonlinear trigonometric or transcendental functions. The 
proposed thrust program is introduced in Section 4.3, and the corresponding exact analytic 
solution of the equations of motion is then immediately supplied in Section 4.4. An 
important quadrature is performed in detail in Section 4.5 taking advantage of the fact that 
the trajectory equation that was obtained in Section 4.4 is identical in form to the trajectory 
equation corresponding to Keplerian motion. The fact that the thrust program proposed in 
Section 4.3 is not ad hoc is documented in Section 4.6, where it is shown that one can 
obtain this thrust program by using the tangential thrust assumption in one costate equation 
corresponding to the problem of power-limited optimization. This problem of power- 
limited optimization is formulated in Section 4.6.1, after which the consequences of the 
tangential thrust assumption are followed in Section 4.6.2. The corresponding 
transversality conditions are derived from scratch in Section 4.6.3, and they are combined 
with the tangential thrust assumption in Section 4.6.4. All this leads to the thrust program 
proposed in Section 4.3, for which the state equations can be solved analytically and 
exactly, leading to the solution supplied in Section 4.4. Section 4.6 concludes with a 
preliminary check of the near-optimality of the tangential thrust assumption (in 4.6.5). 
Section 4.7 casts the boundary conditions in a ready-to-use form. These boundary 
conditions are manipulated further in the next three sections during the treatment of general 
transfer and escape problems. Examples are also given in Sections 4.8, 4.9, and 4.10 
corresponding to specific transfer and escape problems, in which the existence of particular 
solutions for such problems is demonstrated for several cases. The existence of one- 
segment solutions, that is, solutions on which the throttling parameter A (see 
Nomenclature) has a single value is discussed in Section 4.1 1, along with some geometric 
facts pertaining to (one-segment) transfer trajectories. Section 4.12 documents the 
practically important fact that, using multiple segments, a given coplanar transfer can in 
principle be performed in an arbitrarily large preassigned duration. Questions about the 
optimality of the transfer trajectories, both from the point of view of power-limited and 
constant ejection velocity propulsion systems are discussed in Section 4.13. Section 4.14 
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constitutes a short deviation from the main subject and supplies a (hopefully reasonable) 
prediction for the initial values of the costates for the associated rendezvous problem 
pertaining to power-limited optimization. After a brief summary in Section 4.15, the 
chapter concludes with some comments in Section 4.16. Preliminary results on the topic 
presented in this chapter and the Appendices have been published in Ref. 52 (See also Ref 
53). 

4.2 Equations of Motion 

Written in polar coordinates, the equations of two-dimensional motion for a space 
vehicle, viewed as a point mass, and moving in the vicinity of an inverse-square 
gravitational center of attraction (a planet, the sun, etc.) are: 



_d 

dt 




= REe 


(4.2) 


where, R and 0 specify the position of the vehicle (see Fig. 4.1), t is the real time, (i is the 
strength of the inverse-square center, and E R and E e are respectively the sums of the radial 
(vertical) and transverse (horizontal) components per unit mass of any non-gravitational 
forces acting on the vehicle. For the purposes of the present work, the only non- 
gravitational force acting on the vehicle will be the thrust, specified by its components 
(controls) per unit mass Er and E e . 

The goal now is to cast Eqs. (4.1) and (4.2) into a simple state-space form. It turns 
out to be advantageous for this purpose to define the relevant state variables in such a way 
that the resulting state equations do not contain any (highly nonlinear) transcendental 
functions and can be written only in terms of rational functions of the state variables. 
Although there are many ways of accomplishing this, one simple way is to introduce the 
angular momentum per unit mass of the vehicle, M, and the vertical component of its 
velocity X, defined by: 


M = R 2 




(4.3) 


and use M, X, R, and 0 as the governing state variables. 
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It would also be convenient to have the state equations in nondimensional form, so, 
for this purpose, one can select an arbitrary distance R s and denote the gravitational 
acceleration and circular orbital velocity at that distance by gj and V $ , that is: 

g,=£ ; = < 4 - 4 ) 

K s 

By defining now the nondimensional variables: 



and using h, x, r, and 0 as the new (nondimensional) state variables, one can write 
Eqs. (4. 1 ) and (4.2) in the following equivalent nondimensional state-space form (see Fig. 
4.1): 


dh 



(4.6) 


dx 

dt 


(V-0 

< 1-3 ; 


+ £, 


(4.7) 


dr 

dx 


= x 


(4.8) 


d6 = h_ 
dx r 2 


(4.9) 


4.3 The Proposed Thrust Program 


Consider the thrust (acceleration) program given explicitly by the components: 



(4.10) 


where A is an arbitrary constant that will be called the throttling parameter. It is 
straightforward to show that the thrust acceleration corresponding to this program has 
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magnitude equal to A/(rcosy) and acts always tangentially to the vehicle's path. Equation 
(4.10) therefore implies a special type of tangential thrust. Also, it is clear that the 
mechanical energy per unit mass of the vehicle (considered as a point mass) is a strictly 
increasing function of time for A > 0 and a strictly decreasing function of time for A < 0. 
When A = 0 there is no thrust and the motion is Keplerian (see Appendix A). It will be 
shown in the next section, that Eqs. (4.6) through (4.9) afford an exact analytic solution 
when their right-hand-sides are forced with the thrust (acceleration) components given in 
Eq. (4. 10). This chapter shows how one can come up with this thrust program and use the 
corresponding exact analytic solution to satisfy boundary conditions for problems of 
practical interest. 


4.4 An Exact Analytic Solution 

Substituting the thrust acceleration components proposed in Eq. (4. 10) into the right- 
hand-sides of Eqs. (4.6), (4.7), the state equations assume the form: 


dh 

dx 


= A 


(4.11) 


dx 

dx 


^h 2 -r N ) Ax 

r 3 J h 


(4.12) 


dr 

— = x 
dx 


(4.13) 


d0 = _h 
dx r 2 


(4.14) 


The exact analytic solution of the system of Eqs. (4.11) through (4.14) is derived in 
Appendix B. One can verify easily by direct differentiation that this exact analytic solution 
is given by the following four quadratures: 

h = Ax + h 0 (4.15) 
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_ 2rA + Bsin(6-C) 
x _ - - 


(4.16) 


h 2 

1 + Bcos(0-C) 


(4.17) 


e 

J_ 

h 2 

where h,,, B, C, and 0 O are the four (exact) constants of integration associated with the 
system of Eqs. (4.11) through (4.14). The procedure carried out in Appendix B also 
implies the uniqueness of the expressions given in Eqs. (4.15) through (4.18), that is, that 
every particular solution of the system of Eqs. (4.11) through (4.14) corresponds to a 
particular set of values for the constants 1^, B, C, and 0 O . Note that the constants 1^ and 0 O 
are equal respectively to the initial values (at x=0) of the angular momentum h, and the 
argument of latitude 0. Note also that Eq. (4.18) is obtained by combining Eqs. (4.11), 
(4.14) and using Eq. (4.17) to substitute for r, which leads to (the throttling parameter A is 
always assumed nonzero): 


-H 


2Ad0 


[l + Bcos(0-C)f 


(4.18) 


d0 [l + Bcos(0-C)] 2 
dh ” Ah 3 


Equation (4.17), describing the trajectory of the vehicle, is identical in form with the 
trajectory equation corresponding to two-body Keplerian motion (see Appendix A). Thus, 
the trajectory equation corresponding to Keplerian motion can be uncovered in the presence 
of continuous, nonzero thrust! 

There is of course an important difference between the two kinds of motion. 
Keplerian motion takes place in the absence of thrust, and along it the angular momentum h 
is just a constant. On the other hand, the motion described by Eqs. (4.15) through (4.18) 
takes place in the presence of continuous tangential thrust (as specified in Eq. (4.10)), and 
along it h is a linear function of time (see Eq. (4.15)). 
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Based on the above observation, the explicit integration of Eq. (4.18), that depends 
strongly on the value of the constant B, will be performed in some detail in the next section 
by generalizing the concepts of eccentric and hyperbolic anomalies encountered in 
Keplerian motion. The trajectory described by Eqs. (4.15) through (4.18) will from now 
on be called the transfer trajectory, or just the transfer. 


4.5 Generalized Eccentric and Hyperbolic Anomalies 
4.5.1 Preliminaries 

Recall the expression for r given in Eq. (4. 17) and denote 0 - C by £: 

£ = 6-C (4.20) 


Then, Eq. (4.17) can be written as: 


h 2 

1 + Bcos£ 


(4.21) 


The angular momentum h appearing in Eq. (4.21) is the linear function of time given 
by Eq. (4.15). In order to find the relationship between the time T and the argument of 
latitude 0 along the transfer trajectory one must integrate equation (4.14). In this case 
however, since h is a linear function of time, and since the interest is on motion with 
nonzero A, it is more convenient to find the relationship between h and 0. Then, the 
relationship between x and 0 can be deduced. Thus, co mbinin g the equations: 


d0 = _h dh 
dx r 2 ’ dx 


(4.22) 


and using: 


d4 = d0 ; $o=9„tC 


(4.23) 


one finds after substituting from Eqs. (4.21): 


1 

h 2 


_1__ f_2Ad|_ 

h^ J (l + Bcos^) 2 

^0 


(4.24) 
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The constants B and C appearing in Eq. (4.17) play qualitatively the same role as the 
eccentricity and orientation constants e and to do in Keplerian motion 34 (see Appendix A). 
The vehicle, while on the transfer trajectory, can execute full revolutions around a planet 
only if B is less than one. B can always be taken as nonnegative by appropriately adjusting 
the constant C (just as e can always be taken as nonnegative by appropriately adjusting the 
orientation to of a conic). Note that escape from a given orbit can be accomplished for any 
value of B, since the angular momentum h, and accordingly the radial distance r can be 
made to grow without bound if the thrust is kept "on" for a sufficient duration of time. 

Based on the above observations, one can call the constants B and C the generalized 
eccentricity and the generalized orientation of the transfer trajectory respectively. The 
constants B and C are global characteristics of the transfer trajectory, that is, they are 
associated with every point of the transfer trajectory. If at any point along the transfer the 
thrust is suddenly reduced to zero, the eccentricity e and orientation CD of the ensuing 
Keplerian motion are not equal to B and C respectively. This follows by equating the 
expressions for x and r of Eqs. (4.16), (4.17) with the corresponding expressions for x 
and r of Eq. (A.5) (of Appendix A), valid for Keplerian motion, and by remembering that 
the throttling parameter A is nonzero. Note however, that, for small A, the (Non- 
Keplerian) constants B and C are approximately equal to the (instantaneous Keplerian) 
constants e and (0. 

Carrying the above analogy further facilitates the evaluation of the integral appearing 
in Eq. (4.24). Traditionally, the corresponding integral for Keplerian motion is evaluated 
by defining the so-called eccentric and hyperbolic anomalies 34 , depending on whether the 
eccentricity e is less than or greater than one respectively. By analogy, it is possible to 
evaluate the integral in Eq. (4.24) by defining generalized eccentric and hyperbolic 
anomalies, depending on whether B is less than or greater than one. In the borderline cases 
in which B=0, or B=1 the evaluation of the integral simply yields 33 : 

if B=0, then: X = X~ 2A(0-0 O ) ( 4 - 25 ) 

h h 0 


ifB=l,then: -i- = 4---A 

h 2 h 2 


tan 

rq 

-tan 

(l V 

So 

1-f-) 

tan 3 

ill 

- tan 3 

(l Y1 

So 


UJ 


12 J. 

J UJ 

_ 

U ) 


< 2 L 


(4.26) 


Such borderline cases represent only mathematical but not "real" possibilities, and it 
is only for reasons of mathematical completeness that one considers them at all. In practice. 
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when a vehicle is executing an actual transfer, or coasting along a conic, the probability that 
e or B is exactly equal to any particular value is practically zero. Borderline cases can 
therefore also be considered as limiting cases, in which the results are first worked out in 
general (for e or B less than or greater than one) and then the limit is taken in which e or B 
tends to zero, one etc. 


4.5.2 Generalized Eccentric Anomaly 


Consider first those transfer trajectories on which 0 < B < 1. For all these trajectories 
one can define the generalized eccentric anomaly E by: 


cos£ = 


cosE-B 
1 -BcosE 


d$ = 


Vl-B 2 

1- BcosE 


dE 


(4.27) 


The expression on the right in Eq. (4.27) relates the differentials d£ and dE. One can make 
the convention that E is zero when £, is zero and vice versa. Then, E goes through 7t radians 
exactly when £, goes through k radians. This means that one can keep track of the number 
of times a vehicle revolves around a planet (during a given transfer) by just keeping track of 
E and forgetting about ^ (or 0). Using standard trigonometric identities 33 one can show 
that: 


. c Vl-B 2 sinE 

sine = 

1- BcosE 



(4.28) 


Also, using Eq. (4.27) one can write the transfer trajectory equation (4.17) as: 


r = 


h 2 ) 

U-B 2 ; 


(1- BcosE) 


(4.29) 


The inverse of expressions in Eqs. (4.27) and (4.28) can be found as: 


cosE = 


cos^ + B 
1 + Bcos^ 


dE = 


Vl-B 2 

1 + Bcos^ 




(4.30) 
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(4.31) 


sinE = 


Vl-B 2 sinE, 
1 + Bcos2; 



while the initial and final values of E, E 0 and E f , are given by: 


cosE 0 = 


cos^q + B 
1 + Bcos^ 0 


cosE f = 


cos^ f + B 
l + Bcos^ f 


Eq. (4.24) can then be written as: 


1 

h 2 



2A 


O-B-f 


E 

J(l-BcosE)dE 


which after a straightforward integration results in: 

1 i 2A[E-E 0 -B(sinE-sinE 0 )] 

i^ = hl (Wf 5 


(4.32) 


(4.33) 


(4.34) 


4.5.3 Generalized Hyperbolic Anomaly 


For transfer trajectories on which B > 1 one can define the generalized hyperbolic 
anomaly H by: 


B-coshH 
BcoshH- 1 


d$ = 


Vb 2 -i 

BcoshH - 1 


dH 


(4.35) 


The expression on the right in Eq. (4.35) relates the differentials d£ and dH. One can again 
make the convention that H is zero when ^ is zero and vice versa. Using standard 
hyperbolic identities 33 one can show that: 


Vb 2 - IsinhH 
BcoshH -1 ’ 



(4.36) 


Using Eq. (4.35), the transfer trajectory equation (4.17) can now be written as: 
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r = 


( h 2 ^ 

V B 2 - 1 j 


(BcoshH-1) 


(4.37) 


The inverse of expressions in Eqs. (4.35) and (4.36) can be found as: 

cosh H = 7 “—— “r- ; dH-^U 

1 + Bcos^ 1 + Bcos^ 


(4.38) 


sinh H = 


_ Vb 2 - lsin^ 

1 + Bcos^ 


tanh (f) = il^i tan (f) <439) 


while the initial and final values of H, Hq and H,, are given by: 


coshH 0 =-£^±B coshH[ = i2!lL±B 

1 + Bcos^ 0 f l + Bcos£ f 


(4.40) 


Eq. (4.24) can then be written as: 


J___l 

.2 “ U 2 


2A 


h 2 h 2 (B 2 -l) 3/2 J 


H 

J*(BcoshH-l)dH 


(4.41) 


which after a straightforward integration results in: 


11 2A [H - H 0 -B (sinh H- sinh H 0 )] 

vT7 “ 7T + 717 T372 


h h 


(B 2 - 1) 


(4.42) 


4.5.4 Examples of Motion under the Proposed Thrust Program 

At this point one can begin to consider some preliminary examples depicting the kind 
of motion described by the trajectory equation (4.17), in conjunction with Eqs. (4.25), 
(4.26), (4.34), and (4.42) derived in this section. Figures 4.2 through 4.4 show examples 
of transfer trajectories (plots of Eq. (4. 17)) that start from the perigee of the same initial 
elliptical orbit and correspond to different values of the throttling parameter A and to 
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different duration x f . Figure 4.5 shows a similar case for which the initial conic is a 
hyperbola. In each figure the final conic, corresponding to the motion that follows after the 
thrust is turned off at x=x f , is also plotted. It is seen from Eqs. (4.17), (4.25), (4.26), 
(4.34), and (4.42) that a vehicle can escape the gravitational field of a planet even along a 
transfer trajectory for which the generalized eccentricity B is less than one, providing that 
the thrust is kept "on" for a sufficient duration of time. Note that the distance from the 
planet r becomes infinite on a transfer trajectory (for any B) when, and only when the 
angular momentum h becomes infinite. This happens when the right-hand-side of the 
applicable one among Eqs. (4.25), (4.26), (4.34), or (4.42) becomes zero at some 0. 
Figure 4.4 shows an escape trajectory on which B < 1. 


4.6 The Origin of the Thrust Program 


In this section a connection will be developed between the thrust program proposed in 
Section 4.3, and the problem of optimal, power-limited, coplanar motion. It will be shown 
that, the thrust program proposed in Eq. (4.10) can be obtained by making the tangential 
thrust assumption and by imposing the transversality conditions in the costate system 
corresponding to the problem of power-limited, minimum-fuel transfers between two 
coplanar conics. 

A Mayer type formulation of the above problem is the most convenient one, and for 
this, one will have to add to the system of Eqs. (4.6)-(4.9) the equation governing the 
performance index: 

(4 - 43) 

The final value of J (at a fixed final time x f ) is a direct (nondimensional) measure of 
the fuel consumption for power-limited propulsion systems 44 - 51 (see Appendix C for more 
details). Note that the variation for the mass of the vehicle during a maneuver is fully taken 
into account by means of Eq. (4.43) (see Eq. (C.5) of Appendix C). 


4.6.1 The State-Costate System for Power-Limited Optimization 


The Hamiltonian associated with the system of Eqs. (4.6) - (4.9), and (4.43) is: 


H = re 0 P h + 


^h 2 -r^ 


P„ +e,P„ +xP, +£ p, + i(ef +£;)!>, (4.44) 
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In order to have minimum fuel consumption one must minimize J f - J 0 , or maximize -J f + j 
The cost J is an ignorable coordinate, meaning that the corresponding costate Pj is constant 
(Specifically, the Hamiltonian does not depend explicitly on J). Moreover, examination of 
the transversality conditions in Section 4.6.3 will further dictate that (see Eq. (4.98)): 

Pj = -1 (4.45) 

With this result, the optimal controls can be found by setting the partial derivatives of H 
with respect to e r , £g in Eq. (4.44) equal to zero. The result is: 

e; = P* ; e; = rP h (4.46) 

Substituting these controls back into Eq. (4.44) the Hamiltonian can be written as (using 

P 


H‘=i( 

2 X 


r P; 


+P 2 )- 


h 2 -r 


+ xP r +-P e 


(4.47) 


Using this expression for the Hamiltonian, the differential equations governing the 
optima] trajectories can be found from: 


ds _ 3H* dP _ 3H* 

dt 3P ’ dt Bs 


(4.48) 


where s and P denote the state and costate vectors respectively. Explicitly, the result is the 
following tenth order system for the state and costate equations: 


dh 

dt 



(4.49) 


dx f h 2 — r _ 
dt”( r 3 ) ' 

(4.50) 

dr 


3-1 

ii 

X 

(4.51) 

de_ h 


dx r 2 

(4.52) 
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(4.53) 


dJ 

dx 





dt 


2hP 



d^. 

dx 


= - P. 


dP^ 

dx 


-rP h 2 + 


3h 2 P. 2P„ 2hP„ 


+ 


dPe 

dx 


= 0 


(4.54) 

(4.55) 

(4.56) 

(4.57) 


= 0 (4-58) 

dx 

Note that Eqs. (4.49) through (4.52) are just the state equations (4.6) through (4.9), 
with the controls substituted in from Eq. (4.46). It can be verified that this tenth order 
system has the following three first integrals: 

H* = const. ; P e = const. ; P, = const.= -1 (4.59) 

and a fourth integral 48 given by: 

2rP r +hP h -xP x = 3H*x-5J + A 4 (4.60) 

where A 4 is an integration constant. These four integrals can in principle be used to reduce 
the order of the system to six, involving the variables h, x, r, 0, P h , and P x . 

An important point that can be made regarding the state-costate system of Eqs. (4.49) 
through (4.58) is the relative simplicity of the right-hand-sides of these equations. Note that 
these equations contain no trigonometric functions. A quick inspection reveals that the 
number of basic operations (addition, subtraction, multiplication, and division) that need to 
be performed to compute the right-hand-sides of Eqs. (4.49) through (4.58) is about half 
times the number of basic operations that need to be performed to compute the right-hand- 
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sides of the same equations, but written using traditional state variables, like the set (V, y, 
r, 0) or the orbital elements. Pursuing this argument further, one can conclude that the 
computer running time required for the numerical solution of a two-point boundary value 
problem using Eqs. (4.49) through (4.58) should be about half, compared to the running 
time required for the numerical solution of the same problem using the traditional state 
variables (V, y, r, 0) or the orbital elements. 

4.6.2 The Tangential Thrust Assumption and its Implications 

One can now pursue the consequences of the assumption that the thrust along an 
optimal trajectory is approximately tangent to the direction of the optimal flight path. This 
assumption will be referred to as "the tangential thrust assumption". In terms of the state 
and costate variables, the tangential thrust assumption means that along an optimal 
trajectory the relationship: 



rP h 


rx 

— = tany 
h ' 


(4.61) 


is approximately valid. Note that y is the flight path angle, formed by the path of the vehicle 
and the local horizontal direction. It will be assumed, as Eq. (4.61) implies, that the tenth 
order system of Eqs. (4.49)-(4.58) behaves in practice as a ninth order system, meaning 
that one will have to drop one of these equations out. As it turns out, one can combine Eqs. 
(4.54) and (4.61) and solve for P h as a function of the states. Then, P x can be found by 
resorting back to Eq. (4.61), and P r can be found from Eq. (4.55) by differentiating the 
expression for P x . Accordingly, one won’t need Eq. (4.56), which will be the one equation 
to be dropped out. Thus, using the tangential thrust assumption, one can replace the system 
of Eqs. (4.49)-(4.58) by the system: 



rP h h 


(4.62) 


dh 

dt 



dx 

dx 



+ P. 


(4.63) 

(4.64) 
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dr 

— = X 

dx 

(4.65) 

d0 __h 
dx r 2 

(4.66) 

— = -(P. 2 +r 2 P h ! ) 
dx 2' h ' 

(4.67) 

dP h _ 2xP h P e 

dx r r 2 

(4.68) 

dP x 

dt ' 1 

(4.69) 

dP *=0 

dx 

(4.70) 

o 

II 

oT £ 

(4.71) 


Note here that: i) Eq. (4.62) is just the tangential thrust assumption, ii) All the 
remaining equations (4.63)-(4.71) are the same as the original equations (4.49)-(4.55), and 
(4.57), (4.58), except equation (4.68), which was obtained from Eq. (4.54) by changing 
its right-hand side using the tangential thrust assumption, iii) The optimal Hamiltonian H' 
given in Eq. (4.47), and the constant A 4 defined in Eq. (4.60) are no longer first integrals 
of Eqs. (4.62)-(4.71). iv) P e and Pj are still constant and P, = -1. 

One can now show the primary analytical result of this chapter, and also of this 
report, which can be summarized in the following statement: 

For the problem of transfers between two arbitrary conics, if one enforces the 
transversality conditions exactly, then Eqs. (4.62)-(4.71) can be reduced to quadratures 
without any further approximations. 

By defining the transformation: 

Q = r 2 P h (4.72) 
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and then, by using Eqs. (4.65), (4.68) one finds that: 


dQ 

dx 


Pe 


Since P 9 is constant, the solution of Eq. (4.73) is simply: 

Q = -P e x + A 


(4.73) 


(4.74) 


where, A is an integration constant. Combining now Eqs. (4.72), (4.74) and (4.63) one 
obtains: 


dh 

dx 


= Q = -P e x + A 


(4.75) 


the solution of which is: 


h 



+ Ax + h 0 


(4.76) 


where, h,, is another constant of integration, equal to the value of the angular momentum h 
at time x=0. Using Eqs. (4.72), (4.62), and then carrying out the differentiation in Eq. 
(4.69) one can express the costates P x , P h , and P r in terms of the state variables, the time, 
and the three constants P 0 , A, and h^, as: 



The costate variable P 8 represents the sensitivity of the minimum cost to the initial 
value of the state component 0. By examining the first variation of the (augmented) cost in 
Section 4.6.3, it will be shown that for transfers "from" or "to" a circular orbit the exact 
value of P e is zero. This is not surprising, since for such problems, due to the initial 
circular symmetry, the initial value of the state component 0 has no effect on the minim um 
value of the cost. It will also be shown in Sections 4.6.3 and 4.6.4, by enforcing the 
transversality conditions and by using the expressions for the costates given in Eq. (4.77), 
that for transfers between two arbitrary conics, P 6 is zero as a direct consequence of the 
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tangential thrust assumption. With zero P e the variable Q from Eq. (4.75) is just equal to 
the constant A, and the expressions for the costates given in Eq. (4.77) reduce to: 


P = 





U 


(4.78) 


Using the expressions for P x and P h from Eq. (4.78) one can now go back and obtain the 
corresponding explicit thrust acceleration components, by substituting into Eq. (4.46). The 
result is: 






(4.79) 


But this is none other then the thrust (acceleration) program proposed in Eq. (4.10)! 
Therefore, one can state that: 

The thrust program proposed in Eq. (4.10) can be obtained by making the tangential 
thrust assumption in one costate equation and by imposing the transversality conditions for 
the problem corresponding to power-limited, minimum-fuel, coplanar transfers between 
two arbitrary conics. 

Note that the asterisks have been omitted from e r and £ e in Eq. (4.79), since these 
thrust acceleration components are consequences of the tangential thrust assumption and 
can be near-optimal at best. 

Now that the thrust program is "fixed" and the same as the one proposed in Eq. 
(4.10), the exact analytic solution of the state equations is also "fixed" and the same as the 
one given in Section 4.4. For Q = A = constant, Eqs. (4.75), (4.76) lead to: 


— = A ; h = Ax + h 0 (4.80) 


These equations are the same as Eqs. (4.11) and (4.15) of Section 4.4. Using the 
expression for P x given in Eq. (4.78), Eqs. (4.64), (4.65) assume the form: 


dx 

dx 


(\f- 


h -r 


Ax 


dr 

dx 


= x 


(4.81) 
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which are the same as Eqs. (4.12), (4.13) of Section 4.4. The solution of Eqs. (4.81) is: 


2r A + Bsin(6-C) . h 2 

h ’ r_ 1 + Bcos(0-C) 


(4.82) 


which is none other than the one given in Eqs. (4.16), (4.17) of Section 4.4. The last state 
equation that remains to be integrated is Eq. (4.66), which is the same as Eq. (4.14) and 
can be reduced to quadrature as in Eq. (4.18), or integrated as in Section 4.5. 

Denoting by h, the value of the angular momentum h at the final time x f one can write 
the constant A, playing the role of a throttling parameter, and appearing in the linear 
variation of h (Eqs. (4.15), (4.80)) as: 


A = 





(4.83) 


Combining Eqs. (4.79) and (4.83), the two thrust acceleration components along a 
transfer can be written explicitly as: 


(h f - h 0 )x _ (hf-hp) 

r x f h ’ e x f r 


(4.84) 


By replacing 1^ by h and T f by the "time to go" x f -x in Eq. (4.84) one obtains the following 
explicit, closed-loop, finite-time, exact feedback guidance law: 


(hf-h)x 

(Xf-X)h 



(4.85) 


During an actual transfer this law will guide a vehicle from the initial conic to the desired 
final conic, providing that the transfer starts from the "correct" point on the initial conic. 
Since this happens to be a finite -time feedback guidance law, it is bound to suffer from 
saturation problems (in the presence of disturbances) as the final time x = x f is approached. 
In such a case, the reasonable thing to do would be to switch to the open-loop law of Eq. 
(4.84) as late as possible on the transfer, but before the final time x = x f is reached. 

As expected, there is a partial loss of controllability that goes along with the tangential 
thrust assumption, and it will be shown when examining the boundary conditions that, if 
the initial and final conics do not intersect, then, a one-segment solution for a particular 
transfer problem, that is, a transfer trajectory on which the throttling parameter A has a 
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single value may in general exist only for particular values of the (fixed) final time x f . For 
escape problems on the other hand a solution may exist (in principle) for any value of x f 
because the final value of the angular momentum h f , or the orientation difference co f - co 0 
between the initial and final conics may be left free. A one-segment solution however does 
not exist, and a (one-segment) transfer trajectory cannot be found, when the initial and final 
conics intersect (see Theorem 4.1 of Section 4.11, and Appendix F). However, even in 
such a case, a transfer can still be accomplished using two separate transfer trajectories, that 
is, by first transferring to an intermediate conic (that does not intersect either the initial or 
the final conic), and then by transferring to the desired final conic (see Sections 4.1 1, 4. 12 
for more details). 


4.6.3 Derivation of the Transversality Conditions 


For transfer or escape problems the state of the vehicle at the initial and final times is 
only partly specified. Accordingly, the boundary conditions have to be supplemented by 
the so-called transversality conditions, arising from the requirements of optimality, in order 
to provide one with enough conditions for a particular solution. The exact form of the 
transversality conditions will be given in this Section, and then combined with Eq. (4.77) 
of Section 4.6.4, in order to show the result P e = 0, which led to Eq. (4.78) and the thrust 

program proposed in Eq. (4.10). 

There are several methods that can be used to derive the transversality conditions for 
the problems of interest. The one that will be used here is very straightforward, and it is 
based on a direct examination of the first variation of the (augmented) performance index. 
Recall that the performance index is just the negative of the final value of J, that is, -J f . 
Augmenting this by the dynamics of the problem one obtains the augmented performance 
index I: 


I = -J f + 


h 

J P(f-s)dx 


To 


(4.86) 


The costate (Lagrange multiplier) vector P is a row vector with components P h , P x , P r , P e , 
and Pj. The vector of dynamics f is a column vector with components the right hand sides 
of Eqs. (4.6)-(4.9) and (4.43). The vector s is also a column vector with components the 
left-hand sides of Eqs. (4.6)-(4.9) and (4.43). Using the Hamiltonian H=Pf, Eq. (4.86) 
can be written as: 
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(4.87) 


I = -J f + J*[H(s,P,u)-Ps]dx 

*0 

The components of the control vector u are just the two components of the thrust 
acceleration, e r and Eg. The initial time x 0 and the final time x f are always fixed. 

Assume now that the value of the augmented performance index I given in Eq. (4.87) 
is the optimal one, that is, it corresponds to the optimal choice of the control u. Then, if the 
optima] control is perturbed by a small amount 8u to u+8u, the augmented performance 
index I is correspondingly perturbed to I+AI, given by: 

I + AI = -J f -5J f + J[H(s + 8s,P,u + 8u)-P(s + 8s)]dx (4.88) 

T 0 


where, the perturbations in the state components resulting from the perturbation in u have 
been denoted by 8s, and where: 


8s = 8 


ds A _ d(8s) 


,dx 


dx 


From Eqs. (4.87) and (4.88), the total variation of I is: 


AI = -SJ f 



(4.89) 


(4.90) 


The first variation of I is that part of AI corresponding only to the first order change of H 
in the small quantities 8x, 8u: 


81 = -8J f + 



(4.91) 


Integrating the second term in the integral by parts and expressing 8H in terms of 8s and 8u 
one can write Eq. (4.91) as: 
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dT 


(4.92) 


SI = -8J f - P f 5s f + P 0 5s 0 + 



5s + — Su 
3u 


where, the subscripts 0 and f were used to denote the values of the costates P and the 
variations 5s at the initial and final time respectively. The natural choice for the costate 
functions P at this point is: 


P = - 


9H 

3s 


(4.93) 


based on which Si from Eq. (4.92) can be written as: 


Si = — 5J f — Pf5s f + P 0 5s 0 + 



*0 


(4.94) 


Since the original choice of u was assumed to be the optimal one, 51 should vanish for 
arbitrary perturbations Su about u. From Eq. (4.94) a necessary condition for this is: 

— = 0 (4.95) 

3u 

In view of Eq. (4.95) and the requirement that Si be zero, one obtains from Eq. (4.94): 

51 = -8J f - P f 8s f + P 0 8s 0 (4.96) 

The conditions given by Eqs. (4.93) and (4.95) are necessary for optimality. The 
explicit form of Eqs. (4.93) and (4.95) has already been obtained in the text. Eq. (4.95) is 
just the optimality condition used in obtaining Eq. (4.46), while Eq. (4.93) represents the 
costate equations (4.54)-(4.58). Here, Eq. (4.96) will be further manipulated into a form 
that will readily yield the transversality conditions for any problem that one wishes to 
study. First, one can write Eq. (4.96) explicitly, using the components of P and 8s at the 

initial and final times as: 

81 = -(1 + P } )5J f - P hf 5h f - P xf 8x f - P rf 5r f - P e 80 f 


68 



(4.97) 


+P, 5J 0 + P h0 8h 0 + P x0 8x 0 + P r0 6r 0 + P e 80 o 

The initial value of the cost is zero (J 0 =O), meaning that 8J 0 = 0. The final value of the 

cost is free (since the aim in the optimization problem is to find the optimal final value), 
meaning that: 


Pj = -1 (4.98) 

Note that this proves the result claimed in Eq. (4.45). Now Eq. (4.97) is reduced to: 

51 = p h 0 5h 0 + Pxo 5x o + Pro 5r o + P e 500 - P h f 5h f - P xf 8x f - P rf 8r f - P e 80 f (4.99) 


For transfer problems the variations 8h, 8x, 8r, and 80 are not in general independent 
of each other, neither at the initial nor at the final time. However, in all such problems, the 
vehicle, before the initial and after the final times is assumed to coast along conics (circles, 
ellipses, parabolas, or hyperbolas), defined by their orbital elements h, (angular 
momentum), e, (eccentricity), and to, (orientation) (i is 0 or f). For all such problems the 
variables that will be prescribed or left free at the initial or final time are h,, e„ to, and 0,. To 
derive the transversality conditions one will therefore have to cast the remaining part of the 
first variation of I given in Eq. (4.99) in terms of the independent variations 8h,, Se,, 8to„ 
and 80j rather than the dependent variations 8h,, 8x,, 8r,, and 80,. Toward this end one can 
first observe that at any point along a conic (corresponding to Keplerian motion) the radial 
distance r and the radial component of the velocity x are given by (see Eq. (A.5) of 
Appendix A): 


h 2 

l + ecos(0-©) 


esin(0-©) 

h 


(4.100) 


Differentiating Eq. (4. 100) one finds that, to first order, variations in x and r corresponding 
to small variations 8h, 8e, 8©, and 80 are given by: 



r ( h2 ~ r ) 

h 2 e 


8e + 


^-(80-8©) 


(4.101) 



e 2 r 2 -(h 2 -r) 2 

exh 2 r 2 


8e + 


(h 2 -r) 

hr 


(80-8©) 


(4.102) 
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Writing these expressions for 5x and 5r at the initial and final times, one can substitute into 
Eq. (4.99) and obtain the desired form of the remaining part of 51: 


81 - K h0 5h 0 - K hf 8h f + K e0 5e 0 K ef 8e f + K e0 59 0 KgfSQf K^qScOq + K^ScOj (4.103) 


where the variables K are defined as: 


K h =P h + 


2rP -xP 


2 2 

e r 


K e = , , 

e exh 2 r 2 




elf 


(4. 104. a) 


K e = P e + 


hr 


2 u2 

P, + ^-P r ; K = 


If - r „ xf „ h -r„ xr 


P. + — P r 
hr x h r 


(4.104.b) 


and the additional subscript 0 or f represents corresponding values at the initial or final 
time. 

For any transfer problem of interest Eq. (4.103) can be used to pick the desired 
transversality conditions at will. First, for all transfer problems the eccentricities of the 
initial and final conics are always fixed, meaning that, 8e 0 and 8e f in Eq. (4.103) will 
always be zero. Then, whenever one of the six quantities 1^, co 0 , 0 O , h f , co f and 0 f is left 
free in a problem, the coefficient of the corresponding variation of that quantity in Eq. 
(4.103) is set equal to zero, resulting in a transversality condition. For example, if for a 
given problem the angular momentum at the initial and final times is left free, the 
requirement that 81=0 results in the two transversality conditions K^K^O. This result 
has an interesting implication when combined with the fourth integral of motion given in 
Eq. (4.60). It implies that the constant A 4 is zero, and that the exact optimal cost for all such 

problems is given by: 



(4.105) 


In all the (transfer) problems that will be considered here the argument of latitude 9 
will be free both at the initial and final times. From Eq. (4.103), this means that one will 
always have at least the two transversality conditions: 

2 2 

if 0j is free, then: P 9 4- ^ - 1 ' P xi +^-^-P r , = 0 (4.106) 
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where i is 0 or f. Also, the angular momentum of the initial conic will always be specified, 
meaning that 6ho=0. Taking the above considerations into account, using Eqs. (4.106), and 
defining the difference in the orientation between the initial and final conics by ft: 

ft = co f -co 0 ; 8ft = 5o) f - 8co 0 (4.107) 


one can reduce Eq. (4. 103) to (with zero 8e 0 and 8e f ): 


81 = 


P , 2 r f P rf ~ X f P xf 

hf v, 

h r 


8h f - P e 8ft 


(4.108) 


For some of the problems the final value of the angular momentum h will be left free, 
meaning that, for such problems one will always have the additional transversality 
condition: 


if hf is free, then: 


P hf + 


2 r f P rf L Q 

hr 


(4.109) 


Finally, for problems in which the orientation difference £2 between the initial and 
final conics is left free one will also have the transversality condition: 

if ft is free then P e = 0 (4. 1 10) 

At the limiting case in which at least one of the eccentricities associated with the initial 
or final conics tends to zero, this last category of problems includes the very important 
subclass of transfers from or to a circular orbit. For all the problems in this subclass the 
exact optimal value of P e is zero. 

In concluding this section one must note that the transversality conditions obtained 
here were all exact and also necessary for a (power-limited) trajectory to be optimal. 

4.6.4 Application of the Transversality Conditions 

Sections 4.6. 1 and 4.6.3 have supplied the exact first-order necessary conditions for 
optimality for the problem corresponding to power-limited orbital motion. The goal is of 
course to use the solution uncovered in Section 4.4 and try to satisfy as many first-order 
conditions as possible. Note that, along a trajectory violating even a single first-order 
condition, the first variation of the (augmented) performance index fails to be identically 
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zero. The hope therefore is that, if any first-order condition is violated by the solution 
uncovered in Section 4.4, the violation will be small, resulting in nearly-optimal trajectories 
(see also Section 4.6.5 for more on this subject). Thus, this section will combine the 
results of Sections 4.6.2 and 4.6.3 in order to uncover explicitly the form that the 
transversality conditions assume under the tangential thrust assumption. As seen in Section 
4.6.3, there will always be at least two transversality conditions, given in Eq. (4.106). 
Thus, substituting for the costates from Eq. (4.77) into Eq. (4.106) one obtains: 


f gX A 
l K j 


P e =0 ; 


r?Vp 
. hf 


P e =0 


(4.111) 


where V 0 , V f are the (nondimensional) speeds of the vehicle at the initial and final times 
respectively. Note that, in terms of the state variables h, x, and r, the (nondimensional) 
speed of the vehicle, V, is given by (see Fig. 4.1): 


V = 




(4.112) 


Since for obvious reasons the terms within the parentheses cannot be zero at the endpoints 
of a transfer, Eq. (4. 1 1 1) implies that the two transversality conditions in Eq. (4. 106) can 
be satisfied by the solution corresponding to the tangential thrust assumption only if one 
selects P e as: 


P e =0 


(4.113) 


This result was used as the justification in Section 4.4.3 for reducing Eq. (4.77) to 
Eq. (4.78). If one now substitutes for the costates from Eq. (4.78), and also takes into 
account Eq. (4.1 12), then, Eq. (4.109) assumes the form: 


A. 

hf 



= 0 


(4.114) 


But on the final conic, the following two well-known 34 relations of Keplerian motion are 
valid (see Appendix A): 

V f 2 - — = — — ; h f 2 =a f (l-e?) (4.115) 

r, a f 
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a, in Eq. (4.1 15) is the semimajor axis of the final conic. If the final conic is an ellipse, then 
a, is positive and e f is between 0 and 1. If the final conic is a hyperbola, then a, is negative 
and e f is greater than 1 . If the final conic is a parabola, then a, is infinite, e f is 1, and hf is 
finite 34 . Substituting from Eq. (4.1 15) into Eq. (4.1 14), Eq. (4.109) can now be written as: 

77( e f -l) = 0 (4.116) 

h f 

The throttling parameter A cannot be taken equal to zero, since in such a case the 
thrust is zero and there can be no transfer. Also, the angular momentum of the final conic h f 
is always finite. Accordingly, from Eq. (4.1 16), the only way in which the transversality 
condition given in Eq. (4.109) can be satisfied by the solution corresponding to the 
tangential thrust assumption is by selecting the final eccentricity e f as: 

e f = 1 (4.117) 

But problems for which e f =1 are just escape problems! Therefore, leaving the final value 
hf free, combined with the solution found using the tangential thrust assumption, suggests 
an escape problem! (from an arbitrary conic). 

This section now concludes by summarizing the important results from the present 
and the last section: 

(i) For problems in which the angular momentum at the initial and final times is left 
free the exact optimal (power-limited) cost is given by Eq. (4.105). 

(ii) For all the transfer problems in which either the orientation difference between 
the initial and final conics is left free or either one of the conics (or both) is circular the 
exact optimal value of P 6 is zero (see Eq. (4. 1 10)). 

(iii) For transfer problems between two arbitrary conics, in which the orientation 
difference between the conics is fixed, the costate P e is also zero, but as a consequence of 
the tangential thrust assumption (see Eq. (4.1 13)). 

(iv) Leaving the angular momentum at the final time free implies that the eccentricity 
of the final conic should be fixed at 1, and, fixing the eccentricity of the final conic at 1 
implies that the angular momentum at the final time may be left free (see Eq. (4. 117)). 

Note that (i) and (ii) are independent of the tangential thrust assumption, while (iii) 
and (iv) are consequences of the tangential thrust assumption. 
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4.6.5 How Good is the Tangential Thrust Assumption? 


One will sooner or later have to go back and check the validity of the tangential thrust 
assumption if one is to be able to say anything at all about the possible near-optimality of 
the thrust program proposed in Eq. (4.10). Recall that the tangential thrust assumption was 
used in three steps: The first two steps consisted of dropping out Eq. (4.56) representing 
the costate equation for the radial distance r, and changing the right-hand-side of Eq. (4.54) 
representing the costate equation for the angular momentum h. This led to the expressions 
for the costates P x , P h , and P r given in Eq. (4.77). Then, enforcing the transversality 
conditions (see Sections 4.6.3 and 4.6.4) "fixed" the value of the costate P 9 at zero and 
simplified the costate expressions for P x , P h , and P r as in Eq. (4.78). A preliminary check 
of the validity of the tangential thrust assumption should therefore consist of a comparison 
of the left and right-hand-sides of Eqs. (4.54) and (4.56), with P e = 0, and using the 
expressions for the costates P x , P h , and P r found in Eq. (4.78). Differentiating the 
expressions for the costates P x and P r in Eq. (4.78), and using Eqs. (4.78), and (4.49) 
through (4.5 1), one finds that, for P e = 0, Eq. (4.54) is always satisfied, while the left and 
right-hand sides (LHS, RHS) of Eq. (4.56) assume the form: 

LHS of Eq. (4.56): r + — = (4.118.a) 

hr r r hr 

O Ay lAhv A 2 

RHS of Eq. (4.56): + * (4.118.b) 

hr r r 

Thus, Eq. (4.56) is satisfied only when (A 2 /h 2 r 2 ) is zero. In other words, when (A 2 /h 2 r 2 ) is 
very nearly zero, the thrust program proposed in Eq. (4.10) is very nearly optimal (in the 
power-limited sense), because, in such a case, every single first-order necessary condition, 
except Eq. (4.56), is satisfied exactly, while Eq. (4.56) is satisfied very nearly. Note that 
the term -(A 2 /h 2 r 2 ) causing the disagreement in Eq. (4.118.a) is second order in the 
throttling parameter A compared to the remaining terms. Since the throttling parameter A is 
inversely proportional to the transfer duration (see Eq. (4.83)), for long duration transfers 
this term is negligible, the optimal thrust is approximately tangential, and the thrust 
program proposed in Eq. (4.10) is practically optimal! The result in Eq. (4.118) also 
suggests that the assumption that the optimal thrust is approximately tangent to the optimal 
flight path is more consistent with transfers that are performed far away from a planet and 
on which the average levels of the magnitude of the angular momentum remain high. This 
observation is clearly in agreement with one's intuition. 
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As the duration of transferring between two given conics tends to infinity, and the 
thrust levels go to zero (see Eq. (4.83)), the solution supplied in Section 4.4, 
corresponding to the thrust program proposed in Eq. (4.10), asymptotically tends to the 
optimal solution for power-limited transfers. In fact, in such a case, both the optimal 
solution and the solution of Section 4.4 tend to Keplerian motion, for which the cost is 
zero! Note that it has been well-established in the literature, that, for power-limited 
transfers between two given conics the optimal cost is inversely proportional to the transfer 
duration 44 - 51 . Assuming that the angular momentum h 0 , the eccentricity e 0 , and the 
orientation co 0 of the initial conic are always fixed, and that a transfer is performed using a 
one-segment transfer trajectory, that is, a trajectory on which the throttling parameter A has 
a single value, one can distinguish between the following three classes of transfers in 
conjunction with the near-optimality of the thrust program proposed in Eq. (4. 10) and the 
corresponding exact analytic solution given in Section 4.4 (see Sections 4.7 through 4.12 
for a justification): 

(i) For "pure" transfer problems the shape, size, and orientation of the final conic are 
always of interest. This means that the angular momentum h,, the eccentricity e f , and the 
orientation to f of the final conic are always fixed. For such problems, due to the loss of 
controllability suffered because of the tangential thrust assumption, the duration of (a one- 
segment) transfer cannot be preassigned. Solutions exist only for a finite set of values of 
the transfer duration. Accordingly, one is not free to perform the (one-segment) transfer in 
arbitrarily large time intervals, using arbitrarily small thrust levels. Once the initial and final 
conics are fixed, the minimum amount of fuel that is required to perform a one-segment 
transfer using the thrust program proposed in Eq. (4.10) is also fixed. 

(ii) An important exception in case (i) corresponds to the subclass of transfers 
between two coplanar circular orbits. This turns out to be an interesting singular case (see 
Section 4.8). If both the initial and final conics are circles, then, one-segment transfer 
solutions between the conics do not exist. However, if the eccentricity of the final conic is 
slightly off zero, then there exists a very large number of solutions! Specifically, as the 
eccentricity of the final conic tends to zero and the final conic tends to a perfect circle, the 
number of possible solutions increases without bound, and one can in such a case find 
solutions corresponding to arbitrarily large transfer durations! Thus, in this case the amount 
of fuel required for a transfer can (in principle) be reduced by as much as desired (within 
the operating limits of the power-limited propulsion system), as long as one picks the 
solution corresponding to a "large enough" transfer duration. 
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(iii) Escape problems are a special kind of transfer problem in which the interest is 
in escaping the gravitational field of a planet. Accordingly, the size of the final parabolic 
trajectory is of little importance and one can now get away (in practice ) by specifying only 
the eccentricity (e f = 1) of the final parabolic trajectory and by leaving its angular 
momentum h f free. Because of this, one can now freely preassign the duration of 
transferring to tlje final parabola, and thus, for such problems, one can (in principle) make 
the amount of the required fuel very small by picking a large enough transfer duration. A 
variation on these problems is the one for which the angular momentum of the final 
parabola is specified, but its orientation (0 f is left free. This situation is not as practical as 
the previous one because there is usually a need to specify the orientation co f and exit (for 
example) a planet's sphere of influence with a heliocentric speed having a preassigned 
direction. A very interesting property of the thrust program proposed in Section 4.3 is the 
fact that the corresponding exact analytic solution of Section 4.4 satisfies the corresponding 
transversality condition for both of the above cases. This was shown in some detail in 
Section 4.6.4. 

It will be shown in Section 4.12, that the difficulty in preassigning large transfer 
durations, associated with one-segment transfers, can be bypassed by allowing multiple- 
segment transfers, that is, transfers composed of a finite number of segments, on each one 
of which the throttling parameter A assumes a specific value. If multiple-segment transfers 
are allowed, then one can in principle still use the thrust program proposed in Eq. (4.10) 
and preassign arbitrarily large durations for each individual segment of a transfer between 
any two given coplanar conics! This result is of great practical importance, because it 
implies that if one is willing to compromise with respect to the transfer duration and "wait 
longer" one can still reduce the fuel consumption to very small levels using a power-limited 
propulsion system and the guidance corresponding to the exact analytic solution of Section 
4.4 (see Section 4.12 for more details). 


4.7 Boundary Conditions 

The boundary conditions for any transfer problem arise from the requirement that the 
state of the vehicle be continuous at the initial and final times. A "split second" before the 
initial time x = 0 the vehicle coasts along a Keplerian conic, specified by its orbital elements 
ho, e 0 , and co 0 . A "split" second after x = 0 the thrust is turned on, and the vehicle is on the 


76 



transfer trajectory, specified by its integration constants A, B, and C. Similarly, a "split 
second" before the final time x = x f the thrust is turned off, and a "split" second after x = x f 
the vehicle coasts along a Keplerian conic, specified by its orbital elements h f , e f , and co f . 
Combining Eqs. (4.16), (4.17), describing the transfer trajectory, and Eq. (A. 5), 
describing a Keplerian conic (see Appendix A), one finds that the following conditions 
should be satisfied at t = 0 and x = x { : 


*o = 


ho 


l + e 0 cos(6 0 -co 0 ) 1 + Bcos(6 0 -C) 


(4.119) 


x o = 


_ e 0 sin(9 0 -co 0 ) _ 2r 0 A + Bsin(Q 0 -C) 


(4.120) 


r f = 


hf 


hf 


l + e f cos(0 f -(o f ) 1 + Bcos(0 f -C) 


(4.121) 


e f sin(0 f -to f ) _ 2r f A + Bsin(0 f -C) 
h f h f 


(4.122) 


In addition, one always has the explicit expression for the throttling parameter A: 


A _ hf — ho. (4.123) 

'tf 

found in Eq. (4.83) from the requirement that h=h 0 at x=0, and h=hf at x=x { . 

The final condition arises from the requirement that at the final time t=x f , one among 
Eqs. (4.25), (4.26), (4.34), and (4.42) of Section 4.5, relating h with 0 along the transfer 
trajectory, should be satisfied. Only one among these equations applies along a particular 
transfer, according to whether B=0, B=l, 0<B<1, or B>1 respectively. Note that (see 
Section 4.5) there is a one-to-one correspondence between the generalized eccentric or 
hyperbolic anomalies E or H, and the argument of latitude 0 or £, where £=0-C. Note also, 
that the generalized hyperbolic anomaly has been denoted by H, since there is no danger of 
confusion with the Hamiltonian H given in Eq. (4.44) which does not play an important 
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role in this chapter. The conditions given in Eqs. (4.1 19) through (4.123) can only specify 
the final values E f , or 0 f up to a multiple of 2n, that is, one can write: 


E f = Efp + 2kn ; 0 f = ©„> + 2k7t 


(4.124) 


where k=0,l,2,3,4,...etc, and E^ and are between 0 and 2n. Using Eq. (4.124), one 
can write Eqs. (4.25), (4.26), (4.34), and (4.42) of Section 4.5 at the final time X=t f in a 
unified fashion as: 


(i) For B=0 


T f (h f + h 0 ) 
4nh 2 0 h 2 f 


-^(<V- e o) = k 


(4. 125. a) 


(ii) For 0<B<1 


T f (h f +h 0 )b | B(sinE fP -sinE 0 )-(E (P -E 0 ) _^ 
47thohf 2 tc 


(4.125.b) 


(iii) For B= 1 


x f (h f + h 0 ) (y?-yo) + 3(y f -y 0 ) 

47thohf 1271 


(4.125.C) 


(iv) For B>1 


x f (h f + h 0 )b B(sin H f - sin H 0 ) - (H f - H 0 ) _ Q 

47thohf 2tt 


(4.125.d) 


where, b in Eqs. (4.125.b) and (4.125.d) is defined respectively as: 

b = (l — B 2 ) 3/2 or b = (B 2 - 1) 3 ' 2 


and y f , y 0 in Eq. (4.125.C) are defined as: 


y f 


= tan 


rtv 

l 2 )' 


y 0 


= tan 



(4.126) 


(4.127) 


Note that k in Eqs. (4. 125. a), (4.125.b) is a nonnegative integer, related directly to 
the number of revolutions around the planet during the transfer. For the case in which B is 
less than one, among all the particular solutions satisfying Eqs. (4.119) through (4.123) 
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only the ones for which the left-hand-side of Eqs. (4.125.a) or (4.125.b) is a nonnegative 
integer are acceptable. It is clear from Eqs. (4.125) that the time of transfer T f and the 

number of revolutions around the planet during the transfer are intimately connected with 
each other. 


4.8 Transfer from a Circular Orbit to an Arbitrary Conic 

For this class of problems the quantities that ar t fixed at the initial and final times are 
ho, e 0 , hf, e f , and co f (co 0 is not defined). The initial conic is assumed to be a circle, meaning 
that the initial eccentricity is zero. With no loss of generality, due to the circular symmetry 
of the initial conic, one can assume that the orientation of the final conic is zero. It will also 
be assumed that is different than ho- Explicitly, one has: 

e 0 =0; (0 f = 0 ; h f *h 0 (4.128) 

Taking the above into account, and according to whether hf is greater or less than h 0 , Eqs. 
(4.119), (4.120) yield: 


A > 0: r 0 = hg ; 0 O -C = -| ; B = 2h’A (4.129.a) 

A < 0: r 0 = h^ ; 0 O -C = | ; B = -2h*A (4.129.b) 

Using Eqs. (4.129) into Eqs. (4.121), (4.122) one obtains, for both positive or negative 
A, the conditions: 


e f cos0 f =2hoAsin(0 f -0 O ) 


(4.130) 


e f sin 0 f = 2 r f A - 2hg A cos(0 f - 0 O ) (4.131) 

Note that a solution does not exist if e,=0, that is, if the final conic is also a circular 
orbit. It will therefore be assumed that e f >0. The singular case e,=0 will have to be treated 
as a limiting case in which e f is taken as very small. Eqs. (4.130), (4.131) can be combined 
and solved for the sine and the cosine of 0 f : 
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(4.132) 


2A(hf -hoe f cos0 o )(e f +2hoAsin0 o ) 

sin0 f = — — ■ -p - 

(e f + 2hoAsin0 o ) + 4hoA 2 cos 2 0 o 


2A(h 2 - hj;e f cos0 o )(211qAcos 0 o ) 

(e f + 2hoAsin0 o ) 2 + 4hgA 2 cos 2 0 o 

Using now the trigonometric identity: 

sin 2 0 f +cos 2 0 f = 1 

and substituting from Eqs. (4.132), (4.133) into Eq. (4.134) one obtains a condition for 
the starting value of the argument of latitude 0 O on the initial circular orbit. Using standard 
trigonometric identities 33 , this condition can be put in the form of a fourth order polynomial 
equation: 

b 4 y 4 + b 3 y 3 + b 2 y 2 +b,y + b 0 = 0 (4.135) 


(4.133) 


(4.134) 


where the unknown y is defined as: 


y 



(4.136) 


The coefficients of Eq. (4.135) are given explicitly in Appendix D. Note that, due to the 
identity expressed in Eq. (4.134), the expressions for the sine and the cosine of 0 f given in 
Eqs. (4.132) and (4.133) can be considerably simplified and written as: 


hf -h;5e f cos0 o 


(4.137) 


COS0 f = 


hoCos0 o 
h 2 -hoe f cos0 o 


(4.138) 


The solution methodology for a given problem proceeds in the following fashion: 
One assumes a value for the final time x f and evaluates A from Eq. (4.123), and B from Eq. 
(4.129). B cannot be zero. Then, solving the polynomial equation (4.135) one finds all the 
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real roots y, and determines the corresponding starting points 0 O on the initial circular orbit. 
Any such 0 O can be taken (at first) with no loss of generality to be between 0 and 2k, that 
is, O<0 o <2tt. For every such 0 O , Eqs. (4.137), (4.138) fix the corresponding 0 f up to a 
multiple of 2k. That is, 0 f = 0fp+2k7C, with O<0 fP <2jt, and k=0,l,2,3,...etc. If B>1, then 
k=0 and one checks whether the left-hand-side of Eq. (4.125.C or d) is zero. If B<1, then 
one checks whether the left-hand-side of Eq. (4.125.b) is equal to zero or a positive 
integer. If not, one repeats the procedure with a new value of x f . Therefore, one way of 
finding all the possible solutions for a given problem is to plot the left-hand-side of Eq. 
(4.125) as a function of x f , and for all the real roots of Eq. (4.135). The values of x f (for 
0<B<1) for which this plot intersects any positive integer (k=l,2,3,...), plus the values of 
x f (for B>0) for which it intersects zero (k=0) represent the candidate solutions to a given 
problem. For k=0 the starting point on the initial circular orbit is either 0 O , or 0 o -27t. Each 
such candidate value of x f is a real solution only if the transition from the starting point to 
the final point can be made without going through a value of 0 at which the angular 
momentum h, evaluated through the applicable one among Eqs. (4.25), (4.26), (4.34), and 
(4.42), becomes infinite. For each such value of x f representing a real solution the 
corresponding value of k fixes both the value of 0 f and the number of full revolutions 
around the planet during the transfer. The transfer problem is then completely solved, that 
is, everything related to the transfer can be easily calculated. 

One may argue that the problem with the above procedure is that, to find all the 
possible solutions, one will have to plot Eq. (4.125) for x f between 0 and °°! However, 
after working on specific problems for a while, it becomes apparent that, except for the 
case of transfers between a circular and a very nearly circular orbit, the interval of values of 
x f for which Eq. (4.135) has real roots is very limited. However, it has not been proven in 
this chapter that solutions for much higher values of x f do not exist, and this topic may 
deserve some further consideration in the future. An example of a transfer between a 
circular and a very nearly circular (e^O.05) orbit is given in Figs. 4.6 and 4.7. From Fig. 
4.6 it is seen that there are 9 solutions for this example. The actual transfer (Eq. (4.17)) 
corresponding to the solution with k=10 is plotted in Fig. 4.7. It is interesting to remark 
here that if one starts reducing the eccentricity of the final orbit toward zero, then, both the 
number of solutions and the corresponding transfer durations greatly increase. This fact is 
vividly depicted in Figs. 4.8 through 4.11, for which e f is equal respectively to 0.01, 
0.005, 0.001, and 0.0005. Since hf is still equal to 2, Figs. 4.8 through 4.11 basically 
correspond to (practically) the same transfer problem defined in Fig. 4.6. Thus, in practice, 
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for the case of transfers between two coplanar circular orbits, by playing around with the 
(very small) value of e f one can find (one-segment) transfer solutions of arbitrarily large 
transfer duration. This result has important practical implications from the point of view of 
power-limited optimality and fuel consumption and more will be said on it in Section 4.12. 
The case of transferring between two coplanar circular orbits, along with a complete 
mathematical description of what happens as e f tends to zero, is given in Appendix D. Note 
that when e f = 0 Eq. (4.135) has no real solution. In fact, when e f is exactly equal to zero, 
to satisfy this polynomial equation one must perform the transfer in infinite time using zero 
thrust. The case with e f = 0 is thus clearly a singular case. Figures 4.12 through 4.17 
supply three examples of transfers (Eq. (4.17)) to elliptic orbits of low and high 
eccentricities and to a hyperbolic conic. It appears that, as the eccentricity of the final conic 
increases, the number of (one-segment) solutions decreases (very fast), but does not go to 
zero. As long as the final conic does not intersect the initial conic (see Section 4.1 1) there 
always seems to be at least one (one-segment) solution. 


4.9 Escape from a Circular Orbit 

Escape problems are just a special case of transfer problem, in which the primary aim 
is to escape from the gravitational attraction of a planet 34 . This can be done by transferring 
to an open conic for which the eccentricity e f is greater than or equal to one 34 . For e f >l the 
vehicle escapes with finite speed at infinity 34 . For e,=l the vehicle escapes with zero speed 
at infinity 34 . Accordingly, from an energy point of view, the cheapest way to escape is by 
transferring to a conic for which e,=l . The actual final value h, of the angular momentum of 
the vehicle is not important and can be left free. Then, as was seen in Section 4.6.4, as 
long as epl, the corresponding transversality condition (Eqs. (4.109), (4.116)) is 
automatically satisfied. Thus, for an escape problem from a circular orbit, it will be 
assumed that the quantities that are fixed at the initial and final times are hg, e 0 , e,= 1 , and 
(OpO; gj 0 is not defined, while h f is left free. The important difference between such 
problems and the transfer problems considered in Section 4.8 is that now, solutions may 
exist (in principle) for any value of T f . This is so, because the free final value of the angular 
momentum hf now plays (technically) the same role that the time of transfer x f played in 
Section 4.8. Physically, escape can be achieved at any desired final time by appropriately 
selecting the constant A, so that the mechanical energy (per unit mass) of the system is 
increased at the correct (mean) rate. Thus, the solution methodology for a given escape 
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problem is basically identical to the one given in Section 4.8. The difference is that instead 
of plotting the left-hand-side of Eq. (4.125) as a function of x f , one now plots it as a 
function of h,. Figure 4. 18 supplies the plot of the left-hand-side of Eq. (4.125) for such an 
escape problem. The transfer trajectory (Eq. (4.17)) corresponding to point S in Fig. 4.18, 

with k— 3, is plotted in Fig. 4.19, along with the initial (circular) and the final (parabolic) 
conic. 


4.10 Transfer Between two Arbitrary Conics, 
and Escape from an Elliptic Orbit 

The more general problems relating to transfers between two arbitrary conics and to 
escape from an elliptic orbit can be treated basically in exactly the same way as the 
corresponding problems of Sections 4.8 and 4.9. However, the corresponding algebra is 
considerably more involved. 

For this class of problems the quantities that are always fixed at t=0 and at T=t f are 
respectively h 0 , e 0 , co 0 , and e f . For transfer problems h, and to f are also fixed, and (one- 
segment) solutions exist only for particular values of x f . For escape problems the initial 
conic is assumed to be an ellipse (e 0 <l), h f is free, e^= 1 , and (one-segment) solutions (in 
principle) exist for any value of x f . With no loss of generality one can always assume that 
© 0 =0, that is, the orientation of the initial conic is zero. FromEqs. (4.1 19) through (4.122) 
one then obtains: 


e o cos0 o = Bcos(0 o - C) 

(4.139) 

e o sin0 o = 2r o A + Bsin(0 o -C) 

(4.140) 

e f cos(0 f - co f ) = Bcos(0 f - C) 

(4.141) 

e f sin(0 f - co f ) = 2r f A + Bsin(0 f - C) 

(4.142) 


From Eqs. (4.139) and (4.141) B can be expressed as: 
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(4.143) 


e 0 cos8 0 _ e f cos(8 f -co f ) 
cos(6 0 - C) cos(8 f - C) 


With no loss of generality one may assume that 0<C<27t. The equality on the right in Eq. 
(4.143) can be solved for the tangent of C to yield: 


tanC = 


e 0 cos8 0 cos8 f - e f cos(8 f - © f )cos6 0 
e f cos(8 f -co f )sin8 0 - e 0 cos8 0 sin8 f 


(4.144) 


Equation (4. 144) supplies two roots for C in 0<C<2jt. The possible solution(s) for C are 
the one(s) for which B is nonnegative. The remaining steps are algebraically quite intense. 
First, substituting from Eqs. (4.143), (4.144) into Eqs. (4.139)-(4.142) one can solve for 
the sine and the cosine of 8 f as: 

sin8 f =— ; cos8 f =— (4.145) 

n 3 n 3 

where the quantities N,, N 2 , and N 3 are given by: 

Nj = 2hoAsin8 0 + (l + e 0 cos8 0 )(e f cos(O f -e 0 ) (4.146) 


N 2 =2hoAcos8 0 -(l + e 0 cos8 0 )e f sinco f 


(4.147) 


N 3 = (2hfA + e 0 e f sinco f )(l + e 0 cos8 0 )-2hoAe f cos(8 0 -co f ) (4.148) 

Then, using the trigonometric identity: 

sin 2 6 f + cos 2 8 f = 1 (4.149) 

and substituting from Eqs. (4. 145)-(4. 148) into Eq. (4.149) results in a condition for the 
starting point 6 0 on the initial conic. Using standard trigonometric identities 33 , this condition 
can again (see Eq. (4. 135)) be put in the form of a fourth order polynomial equation: 

b 4 y 4 + b 3 y 3 + b 2 y 2 + b,y + b 0 = 0 
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(4.150) 



where the unknown y is (see Eq. (4.136)) defined as: 


y = tan 



(4.151) 


The coefficients of Eq. (4.150) are given explicitly in Appendix E. 

Note that for e 0 = co f = 0 the results of the present section reduce to the corresponding 
results of Section 4.8, valid for transfers and escapes from a circular orbit. 

The solution methodology for a given transfer or escape problem is more or less 
identical to the one given at the end of Section 4.8, with the analogies involved being too 
obvious to go through in detail once again. The left-hand-side of Eq. (4.125) can always be 
plotted as a function of x f or h f to determine the candidate solutions. However, due to the 
well-defined orientation difference between the initial and final conics, there arises in this 
case another possibility for escape problems that is worth exploring. If one chooses to 
leave this orientation difference free, then one may consider escape problems with fixed T f 
and hj. Recall that in such a case the exact optimal value of P e is zero (see Eq. (4.1 10)), 
that is, P e is not zero just because of the tangential thrust assumption. For escape problems 
in which co f is left free the candidate solutions are determined by plotting the left-hand-side 
of Eq. (4.125) as a function of co f , and only for 0<to f <2re. 

A transfer example between two elliptic orbits, involving no orientation change, is 
supplied in Figs. 4.20, 4.21. One must note here that for transfers between two coplanar 
elliptic orbits having the same orientation the number of (one-segment) solutions and the 
corresponding transfer durations greatly increase as the final eccentricity e f approaches the 
initial eccentricity e 0 . This situation is similar to the one corresponding to transfers between 
two coplanar circular orbits (see Section 4.8), and is depicted clearly in Figs. 4.22, 4.23. 
The solution search conducted in Figs. 4.22, 4.23 corresponds practically to the same 
transfer problem. An actual transfer corresponding to either Fig. 4.22 or Fig. 4.23 is not 
shown. Figures 4.24, 4.25 supply a transfer example between two elliptic orbits, involving 
a large orientation change. Figures 4.26, 4.27 supply a transfer example from an elliptic 
orbit to a hyperbolic trajectory, while Figs. 4.28, 4.29 supply a transfer example between 
two hyperbolic trajectories. The jump of magnitude 1.0 that appears whenever the left- 
hand-side of Eq. (4.125) is plotted is due to the fact that all such plots are done with 
O^6 0 < 2^, and O<0fp<27t. A maneuver such as the one given in Fig. 4.29 could play an 
important role during a planetary encounter, in which, having control over the vector of 
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heliocentric velocity of a vehicle at the exit of a planet's sphere of influence (Jupiter for 
example) is most desirable for setting up the next (planetary) encounter (with Saturn, 
Uranus, Neptune, etc.). 

Figures 4.30 and 4.31 supply two solutions for an escape problem from an elliptic 
orbit for which the orientation of the Final parabolic trajectory is left free. The solution that 
is cheaper, both in the power-limited and the constant ejection velocity sense, is the one 
given in Fig. 4.30, corresponding to (0 f « 0 (see Table 4.1, Section 4.13). Note that the 
plot of the left-hand-side of Eq. (4.125) as a function of 0) f is not shown for this problem. 

This section concludes on the note that, as long as the final conic does not intersect 
the initial conic, there always seems to be at least one (one-segment) solution for the 
corresponding (transfer or escape) problem. 


4.11 Existence of One-Segment Solutions, 
and some Geometric Considerations 

The search in the previous three sections was for one-segment solutions to given 
transfer or escape problems, that is, trajectories on which the throttling parameter A had a 
single value. An n-segment solution (n=l,2,3,...) will henceforth denote a trajectory on 
which the value of A switches n-1 times in 0<T<x f . Note that there is really no reason to 
refer to escape problems explicitly, since escape problems are just a special kind of transfer 
problem. Thus, for simplicity, one-segment solutions will henceforth be referred to just as 
"transfers", or as "one-segment transfers". Before discussing the existence of such 
transfers it will be helpful to adopt some notational conventions. For any transfer problem 
the word "from" will always be associated with the initial conic (at the initial time) and the 
word "to" with the final conic (at the final time). The quadruplet (h,e,a>,0) will be used to 
denote a (Keplerian) conic with orbital constants h, e, and to, with 0 fixing the point of 
departure "from" or arrival "to" the conic. If the conic is a circle the second entry will be 
zero and the third entry will be a dash. Sometimes the fourth entry may be omitted and a 
conic be denoted by the triplet (h,e,fo). Similarly, a one-segment transfer with constants 
A *0, B, C, and time duration T f > 0, will be denoted by the quadruplet (A,B,C,x f ), or by 

the triplet (A,B,C). 

Let now P be any point on a transfer trajectory (A,B,C), and let h, x, r, 0 be the state 
components at P. If the thrust were suddenly turned off at P, the motion that would follow 
after P would be Keplerian. The orbital elements of this Keplerian motion define the so- 
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called 34 instantaneous conic (h,e,co) at P. One can also define the associated conic at P to be 
the conic (h,B,C). The associated conic (h^B.C) at x=0 will be called the initial associated 
conic and the associated conic (h f ,B,C) at x=x f will be called the final associated conic. 
Figures 4.32 and 4.33 show the initial (ACO) and final (ACf) associated conics for the 
transfer examples given in Figs. 4.21 and 4.25. 

It is now possible to supply a few lemmas, related to the geometry of motion and the 
existence of one-segment transfers. It will henceforth be tacitly assumed that the motion 
along any transfer is performed under the thrust program proposed in this chapter (Eq. 
(4.10)), and that the initial and final conics are different. A detailed proof for a lemma will 
be supplied only when absolutely necessary. The proof of the following four lemmas is 
elementary. 

Lemma 4.1. A one-segment transfer trajectory (A,B,C,x f ) can never intersect 
itself. 

Lemma 4.2. If (A,B,C,x f ) is a transfer from conic (h o ,e o ,(O o ,0 o ) to conic 
(h f ,e f ,co f ,0 f ), then (A,B,C,x f ) is also a transfer from conic (-h f ,e f ,to f ,0 f ) to conic 
(~ho,e o ,co o ,0 o )- 

Lemma 4.3. If a vehicle takes off at x = 0 from a point K of a circular orbit 
(1^,0,-, 0 O ) along a one-segment transfer trajectory (A > 0,B,C,x f ), and if at time x > 0 the 
vehicle is at point L, then, the distance KL is equal to (e/2A), where e is the eccentricity of 
the instantaneous conic at L (see Fig. 4.34). 

Lemma 4.4. Consider a vehicle that takes off at x = 0 from a conic (h o ,e o ,co o ,0 o ) 
along a one-segment transfer trajectory (A > 0,B,C,x f ), and draw a vector with origin at the 
center of the planet, of magnitude (B/2A) and direction C-(nI2). Let K be the endpoint of 
this vector. Then, 

(i) If M is the point of departure of the vehicle from conic (h^e^CD^©^, then the 
distance MK is equal to (e(/2A) (see Fig. 4.35). 

(ii) If at time x > 0 the vehicle is at point L, then the distance KL is equal to (e/2 A), 
where e is the eccentricity of the instantaneous conic at L (see Fig. 4.35). 

The proof for Lemma 4.1 is obtained by considering the trajectory equation (Eq. 
(4.17)). The proof for Lemma 4.2 is obtained by considering Eqs. (4.1 19)-(4.122). The 
proof for Lemmas 4.3 and 4.4 is obtained by considering Eqs. (4.130), (4.131), (4.139)- 
(4.142), dropping the subscript f from Eqs. (4.141), (4.142), and using the law of 
cosines 33 . 
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Lemma 4.5. A one-segment transfer (A,B,C,x f ) from conic (h^e^©^) to conic 
(h^e^CDpOf) does not exist if at least one among the following two equalities is valid: 


h 0 = h f ; 


hf 

l" e f 


(4.152) 


Proof 4.5. If the equality on the left is valid then A=0 and the motion is Keplerian. 
If the equality on the right is valid then the orbital energy of the two conics is the same 34 
(see Eqs. (A.7), (A.9)), contradicting the fact that for nonzero A the orbital energy either 
increases or decreases (strictly monotonically) on a one-segment trajectory (Q.E.D.). 

Note that if at least one among the equalities given in Eq. (4.152) is true, then the 
initial and final conics intersect. This fact, combined with one's previous experience with 
specific examples, leads one to suspect that Lemma 4.5 is most probably a particular case 
of a more general result: One-segment solutions actually don t exist when the initial and 
final conics intersect. For example, consider the problem of transfer between two 
hyperbolic conics (see Fig. 4.36) having the same eccentricities, different (but both 
positive) angular momenta, and opposite orientations (co f = © 0 - k). It appears that it is 
impossible to construct a one-segment transfer trajectory that joins these two conics. Due to 
the tangential character of the thrust, the vehicle, approaching from leg a of the initial 
hyperbola cannot transfer to leg d of the final hyperbola unless it executes at least one full 
revolution around the planet. But this implies an initial capture, energy loss, and a negative 
value of A. The subsequent transfer to leg d of the final hyperbola implies a final escape, 
energy gain, and a positive value of A. 

Lemma 4.6. The instantaneous conic (h,e,©) and the associated conic (h,B,C) (at 
the same instant) have always two and only two points in common. At any instant of time 
along a transfer (A,B,C,x f ) the vehicle occupies one of those two points and has velocity 

that is tangential to the instantaneous conic. 

Proof 4.6. (h,e,co) and (h,B,C) are two conics having the same angular momentum 

but different eccentricities and orientations (Q.E.D.). 

The following lemma can be proved by considering Eq. (4.17). 

Lemma 4.7. A one-segment transfer (A,B,C,t f ) from conic (h^e^©^,,) to conic 
(h^epW^Qf) has no common points with either the initial, (h^B.C), or the final, (h^B.C), 
associated conic for 0 < t < T f . 

A fact that is intuitively obvious, but the proof of which is not trivial, is provided by 
the lemma that follows. 
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Lemma 4.8. Let (A,B,C,x f ) be a one-segment transfer from conic (h^e^co^©^ to 
conic (h f ,e f ,(O f ,0 f ), and let x, and x,+dx be two infinitesimally separated time instants in 
0<x<x f . Then, the instantaneous conic (hpe,,©,,©,), and the instantaneous conic 
(h,+Ah,e,+Ae,co,+Ato,0,+A0), associated with the transfer trajectory at the instants x, and 
x,+dx respectively, can have at most one common point. 

Proof 4.8. The trajectory equations corresponding to the two instantaneous conics 
are given by: 


r = 

1 +e,cos(0-co ] ) 

(4.153) 

r _ (h, + ih) 2 

1 + (e, + Ae)cos(0 — CO, — Ato) 

(4.154) 

To first order in dx the changes in the orbital elements h,, e,, and to, 

can be found as 34 : 

Ah = Adx 

(4.155) 

Ac _ 2A h +cos ( 0 i-®.)] dT 

h, 

(4.156) 

Affl= 2Asm( e l - a ,|)^ 

e ihj 

(4.157) 

Equating the right-hand-sides of Eqs. (4.153), (4.154), using Eqs. (4.155)-(4.157), 

and expanding in Taylor series the terms sin(Ato) and cos(Aoa) one finds that to first order 

in dx the following equality must be valid: 


[cos(0 - 0, ) - l]dx = 0 

(4.158) 


for arbitrary, but infinitesimally small dx. This implies that the two conics can have at most 
one common point, corresponding to 0=0, in O<0<2rc just as the lemma claims (Q.E.D.). 

Along a one-segment transfer trajectory the mechanical energy and angular 
momentum (per unit mass) associated with the instantaneous conic either increase (A > 0) 
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or decrease (A < 0) strictly monotonically with time. Thus, Lemma 4.8 implies the lemma 
that follows, which is basically a companion to Lemma 4.7. 

Lemma 4.9. A one-segment transfer (A,B,C,x f ) from conic (ho,e o ,© o ,0 o ) to conic 
(h f ,e f ,co f ,0 f ) has no common points with either (ho,e o ,© o ,0 o ) or (h f ,e f ,CO f ,0 f ) for 0 < x < x f . 

Lemmas 4.7 and 4.8 suggest that for A > 0 a one-segment transfer trajectory always 
stays "outside" the initial and the initial associated conics, and "inside" the final and the 
final associated conics. Analogous observations can be made also for A < 0. When the 
initial and final conics do not intersect the allowable space for the transfer trajectory has the 
shape of a "ring" (see Figs. 4.32, 4.33). However, as the initial and final conics come 
closer and closer toward intersecting each other, this allowable space is constrained 
considerably (see Fig. 4.33), and after the intersection it is cut-off, that is, it looses its 
"ring" character. 

The above considerations, and in particular Lemmas 4.7, 4.8, and 4.9, point toward 
the following underlying theorem: 

Theorem 4.1. Consider the problem of transferring from conic (ho,e 0 ,a) 0 ,6 0 ) to 
conic (hpepCDpOf) using the thrust program proposed in Eq. (4.10) and the corresponding 

exact analytic solution of Section 4.4. Then, 

(i) If (h^e^CDo^o) and (hf,e f ,© f ,0 f ) have more than one common point (that is, if 

they intersect), then a one-segment transfer solution does not exist. 

(ii) If (h^eo,©,,,©^ and (h f ,e f ,a> f ,0 f ) have only one common point (a point of 
tangency) then there exists at most a single impulsive one-segment transfer solution, 
performed at the point of tangency of the two conics, for which A is infinitely large, x f is 
infinitely small, and Ax f is finite and equal to hf-h^. 

The proof of Theorem 4.1 is rather lengthy and is given in Appendix F without using 
any of Lemmas 4.7, 4.8, or 4.9. In fact, the proof given in Appendix F can also be 
considered to be an indirect proof for Lemmas 4.8 and 4.9. 

Before closing this section one can also state the following conjecture, that the author 
feels is true, having to do with the existence of one-segment and two-segment solutions. 

Conjecture. Consider the problem of transferring from conic (h o ,e o ,© o ,0 o ) to conic 
(h,,e f ,© f ,0 f ) using the thrust program proposed in Eq. (4.10) and the corresponding exact 
analytic solution of Section 4.4. Then, 

(i) A two-segment transfer solution always exists. 

(ii) A one-segment transfer solution (A,B,C,X f ) with finite A and x f > 0 exists if and 
only if (h^e^©^©^ and (h^e,.©,,©,) have no common point. 
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Note that in case (i) above the selection of the intermediate instantaneous conic at the 
point where the switching of the value of A occurs is not unique but will most probably 
have to be decided on the basis of optimality. Figures 4.37 and 4.38 supply two examples 
of two and three-segment transfers respectively. Figure 4.37 corresponds to a situation 
similar to the one given in Fig. 4.36. In Fig. 4.38 there is an intermediate transfer segment, 
KS, which is shown with dotted line, which corresponds to zero thrust (A=0) and 
Keplerian motion along a circular orbit. 


4.12 When is it Possible to Preassign Arbitrarily Large Transfer Durations? 


If multiple-segment transfers are allowed, the answer to the above question is: 
Always. Before expanding on this issue however, and explaining its practical importance, 
it is insightful to make a few comments concerning the behavior of the thrust acceleration 
along a one-segment transfer. Recall that the angular momentum of the vehicle is given by 
h = rVcosy, where, V is the speed of the vehicle, given by Eq. (A.8) in Appendix A. Using 
the components given in Eq. (4.10), the thrust acceleration corresponding to the proposed 
(tangential) thrust program can be written as: 



A 

rcosy 


(4.159) 


Squaring and differentiating the right-hand-side of Eq. (4.159), and then using Eqs. (4.1 1) 
through (4.14), the first and second time derivatives of (e 2 /2) can be found as: 


djvy a 2 x 

dx ^ 2 j h 2 r 2 



(4.160) 


where, the second derivative was evaluated at the points where x is zero. From Eq. (4.160) 
one can conclude that the thrust acceleration has a (local) maximum at the successive 
perigees and a local minimum at the successive apogees of an (osculating) one-segment 
transfer trajectory. Figures 4.39 and 4.40 supply the variation of (e 2 /2) with x and 0 
respectively for the transfer example of Fig. 4.2. Note how the variation in Fig. 4.39 
resembles a series of impulses delivered at the successive perigees of the (osculating) 
transfer trajectory. 

Using Eq. (159), and the expression for A found in Eq. (4.83), the power-limited 
cost associated with a one-segment transfer can be written as: 
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'PL 




- a2 f dT _ ( h f ~ h o ) 2 f 


2 J r 2 cos 2 y 2 t 2 J r 2 cos 2 y 


(4.161) 


This expression can be further manipulated into: 


T - (hf h o) 

J PL — o _ 


2x f 


average value of 


2 2 

r cos y 


along the transfer 


(4.162) 


The average value of (l/r 2 cos 2 y) along a one-segment transfer remains a well-behaved 
quantity as the duration of transfer x f tends to infinity. Thus, as x f tends to infinity, the 
power-limited cost associated with a one-segment transfer tends to zero. This behavior 
compares well with the behavior of the optimal cost, since it is well-known 42 51 , that the 
optimal power-limited cost for a given transfer problem is inversely proportional to the 
transfer duration, and tends to zero as the transfer duration tends to infinity. For a multiple- 
segment transfer the power-limited cost is given by a summation of terms similar to the one 
appearing on the right-hand-side of Eq. (4.162). Thus, for a multiple-segment transfer, if it 
is still possible to preassign large durations for each individual segment, then, the total fuel 
consumption can still be made very small, but at the expense of waiting longer for 
completing the (overall) transfer. Thus, whether it is possible to preassign arbitrarily large 
transfer durations or not is an important practical issue. By recalling the results of Sections 
4.6.5, and 4.8 through 4.10, one can now distinguish between the following three cases: 

(i) Escape problems: For all such problems it is possible to use a one-segment 
transfer trajectory. Because the final value of the angular momentum hf is left free, the time 
duration x f for transferring to the final parabolic trajectory can in principle be freely 
preassigned. As this time duration tends to infinity the one-segment transfer and the 
corresponding power-limited cost asymptotically tend to the optimal transfer and the 
optimal power-limited cost respectively. 

(ii) Transfer between two coplanar circular orbits: For this problem it is 
still possible to use a one-segment transfer trajectory. A small transfer duration x f in this 
case does not exist. However, a large transfer duration x f can still be preassigned (see 
Section 4.8 and Appendix D). Again, as this time duration tends to infinity this one- 
segment transfer and the corresponding power-limited cost asymptotically tend to the 
optimal transfer and the optimal power-limited cost respectively. 
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(iii) Transfer between two arbitrary conics: For this problem it is not in 
general possible to use a one-segment transfer trajectory. However, it is in principle 
possible to perform the overall transfer using a multiple-segment transfer trajectory, for 
which one can still preassign arbitrarily large durations for each segment. For example, a 
vehicle can first transfer to an intermediate closed orbit in time x n , and then, after coasting 
for some time x n on this intermediate closed orbit, it can transfer to the desired final conic 
in time x n , where, the sum of x n , x n , and x D is equal to x f . The fact that the angular 
momentum, the eccentricity, and the orientation of the intermediate closed orbit may be left 
free is what makes it possible to preassign arbitrarily large durations x fl and x n . This of 
course is not the only way to perform a transfer between two arbitrary conics in an 
arbitrarily large duration, and, with a little imagination, the reader can come up with some 
more ways for himself. 


4.13 Optimality of the Trajectories 

In this section, hints will be provided, suggesting that there should be cases in which 
the (tangential) thrust program proposed in this report (see Eqs. (4.10), (4.159)) results in 
sufficiendy opdmal trajectories, both for power-limited (PL) and constant ejection velocity 
(CEV) systems. 

The thrust program was obtained by using the tangential thrust assumption in the 
problem corresponding to minimum-fuel PL transfers. Numerical optimizations 44 45 47 49 51 
for such problems tend to suggest strongly, that there exists a large subclass of PL 
transfers between two coplanar elliptic orbits, corresponding mainly to cases where the 
changes in orientation and eccentricity are small, for which, as the duration of transfer 
increases, the thrust levels decrease, and the direction of the optimal thrust acceleration 
tends more or less to coincide with the direction of the tangent to the optimal flight path 
throughout the transfer. This leads one to expect that this thrust program (Eqs. (4.10), 
(4.159)) should turn out to be sufficiently optimal for at least some of these cases. An 
important subclass of such problems is the one corresponding to transfers between two 
circular orbits, and in this class, as was seen in Section 4.8, this thrust program is capable 
of resulting in solutions of very large duration. Also, as was seen in Sections 4.9 and 4.10, 
for escape problems from an arbitrary elliptical orbit, if one leaves either the angular 
momentum or the orientation co f of the final parabolic trajectory free, then one has the 
freedom to preassign the duration of the maneuver, making the thrust levels as small as one 


93 



wishes. To summarize, it appears that the best starting point for evaluating the performance 
of the thrust program (Eqs. (4.10), (4.159)) in the PL sense is to compare the exact optimal 
PL cost with the PL cost corresponding to the thrust program given in Eqs. (4.10), 
(4.159), first, for long duration escape problems from an arbitrary elliptical orbit, then, for 
long duration transfers between two coplanar circular orbits, and then, for transfers 
between two copjlanar, nonintersecting elliptical orbits having the same orientation. 

Based primarily on the results of Section 4.6.5, one could go as far as to conjecture 
here that, for power-limited, long duration (more than three or four revolutions around the 
planet) escape problems, and for long duration transfers between two coplanar circular 
orbits, the exact analytic solution of Section 4.4 must be very close to the exact optimal 
solution. 

A really interesting question is the optimality of the thrust program given in Eqs. 
(4.10), (4.159) with regard to other than power-limited propulsion systems. The costs 
corresponding to PL and CEV systems along any maneuver lasting for a time duration x f 

are given respectively by 4648 : 


T f T f 

JpL = J 2 ; Jcev = Jedx (4.163) 

0 0 

Table 4.1 supplies the PL and CEV costs for all the transfer trajectories given in the 
examples, as well as the corresponding initial and final values of the (nondimensional) 
thrust acceleration (Eq. (4.159)) £<, and e f for each trajectory. These costs were obtained by 
numerically evaluating the integrals in Eq. (4.163) along the (exact analytic) trajectories 
using a Simpson's Multiple-Segment 1/3 routine 54 . 

The PL cost in Table 4.1 is given for any future comparisons with the optimal costs 
corresponding to the exact optimal PL transfers. The CEV cost is given because it will now 
be compared with the cost corresponding to Hohmann (HM) and Biparabolic (BP) 
transfers 34 - 48 between the same initial and final conics for some of the examples. The 
Hohmann transfer is basic in the theory of impulsive transfers 34 - 48 and for coplanar orbits 
having the same orientation is performed (see Fig. 4.41) via two impulsive, tangential 
speed changes (AV's) at the perigee of the initial orbit and the apogee of the final orbit 34 - 48 . 
The perigee and apogee of an HM transfer orbit thus coincide with the perigee of the initial 
and the apogee of the final orbit, respectively. The BP transfer (see Fig. 4.42) is performed 
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Table 4.1 PL an d CEV costs along the transfers 

Transfer % ^ J PL 


Fig. 4.2 
Fig. 4.3 
Fig. 4.4 
Fig. 4.5 
Fig. 4.7 
Fig. 4.13 
Fig. 4.15 
Fig. 4.17 
Fig. 4.19 
Fig. 4.21 
Fig. 4.25 
Fig. 4.27 
Fig. 4.29 
Fig. 4.30 
Fig. 4.31 
Fig. 4.37 
Fig. 4.38 


5.9500e-3 

2.5500e-2 

5.9500e-2 

1.5610e-2 

6.1 139e-3 

3.6210e-2 

1 .0542e- 1 

1.0131e-l 

2.9601e-2 

5.9693e-3 

7.7549e-3 

7.2092e-3 

2.8987e-2 

5.0729e-3 

1.5905e-2 

1.5000e-l 

9.8358e-3 


2.9355e-3 

9.6038e-3 

7.007 le-3 

2.4088e-2 

1.5127e-3 

9.0432e-3 

5.5871e-2 

3.4644e-2 

2.6492e-3 

1.995 le-3 

1.1576e-3 

6.4353e-3 

3.4832e-2 

1.3587e-3 

7.4914e-4 

7.6719e-2 

5.7974e-3 


1.5369e-4 

7.3223e-4 

3.9537e-3 

3.7966e-3 

8.9225e-4 

5.4860e-3 

1.7064e-2 

2.0457e-2 

5.0593e-3 

1.4626e-4 

3.0947e-4 

2.6292e-3 

1.3490e-2 

5.7875e-4 

6.4252e-4 

6.7975e-2 

1.6055e-3 


•^CEV 

1.2861e-l 

1.6795e-l 

4.1077e-l 

3.0343e-l 

5.0001e-l 

5.0993e-l 

3.9424e-l 

5.9559e-l 

7.8669e-l 

1.6022e-l 

1.7577e-l 

4.5698e-l 

8.1437e-l 

3.7680e-l 

3.7963e-l 

1.1451e+0 

4.7716e-l 


by transferring to an (outgoing) parabola with an impulsive, tangential AV at the perigee of 
the initial orbit, then transferring to an (incoming) parabola at infinite distance from the 
planet (which takes zero cost) and then transferring to the final orbit with an impulsive, 
tangential AV applied at the perigee of the final orbit 48 . The BP transfer is only of 
theoretical importance because it takes infinite time to perform. Due to the impulsive 
character of the thrust, from Eq. (4.163), the cost for a HM or BP transfer is just the sum 
of the two impulsive speed changes required to perform the transfer. Table 4.2 compares 
the HM and BP costs for the transfers corresponding to Figs. 4.7, 4.13, 4.15, 4.21, and 
4.38 with the CEV cost corresponding to the present thrust program (second column), 
which is reproduced from Table 4.1. For Figs. 4.7, 4.13, 4.15, and 4.21 the HM cost is 
the minimum CEV cost possible 48 and one can see that the cost corresponding to the present 
thrust program is about 10 to 20 percent more than the HM cost for Figs. 4.7, 4.13, and 
4.15, 120 percent more for Fig. 4.21, but mostly less than the BP cost, except for Fig. 
4.21. For Fig. 4.38 the BP cost is the minimum CEV cost possible 48 and, instead of the 
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HM cost. Table 4.2 supplies the symmetric AD cost which corresponds to the sum of two 
impulsive, tangential AV's, the first accelerating, the second decelerating, applied at the 
apogees of the initial and final elliptic orbits respectively. The transfer trajectory in this case 
is a circular orbit. This symmetric AD transfer, contrary to the BP transfer, is quite a 
practical one 48 , and one can see from Table 4.2 that the cost corresponding to the present 
thrust program for Fig. 4.38 is about 10 percent more than the corresponding symmetric 
AD cost. 


Table 

4.2 Hohmann and 

Biparabolic 

CEV costs* 

Transfer 

Present J CEV 

Hohmann 

Biparabolic 

Fig. 4.7 

5.0001e-l 

4.4868e-l 

6.2132e-l 

Fig. 4.13 

5.0993e-l 

4.235 le-1 

5.6080e-l 

Fig. 4.15 

3.9424e-l 

3.5130e-l 

5.1015e-l 

Fig. 4.21 

1.6022e-l 

7.3568e-2 

9.8042e-2 

Fig. 4.38 

4.7716e-l 

4.3246e-l* 

9.8718e-2 


The above comparisons were supplied with the hope that they may convince some 
investigators that the near-optimality of the proposed thrust program in the CEV sense may 
deserve some further consideration. 


4.14 The Rendezvous Problem 

This last section constitutes a small deviation from the main subject (of transfers) and 
discusses the possibility of using the results of Section 4.6.2 for predicting the initial 
values of the costates corresponding to the problem of power-limited, minimum-fuel, 
coplanar rendezvous. For these problems the position of the vehicle on the initial and final 
conics is fixed. This is equivalent to specifying the state of the vehicle at the initial and final 
times completely, that is, the quantities Xq, r 0 , 0 O , h*, x f , r f , and 0 f , are fixed. Thus, for 
this problem, there are no transversality conditions, and the value of the (constant) costate 
P e is in general not zero, since the initial position of the vehicle can be expected to have an 
effect on the optimal value of the cost. The tangential thrust assumption was 
(mathematically at least) consistent with the transfer and escape problems, because, 
although the order of the original system of state-costate equations was reduced by one , 
two transversality conditions were satisfied by fixing the value of one constant of 
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integration (namely P e ) at zero. The tangential thrust assumption is, however, inconsistent 
with the rendezvous problem, since the eight boundary conditions are in general completely 
arbitrary, and one cannot expect to satisfy all of them simultaneously with fewer than eight 
constants of integration. However, one can assume, that if the duration of the maneuver is 
long, then the optimal thrust for the most part of the trajectory (away from the endpoints) is 
approximately tangential, and that for such long duration problems, the initial position of 
the vehicle probably has little effect on the optimal value of the cost. This means that one 
can combine Eq. (4.78), which was obtained using the tangential thrust assumption and 
zero P e , with Eq. (4.83), and predict the initial values of the costates: 


_ x 0 (h f -h 0 ) 

Vr 


p _( h f — ho) 
rh0 " *.2— 


( 4 . 164 ) 


P r0 = 


0 

1 

CN O 


{ <1 J 1 

l h 0 x f j 


P(> = 0 


( 4 . 165 ) 


These expressions for the costates may be used as reasonable guesses during a numerical 
optimization scheme aiming at the exact optimal solution of such problems. 


4.15 Summary 

This chapter has documented a case in which a mathematical model representing a 
meaningful physical system affords a closed-form solution. Recall that Keplerian two-body 
motion that takes place in the absence of thrust gives rise to a trajectory equation describing 
a conic (section). It was shown in the present chapter that the same trajectory equation (Eq. 
(4-17)) can be uncovered even in the presence of (a particular form of) nonzero, continuous 
thrusting terms. The difference is that in the former type of motion the angular momentum 
is a constant, while in the latter type of motion it is a linear function of time, with slope 
equal to a throttling parameter. Moreover, the corresponding thrusting terms are not really 
ad hoc, but arise from the optimization problem corresponding to power-limited, minimum- 
fuel, coplanar orbital motion. The assumption of tangential thrust plays a remarkable 
catalytic role in this problem, since it allows one to change the right-hand-side of only one 
costate equation (which affects only optimality), eliminate completely the costates by 
expressing them as functions of the states and a throttling parameter, and obtain the thrust 
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program for which the equations of motion can be solved exactly and analytically. It has 
been demonstrated in this chapter that it is, in principle, possible to solve arbitrary 
(coplanar) transfer and escape problems using this solution, for which one can preassign 
arbitrarily large durations. Transfers between two conics that intersect cannot be performed 
using a one-segment transfer trajectory. They can always be performed however using 
multiple-segment transfer trajectories, although the corresponding procedure is not unique. 
Hints were also given in the chapter suggesting that there should be cases in which the 
thrust program results in sufficiently optimal trajectories, both for power-limited and 
constant ejection velocity propulsion systems. 


4.16 Concluding Remarks 

Although the coverage in this chapter appears to be rather extensive, the work done is 
far from being complete. Above all, how the thrust program will perform in each and every 
particular case is far from having been documented, or even understood. Of primary 
interest for the future would be a more complete investigation of the transfer and escape 
problems studied above, and especially of the near-optimality of the thrust program in a 
general setting. It appears that, both Keplerian motion (zero thrust), and the kind of Non- 
Keplerian motion (nonzero thrust) documented in this chapter afford a common 
mathematical description. Physically, however, and philosophically, there is a significant 
difference between the two kinds of motion. Keplerian motion is guided by Nature , and 
requires no human effort. It represents a fundamental behavior of Nature. The kind of 
motion uncovered in this chapter on the other hand will have to be guided by humans, and 
there is nothing fundamental about it. This brings in a host of practical problems that need 
to be taken care of before a vehicle can actually be made to trace a path corresponding to 
such motion. With no misconceptions about this point, the main purpose of the present 
chapter on the practical side, was only to demonstrate, beyond a reasonable doubt, that the 
analytical results presented here may be of some help during actual planetocentric or 
heliocentric orbital operations. Thus, this report presents this thrust program and the 
corresponding exact analytic solution of the governing -equations of motion, with the hope 
that its complete features will be more extensively investigated, its near-optimal or non- 
optimal aspects as relate to different types of propulsion fully uncovered, and its full 
domain of practical applicability, if any, clearly identified. 
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Chapter V 


Conclusions and Recommendations 

The research in this report concentrated in three areas related to aircraft and 
spacecraft trajectory optimization and optimal guidance. Specifically, the attention was 
focused on the areas of unconstrained and constrained aircraft energy-state modeling and 
spacecraft motion under continuous thrust. 

5.1 Aircraft Unconstrained Energy-State Modeling 

The research in the area of aircraft unconstrained energy-state modeling has 
uncovered a systematic procedure for identifying the singular perturbation parameter in the 
differential equations of motion governing both conventional (subsonic-supersonic, flat 
Earth) and transatmospheric (hypersonic, spherical Earth) flight. A set of arbitrary scaling 
constants was used during the procedure to nondimensionalize all the variables of interest. 
Then, aided by a useful choice of the scaling constants, the conclusion was reached that 
two-time-scale behavior of the corresponding aircraft can be expected when the maximum 
longitudinal load factor during a maneuver remains sufficiently less than one. The 
important point was the validity of this statement regardless of the performance index being 
optimized. This explicit identification of the singular perturbation parameter also explained 
the past successes of singular perturbation treatments of aircraft energy climbs. 

A possible extension of this work would be a similar investigation of the existence of 
conditions under which the aforementioned differential equations exhibit three-time-scale 
behavior. It is straightforward to extend the procedure used in the second chapter by 
introducing two singular perturbation parameters instead of one, of which, one may be 
designated as £, and the other as e 2 . A useful choice for the scaling constants, combined 
with the values of these singular perturbation parameters, would then serve again as a 
means to uncover conditions under which three-time-scale behavior can be expected. For 
example, one could use such a procedure to investigate under what conditions, if any, are 
the energy, altitude, and flight-path angle dynamics of aircraft separated from each other, 
giving rise to three-time-scale behavior. Since it is well-known that usually the altitude and 
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flight-path angle dynamics are highly coupled, such a study would be of some importance 
if it succeeded in uncovering the underlying reason for this coupling in terms of explicit 
physical parameters, or in predicting the existence of possible cases for which the coupling 
is absent, giving rise to three-time-scale behavior. 


5.2 Aircraft State-Constrained Energy-State Modeling 

The research in the area of aircraft state-constrained energy-state modeling has so far 
shown that a transformation technique can be used to isolate and describe completely the 
class of asymptotic controllers that track a given state-constraint boundary. This result 
provided the incentive for proposing a reformulation for optimal control problems 
involving active state-variable inequality constraints, so that, in the reformulated problem 
the optimization is carried out only over the class of asymptotic controllers. The original 
problem leads to optimal controllers that are finite-time and one-sided. The reformulated 
problem leads to controllers that are approximately optimal, asymptotic, but still one-sided. 
If however the state constraint is regarded as a soft constraint, then one can find controllers 
that are asymptotic, two-sided, and result in the same optimal value of the performance 
index corresponding to the original problem, that is, they are practically optimal, but at the 
expense of violating the state constraint. From a singular perturbations point of view this 
suggests that such controllers can be used in a boundary-layer system, to track the reduced 
solution corresponding to a specific problem, when this reduced solution happens to ride a 
state-constraint boundary. However, such controllers do not correspond to stationary 
solutions of the optimization problem, so at the present, a systematic procedure for finding 
them does not exist. 

This last remark also suggests the primary recommendation for future research with 
regard to this area. Although a procedure was introduced in Section 3.3.2 that would 
presumably result in controllers that are asymptotic and two-sided, it is still not clear how 
to pick the two functions 'F and M (see Section 3.3.2) when formulating the related 
companion problem. If a meaningful way of picking these functions can be found, this 
formulation can be applied to supersonic or hypersonic energy climbs, for which part of the 
trajectory lies on a state-constraint boundary representing a constant value of dynamic 
pressure, aerodynamic heating rate etc. The resulting controllers will be nonlinear, 
practically optimal, asymptotic, and capable of tracking such boundaries from both sides. A 
first step toward this direction would be to try to find (by trial and error) an asymptotic. 
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two-sided controller, that for a simple model representing an aircraft energy climb will 
guide a vehicle along a constant dynamic pressure constraint boundary, and will result in 
the same optimal value of the performance index (minimum time for example) 
corresponding to the problem in which the dynamic pressure constraint is explicitly 
imposed. 


5.3 Spacecraft Motion under Continuous Thrust 

In the area of spacecraft motion under continuous thrust the primary result of this 
report has been a generalization of the exact analytic solution of the orbital equations of 
motion corresponding to Keplerian two-body motion (that takes place in the absence of 
thrust) to a type of motion performed under a special kind of continuous thrust. The most 
significant result of this report is most probably the fact that the trajectory equations 
corresponding to the two types of motion have identical form. The difference is that, in 
Keplerian motion the angular momentum is a constant, while in the type of motion 
uncovered in this report it is a linear function of time, with slope equal to a throttling 
parameter. A very important aspect of this work was the fact that the thrust program used 
for the exact analytic solution of the equations of motion was not really ad hoc, but arose 
from the optimization problem corresponding to power-limited, minimum-fuel, coplanar 
orbital motion. It has been demonstrated in the fourth chapter that it is in principle possible 
to solve arbitrary (coplanar) transfer and escape problems using this solution, for which 
one can preassign arbitrarily large durations. Transfers between two conics that intersect 
cannot be performed using a one-segment transfer trajectory. They can always be 
performed, however, using multiple-segment transfer trajectories, although the 
corresponding procedure is not unique. Hints were also given in the fourth chapter 
suggesting the existence of cases for which the thrust program results in sufficiently 
optimal trajectories, both for power-limited and constant ejection velocity propulsion 
systems. 

There are several suggestions for future research in this area. First, one may try to 
establish the near-optimality or non-optimality of the exact analytic solution uncovered in 
Section 4.4 in a full setting, by evaluating the exact optimal cost for a given maneuver and 
by comparing it to the cost corresponding to this exact analytic solution. 

For power-limited propulsion systems the above would be a rather formidable task, 
since it would require repeated numerical solution of two-point boundary value problems 
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for evaluating the power-limited optimal trajectories between many pairs of conics, and the 
corresponding (exact optimal) costs. At this point however, this task appears to be more 
significant in a mathematical, rather than in a practical sense, because, as argued upon in 
Section 4.12, the time duration for any type of (coplanar) escape or transfer problem, using 
the proposed thrust program, can be preassigned to be arbitrarily large. This implies (see 
Section 4.12) that the corresponding fuel consumption can (in principle) be made arbitrarily 
small if one is willing to compromise with regard to the time duration it takes to perform 
such a maneuver! Such a situation may appear as overly restrictive but it really isn’t! 
Balancing the fuel consumption for a maneuver against the time it takes to perform the 
maneuver is nothing new in spaceflight. The Hohmann transfer, which has been routinely 
used during the past thirty years in actual orbital operations, is a good example. Once the 
initial and final elliptical orbits are fixed, the time duration for a Hohmann transfer is fixed 
and cannot be preassigned! The Bielliptic transfer, which in practice can be used as an 
approximation to the Biparabolic transfer, represents just another such example, in which 
"reducing the fuel consumption further” means that one "has to wait longer" for completing 

the transfer. 

For constant ejection velocity propulsion systems the situation may be somewhat 
easier to deal with, since there are many cases of transfers for which the optimal constant 
ejection velocity cost, impulsive in character, can be calculated with relative ease. 

Another recommendation for future research would be a continuation of the work 
started here on the existence or non-existence of one-segment solutions. Theorem 4.1 of 
Section 4.1 1 started at first as a conjecture, and then was rigorously proven in Appendix F. 
Similarly, one now has the conjecture on the existence of one or two-segment solutions 
given at the end of Section 4.1 1 that awaits to be shown true or false. 

Still other possible recommendations for future research would include an extension 
of the work to the case of three-dimensional motion, and a regular perturbations point of 
view of the solution uncovered in Sections 4.3 and 4.4. Although three-dimensional 
transfers cannot be performed by using exclusively tangential thrust, the set of state 
variables introduced in Section 4.2 can be generalized to the case of three-dimensional 
motion, after which, the three-dimensional power-limited optimization problem can be cast, 
and the structure of the corresponding costate equations examined. Viewing the solution 
uncovered in Sections 4.3 and 4.4 as the lowest order solution of the state-costate 
equations corresponding to a regular perturbation expansion would first require the 
nontrivial task involving a clear identification of a regular perturbation parameter. 
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Appendix A 


Main Features of Two-Body Keplerian Motion 

This Appendix is supplied for a direct comparison between the main features of 
Keplerian motion, that takes place in the absence of thrust, and the motion uncovered in 
Section 4.4, that takes place in the presence of continuous thrust, specified by the thrust 
acceleration components proposed in Eq. (4.10). 

The term "two-body Keplerian motion" refers to the translational motion that a space 
vehicle (considered as a point mass) executes in the vicinity of a spherical, homogeneous 
planet, in the absence of any thrusting forces. The only force acting on the vehicle during 
such motion is the inverse-square gravity from the planet. This translational motion is 

confined to a plane, and is fully described by Eqs. (4.6)-(4.9), with the thrust components 
£ r and Eg set equal to zero: 



(A. 1) 


dx _ h 2 - r 
dx r 3 


(A.2) 


dr 

dx 


= x 


(A.3) 


de = _h 

dx r 2 


(A.4) 


The solution of the differential equations describing such motion was first given, of course, 
by Isaac Newton 37 , and in terms of the variables used in Eqs. (A.l) through (A.4) can be 
summarized as: 


h = const. 


esin(0 - to) 
~h 


r = 


h 2 

l + ecos(0-co) 


(A.5) 
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In Eq. (A.5) e and 0) are two integration constants, and the angular momentum h, 
which stays constant, plays the role of a third integration constant. The fourth integration 
constant relates the time x to the argument of latitude 0 and can be uncovered by using the 
expression for r from Eq. (A.5) and reducing Eq. (A.4) to a quadrature: 


=J 


h 3 d0 


[l + ecos(0-u>)]' 


• + const. 


(A. 6) 


According to the trajectory equation for r given in Eq. (A.5), the motion takes place 
along a conic 34 . The constant e is just the eccentricity of the conic. For e = OorO<e< 1 
the conic is respectively a circle or an ellipse, and the corresponding trajectory is closed and 
usually called an orbit. For e = 1 or e > 1 the conic is respectively a parabola or a 
hyperbola, and the corresponding trajectory is open 34 . The evaluation of the integral in Eq. 
(A.6) depends strongly on whether e is less than or greater than one, or equal to zero or 
one. The constant CD represents the orientation of the conic. If one imagines a vector 
pointing from the center of the planet toward the point on the conic where r is a minimum 
(the perigee), then cd is just the angle, measured anticlockwise, from the fixed direction in 
space (where all the latitude angles are measured from - see Fig. 4.1) to that vector 34 . 

An important parameter associated with a conic is its semimajor axis a. The 
semimajor axis is a direct measure of the mechanical (kinetic plus potential) energy per unit 
mass of the vehicle 34 . Because the inverse-square gravity field generated by the planet is 
conservative the mechanical energy per unit mass of the vehicle stays constant during the 
motion 34 . Specifically, this mechanical energy is related to the semimajor axis a by the well- 
known energy integral 34 : 


V 2 -- = -- (A.7) 

r a 

In Eq. (A.7) V is the speed of the vehicle, which in terms of the state variables used in Eqs. 
(A.l) through (A.4) is given by: 


V 2 



(A. 8) 


The semimajor axis a is related to the angular momentum h and the eccentricity e through 34 : 
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h 2 = a(l -e 2 ) 


(A.9) 


When the conic is a circle or an ellipse e is less than one and a is positive. Specifically, 
when e is zero the conic is a circle and a is just the radius of the circle. When the conic is a 
hyperbola e is greater than one and a is negative. When the conic is a parabola e is equal to 
one and a is infinite 34 (but the angular momentum h is always finite). 
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Appendix B 


Derivation of the Exact Analytic Solution 
of the System of Eqs. (4.11) - (4.14) 


The development in this Appendix is valid for nonzero A (the case with zero A 
corresponds to Keplerian motion and was summarized in Appendix A). The solution of Eq. 
(4. 1 1 ) is obviously the linear variation for the angular momentum h supplied in Eq. (4. 1 5). 
Consider now the system of Eqs. (4. 12), (4. 13), and define the transformation: 


w _ h dw_ l^dh^_ 1 

W_ A ’ dx ~ AAdxJ" A 


(B.l) 


Using w rather than X as the independent variable, this system assumes the form: 


dx _ A 2 w 2 - r } x 
dw r 3 w 


(B.2) 


dr 

dw 


= x 


(B.3) 


By switching from the variable x to a new variable P, defined by: 


P = 


w 


(B.4) 


the system of Eqs. (B.2), (B.3) can be written as: 

dP _ A 2 w 2 — £ 
dw wr 3 


(B.5) 


_dL 

dw 


= wP 


(B.6) 


One can now define two new variables K and z, by: 
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(B.7) 



z = 


r 



Replacing r by z, and using z as the independent variable, the system of Eqs. (B.5), (B.6) 
transforms into: 


dP_ (z-A 2 )K 2 
dz (2z-P)z 3 


(B.8) 


dK _ 2K 
dz 2z - P 


Finally, using the transformation: 


S = 


2K 

2z-P 


the system of Eqs. (B.8), (B.9) assumes the analytically soluble form: 


dS 

dz 




z- A 
4z 3 



(B.9) 


(B.10) 


(B.ll) 


dK 

dz 


= S 


(B.12) 


The general solution of Eq. (B.l 1) can be found by a simple integration 33 . The result is: 


S 2 = 


4z 2 

Dz 2 +2z- A 2 


(B.13) 


where D is an integration constant. Using this result, Eq. (B.12) can be written as: 


dK V _ 4z 2 

dzj Dz 2 +2z-A 2 


(B.14) 


More will be said about Eq. (B.14) at the end of this Appendix. First, Eq. (4.14) in 
the main text can be rewritten as a simple quadrature: 
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6 =J(7> +f 


(B.15) 


where F is an integration constant. Tracing the transformations back, Eq. (B.15) can 
further be written as: 



A 




zVDz 2 + 2z- A 2 J 


dz + F 


(B.16) 


Note that the double sign in front of the integral is due to the fact that both signs have been 
kept when taking the square root of S in Eq. (B.13). Evaluation of the integral 33 in Eq. 
(B.16) for any value of D, and for nonzero A (zero A corresponds to Keplerian motion) 
leads, for both signs , to the expression: 


A 2 

1 + Bcos(0-C) 


(B.17) 


where B is a nonnegative constant defined by: 

B = Vl + DA 2 ( B18 ) 


and C is a constant which can be written always as a linear function of the constant F. 
Since z is just equal to rA 2 /h 2 , Eq. (B. 17) can be written as: 


h 2 

l + Bcos(8 — C) 


(B.19) 


which just happens to be the transfer trajectory equation, Eq. (4.17). Differentiating Eq. 
(B.19) one can easily find the corresponding expression for x given in Eq. (4.16). 

The last equation that remains to be integrated is Eq. (B.14). However, a simple 
inspection reveals that this integration has in fact been already carried out (in disguise) in 
Section 4.5 of the main text, by defining more useful variables, such as the generalized 
eccentric and hyperbolic anomalies, etc. There is therefore no need to integrate Eq. (B.14). 
This completes the exact analytic solution of Eqs. (4.1 1) through (4.14). 
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Appendix C 


Power-Limited Propulsion Systems 

Power-Limited propulsion systems basically correspond to electric propulsion, and 
generate thrust by accelerating particles through an electromagnetic or electrostatic 
field 46 ' 4851 . There is usually a propellant feeding mechanism that supplies the particles, and 
an electrical generator that generates the electric field that accelerates the particles. The 
primary characteristic of such systems is that the above two mechanisms operate (ideally) 
independently of each other. For example 51 , if a nuclear reactor heats a working fluid, as in 
an electrothermal device, changing the operating temperature of the reactor regulates the 
mass flow rate of the particles -dm/dt. If a solar cell accelerates ions, as in an electrostatic 
device, changing the operating voltage of the cell controls the ejection speed of the ions c. 
Thus, it is possible in principle to modify both the ejection speed c and the mass flow rate 
-dm/dt, and therefore to control independently the thrust T and the power P, given by the 
ixpressions: 





(C.l) 


The above expressions can be obtained by considering the state of the vehicle and the 
propellant at times t and t+dt, and by performing a momentum and energy balance. That the 
power in Eq. (C.l) is indeed the power supplied by a power-limited thruster can be further 
argued upon as follows. The power P supplied by such a thruster can be written as 46 48 : 

P = UI (C.2) 

A'here U is the beam voltage and I is the beam current. If a particle of mass m, charge q, 
md negligible initial velocity is accelerated through a potential difference U it acquires 
dnetic energy equal to (mc 2 /2), where c is the ejection speed. Accordingly, the following 
wo relations are valid: 
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Substituting from Eq. (C.3) into Eq. (C.2), one obtains the expression for the power given 
in Eq. (C.l). As the name implies, the power that can be supplied by such propulsion 
systems is limited, that is, there is a maximum level of power P max that cannot be 
exceeded 48 51 . In the T vs (-dm/dt) operating domain constant power is represented by the 
parabola: 



(C.4) 


The two relations in Eq. (C. 1) lead to: 



*0 


(C.5) 


where a=T/m is the thrust acceleration, given by: 


a 2 = Er + E 


2 

e 


(C.6) 


The initial mass is given, and in order to minimize the fuel consumption one must 
maximize the final mass m,. Therefore, it is required to operate with maximum power 
P=P max during the entire maneuver, which leads to the minimization of the quadratic 
performance index: 



(C.l) 


In practice 48 , there is usually a minimum ejection speed c mri and a maximum mass 
flow rate (-dm/dt)^. These constraints however need not be taken into account explicitly 
during a formulation of the corresponding optimal control problem, since by selecting the 
time duration for a maneuver appropriately one can ensure that these limits are not violated. 
The primary reason that one can usually do this is the fact that the optimal thrust 
acceleration levels and the optimal cost are inversely proportional to the time duration 
of a maneuver. 
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Nondimensionalizing a and K as (see Eqs. (4.4) and (4.5)): 


J f = 


K 


g s v s 



£r+£e 


(C.8) 


where J f represents the final value of the variable J (the initial value of J can be taken as 
zero), leads to Eq. (4.43), describing the time evolution of the cost J(t) along a maneuver. 
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Appendix D 


Coefficients of Polynomial Eq. (4.135), and an Explicit 
Solution^ for Transfers Between Coplanar Circular Orbits 


The coefficients of the fourth order polynomial equation (4.135) 
following expressions: 

are given by the 

b 4 =e 2 +4hoA 2 (l-e 2 )-4hfA 2 -8h 2 h 2 A 2 e f 

(D.l) 

b 3 = b, = 8hgAe f 

• (D.2) 

b 2 = 2[e 2 + 4h* A 2 (l + e 2 ) - 4h? A 2 ] 

(D.3) 

b 0 =e 2 +4hoA 2 (l-e 2 )-4hfA 2 +8hoh 2 A 2 c f 

(D.4) 


D.l Long Duration Transfers Between two Coplanar Circular Orbits 

Consider the case in which the final orbit is circular, that is, e f = 0. In such a case the 
polynomial Eq. (4.135) has no real roots, and can be satisfied only if A is zero. But when 
A is zero the motion is Keplerian and (assuming that the initial and final circular orbits are 
different) there can be no transfer. Thus, mathematically, the case with e f = 0 is a singular 
case. In practice however, the transfer between two circular orbits can be performed by 
selecting a very small e f . This immediately implies that the time duration for such a transfer 
will be large, that is, for the problem of transferring between two coplanar circular orbits 
long duration is implied if one uses the thrust program proposed in Eq. (4.10) and the 
corresponding exact analytic solution of Section 4.4. Small e f implies large x { and small A. 
Thus, keeping only the lowest order terms in A and e f , the coefficients given in Eqs. (D.l) 
through (D.4) assume the form: 
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b 3 = b, = k 2 = 8h 2 Ae f (D.5) 


b.=^ = b.=k,=e f 2 + 4(h;-h;)A 2 ; 

and now, the quartic polynomial Eq. (4. 135) can be simplified as: 

kj y 4+ k 2 y 3 + 2k,y 2 + k 2 y + k, =0 (D.6) 

Equation (D.6) can be easily factored as: 

(y 2 + l)(kj y 2 + k 2 y + k,) = 0 (D.7) 

from which one concludes that in such a case there can be at most two real solutions for y, 
and this happens only when the discriminant: 

A, = k 2 -4kf (D.8) 


is nonnegative. Substituting from Eq. (D.5) into Eq. (D.8), this condition can be stated 
explicitly as: 


A, =-4[l6(h^-h f 4 ) 2 A 4 -8(h 4 + h 4 )A 2 e f 2 +e f 4 ]>0 (D.9) 


Whenever the above condition is satisfied there are two real roots for y, corresponding to 
the initial value 0 O of 0, which determine the starting point on the initial circular orbit. 
These two real roots are given by: 


y .. 2 = 


-k 2 ±^k\ -4k)* 
2k, 


(D.10) 


For a given (small) e f , real solutions for y exist only for x f) < T f < x n , where x fI and x n are 
the two values of the final time x f for which A, becomes zero. Figures 4.6 and 4.8 through 
4.11 depict this situation very clearly. Using the expression supplied in Eq. (D.9), and 
assuming that the angular momenta fig and h, of the initial and final circular orbits are 
always strictly positive, the two values x n and x n can be found as: 

(i) For hj > ho and A > 0: 

, . _2(h; + h2)(h,-h„) 

c fi » x f2 (D.ll.a) 

e f e f 
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(ii) For hf < ^ and A < 0: 


2(h;-h;Xh,-h t ) 


2(h;+h;Xh»-h r ) 


(D.ll.b) 


Equation (D.l l.a) is in complete agreement with the situation depicted in Figs. 4.6, 
and 4.8 through 4.1 1. With 1^ = 1 and h, = 2, one obtains from Eq. (D.l l.a) x fl = 6/e f and 
x n = 10/e f . With e f = 0.05 these relations yield x f , = 120 and x n = 200, in complete 
agreement with Fig. 4.6. Similar results can be validated for Figs. 4.8 through 4.1 1. 

If now one preassigns a (large) value for the transfer duration x f , in order to make the 
final eccentricity as small as possible, one should obviously choose the (double) root for y 
corresponding to x n . Thus, in such a case, the final eccentricity e f and the corresponding 
value of the throttling parameter A are given by: 


2(h;-h;Xh,-h,) 



(D.12) 


Using this expression for A, and keeping in mind that A, is zero, the (double) root for y 
from Eq. (D. 10) can be found as: 


y = - 



(D.13) 


Since y is equal to tan(0(/2) (see Eq. (4.136)), the initial value of the argument of latitude is 
simply: 


e 0 = - (D.14) 

0 2 

The final value of the argument of latitude can be found substituting the above results for A 
and 0 O into Eqs. (4.137) and (4.138). The result is: 

0 f =0 fP + 2k7t ; e fr = f (D-15) 

Now, what about (the integer) k? k fixes the number of revolutions about the planet during 
the transfer. The value of k is the one for which Eq. (4.125.b) is satisfied. For (long 
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duration) transfers between two coplanar circular orbits Eq. (4.125) decouples from the 

rest of the problem. Using Eqs. (4.129) and (D.12), the generalized eccentricity of the 
transfer trajectory is: 


B = ± 



(D.16) 


Using now Eqs. (D.14) through (D.16), and (4.126) one concludes from Eq. (4.125. b) 
that k is: 


k = 


T f (h f + h 0 ) 

47thp hf 


T3/2 


( h f _h o ) 2 


(D.17) 


Any nonnegative integer value of k is acceptable. Thus, for large t f , if x f is arbitrarily 
preassigned, then, Eq. (D. 1 7) is approximately satisfied. In fact, by playing around with 
the values of x f and e f , connected through Eq. (D.12), one can satisfy Eq. (D.17) exactly, 
that is, one can make Eq. (D.17) result in an integer Value for k. For very large t f one can 
also simplify Eq. (D.17) as: 


__ x f (h f + h 0 ) 
47t h^ hf 


(D. 18) 


D.2 Summary 

To summarize, consider the problem of transferring from an initial circular orbit 
(e 0 =0) with angular momentum 1^ to a final coplanar circular orbit (e^O) with angular 
momentum h f , using the thrust program proposed in Eq. (4.10) and the corresponding 
exact analytic solution supplied in Section 4.4. The solution can be summarized as follows: 

(a) One preassigns a large transfer duration t f . 

(b) Eq. (D.12) then fixes the final (small) eccentricity e f and the throttling parameter 
A to be used for the transfer. 

(c) The value of the right-hand-side of Eq. (D. 1 7), truncated to the nearest integer, 
supplies the number of revolutions k about the planet during the transfer. By selecting t f 

appropriately one can satisfy Eq. (D.17) exactly by making its right-hand-side come out 
equal to an integer. 
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(d) The transfer starts and ends at the same argument of latitude after k full 
revolutions about the planet. 

(e) The generalized eccentricity and orientation constants B and C of the transfer 
trajectory are given by Eq. (4.129). 

(f) The transfer trajectory is given by Eqs. (4.17) or (4.29), where the angular 
momentum depends on the generalized eccentric anomaly E through Eq. (4.34). 

(g) If x f is large enough, further simplifications are obtained by considering B as 
being approximately equal to zero. In this case the transfer trajectory is simply r = h 2 , while 
the angular momentum h along the transfer is given by Eq. (4.25). Thus, for large enough 
x { one obtains the following spiral, representing the one-segment transfer trajectory: 


r = h = 


i-2Ah*(e-e 0 ) 


(D.19) 


(h) Also, when x f is large enough and B can be approximated by zero, the power- 
limited cost along the transfer is approximately given by Eq. (G.19.c) of Appendix G: 


J 


PL “ 


' hf-h 0 

6t f 


V J_ 
A^o 



(D.20) 


This expression for the power-limited cost compares well with the exact optimal 
power-limited cost supplied in Table 5.7 of Ref. 51. Based primarily on the results of this 
Appendix, and of Section 4.6.5, one could conjecture that for long duration transfers 
between two coplanar circular orbits the solution delineated in steps (a) through (h) above 
is the optimal power-limited solution, and the cost given in Eq. (D.20) is the optimal 
power-limited cost for all practical purposes. 
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Appendix E 


Coefficients of Polynomial Eq. (4.150) 

The coefficients of the fourth order polynomial equation (4.150) are given by the 
following expressions: 

b 4 = (l - e 0 ) 2 (e 2 + e f 2 - 2e 0 e f cos co f ) 

“0 ~ e o) 2 ( 2h f. A + e 0 e f sin© f ) 2 +4hgA 2 (l - e 2 cos 2 to f ) 

+2h 2 Ae f (l - e 0 )(2 sin co f - 4h 2 A cos© f - e 0 e f sin 2co f ) (E. 1 ) 

y = 4h 2 A(l - e 0 )(e f cosa> f - e 0 ) + 4hoA 2 e 2 sin2a> f 

+4h 2 Ae f (l - e 0 )(2h 2 A + e 0 e f sinco f )sino) f (E.2) 

2 = ( ! “ e o)( e o + e f “ 2e o e f coscOf)- (l -e 2 )(2h 2 A + e 0 e f sin(o f f 
+2hoAe 0 e f (2sinco f -4h 2 Acosco f -e 0 e f sin2to f ) 

44hgA 2 (l + e 2 cos 2 to f -2e 2 sin 2 © f ) (E.3) 
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= 4hoA(l + e 0 )(e f cosco f -e 0 )-4hoA 2 e 2 sin2co f 
+4hoAe f (l + e 0 )(2h 2 A + e 0 e f sinco f )sinco f (E.4) 

b 0 = (l 4- e 0 ) (cq e f — cos(D f j 

-(l + e 0 ) 2 (2h 2 A + e 0 e f sinto f ) 2 +4hoA 2 (l -e 2 cos 2 a> f ) 

-2hoAe f (l + e 0 )(2sin© f -4h 2 Acosa) f -e 0 e f sin2to f ) (E.5) 

Note that for e 0 = co f = 0 the expressions in Eqs. (E.l) through (E.5) reduce to the 
corresponding expressions given in Eqs. (D.l) through (D.4), valid for transfer and escape 
problems from a circular orbit. 
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Appendix F 


Proof of Theorem 4.1, Section 4.11 

This Appendix supplies the detailed proof of Theorem 4.1 (page 89) of Section 
4. 1 1 , which is reproduced here for convenience to the reader: 

Theorem 4.1. Consider the problem of transferring from conic (ho.eo.cOo.Go) to 
conic (h^epCO^Gf) using the thrust program proposed in Eq. (4.10) and the corresponding 
exact analytic solution of Section 4.4. Then, 

(i) If (h o ,e o ,(O o ,0 Q ) and (hpepCOpGp) have more than one common point (that is, if 
they intersect), then a one-segment transfer solution does not exist. 

(ii) If (h 0 ,e 0 ,co 0 ,G 0 ) and (hf,e f ,(o f ,0 f ) have only one common point (a point of 
tangency) then there exists at most a single impulsive one-segment transfer solution, 
performed at the point of tangency of the two conics, for which A is infinitely large, x f is 
infinitely small, and Ax f is finite and equal to hf-h^ 

Proof 4.1. Since only the relative (and not the absolute) orientation of the two conics 
on the plane affects the problem, there is no loss of generality in taking co 0 as zero. First, 
one can show part (i) using contradiction. Assume that, contrary to the claim made in the 
theorem, the initial and final conics intersect and a one-segment transfer trajectory 
(A,B,C,x f ) joining the conics with x f > 0 does exist. Then, by equating the right-hand-sides 
of the trajectory equations: 

r = h o . r _ hf /p 

l + e o cos0 ’ 1 + e f cos(0-to f ) 

corresponding to the two conics, one obtains the equation: 

hf - hg + (h f 2 e 0 - hge f cosco f )cos0 - h£e f sinco f sin0 = 0 (F.2) 

Using now the definition for y, and the identities 33 : 

y = tan (T) : sin 0 = — ; cos0 = -L^- (F.3) 

W 1 + y 1 + y 2 
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Eq. (F.l) can be manipulated into the form: 

a 2 y 2 + 2a,y + a 0 = 0 (F.4) 

where the quantities a^ a,, and af, are given by: 

a 2 = hf(l-e 0 )-ho(l-e f cosco f ) ; a, =-h 2 e f sino) f ; a 0 = hf(l + e 0 )-h 2 (l + e f costo f ) 

(F.5) 

Equation (F.4) is a quadratic polynomial equation in the unknown y (corresponding 
to 0), and since by assumption the initial and final conics intersect, it should have at least 
one real solution. Accordingly, the discriminant A of Eq. (F.4), given by: 

A = af - a 2 a 0 (F.6) 

should be nonnegative. Substituting from Eq. (F.5), this necessary condition can be 
written explicidy as: 


A = h^(e 2 -l) + h^(e 2 -l) + 2h 2 h 2 (l-e 0 e f costo f ) > 0 (F.7) 

Keeping in mind that the initial and final conics have a common focus, one can now 
conclude that: 

(i) If the initial and final conics have only one point in common (a point of 
tangency), A should be zero, resulting in one real solution for y (of multiplicity two). 

(ii) If the initial and final conics have two points in common (that is, if they 
intersect), A should be strictly positive, resulting in two distinct real solutions for y. 

(iii) Two (different) conics (with a common focus) can have no more than two points 
in common. 

Since by assumption a one-segment transfer trajectory (A,B,C,x f ), with x f > 0, 
joining the initial and final conics exists, the quantities appearing in Eq. (F.7) are all 
dependent on the characteristics A, B, C, and x f > 0 of this one-segment transfer trajectory, 
through the boundary conditions given in the main text by Eqs. (4.139) through (4.142) 
and also through one among Eqs. (4.25), (4.26), (4.34), or (4.42), written at the final time 
x f . To prove the theorem one must uncover this dependence explicitly and show that A is 
strictly negative, contradicting the necessary condition given in Eq. (F.7). With this in 
mind, one can first use the definitions: 
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^o=0o-C 


— 0 f C 


(F.8) 


and write Eqs. (4.141), (4.142) as: 

e f cos(£ f + C - co f ) = Bcos^ f (F.9) 

e f sin(^ f + C - to f ) = 2r f A + Bsin (F.10) 

Now, Eqs. (F.9) and (F.10) can be solved for e,sinco f , e,cosa) f , and the square of e f to 
yield: 

e f sin ( 0 f = BsinC-2r f Acos(£ f +C) (F. 11) 

e f cos(O f = BcosC + 2r f Asin(£ f +C) (F.12) 

e 2 = B 2 + 4A 2 r 2 + 4ABr f sin£ f (F. 13) 

Using Eqs. (F.12), (F. 13), and the definitions: 

Po = h o ; Pf =hf (F.14) 

one can write the discriminant A given in Eq. (F.7) as: 

A = Pf ( e o ~ 0 + Po(B 2 - 0 + 2 PoPf + 4 ABp 2 r f sin £ f - 2B(e 0 cos C)p 0 p f 

+ 4 A 2 p 2 r 2 - 4 A(e 0 cos C)p 0 p f r f sin - 4 A(e 0 sin C)p 0 p f r f cos % f (F. 15) 

In an analogous manner, one can write Eqs. (4.139), (4.140) as: 

e 0 c°s(£o + C) = Bcos£ 0 (F.16) 

e oSin(£o + C) = 2r 0 A + Bsin£ 0 (F.17) 
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and solve Eqs. (F. 16) and (F. 17) for e 0 sinC, e 0 cosC, and the square of e 0 to obtain: 


e 0 sin C = 2r 0 A cos 

(F.18) 

e 0 cosC = B + 2r 0 Asin 

(F.19) 

eg = B 2 + 4 A 2 t 2 0 + 4 ABr 0 sin 

(F.20) 


Using Eqs. (F.18) through (F.20), one can further write the discriminant A given in Eq. 
(F.15) as: 

A = (p f - Po) 2 (B 2 - 1) + 4 A 2 (p 2 r 2 + p 2 r 2 - 2p 0 p f r 0 r f cos(£ f - £ 0 )) 

+ 4 AB(p f - p 0 )(p f r 0 sin £ 0 - p 0 r f sin £ f ) (F.2 1 ) 

Note that the radial distances r 0 and r f at the initial and final times can be expressed more 
explicitly using the trajectory equation (Eq. (4.21)) and the definitions given in Eq. (F. 14) 
as: 

r = Po • r = El CF 22) 

0 1 + Bcos^ 0 ’ f l + Bcos^ f 

To proceed further with the proof one must now express more explicitly the six 
quantities r 0 , r f , r 0 sinq 0 , r 0 cos£ 0 , r,sin£ f , and r,cos^ f . This step however depends strongly 
on whether the generalized eccentricity B of the transfer trajectory is zero, one, between 
zero and one, or greater than one. One therefore has to consider the following four possible 
cases: 

F.l First case: B = 0 

In this case, from Eqs. (4.139) and (4.140) one can deduce the following two 
possibilities for the starting point 0 O on the initial conic: 

if A > 0 then 0 O = — ; if A < 0 then 0 O = -^- (F.23) 

0 2 0 2 
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From Eq. (F.22), with B = 0 one obtains: 


r o = Po ; r f =Pf (F.24) 

Also, one can define the change A0 or A t, in the argument of latitude during the one- 
segment transfer as: 


A6 = e, -8„ =8, -C + C-e 0 = %, -5, = A| (F.25) 

Note that, since the transfer duration is strictly positive, that is, T f > 0, one must necessarily 
have A0 > 0 for the transfer. Along such a one-segment transfer, with B = 0, Eq. (4.25) is 

valid. Using the definitions given in Eqs. (F.14) and (F.25), one can write Eq. (4.25) at 
the final time x f as: 


— = — — 2AA0 (F.26) 

Pr Po 

By defining the quantity D as: 

D = 1 - 2Ap o A0 (F.27) 

p f can be written from Eq. (F.26) as: 

Pf=*j (F.28) 

The last step consists of substituting from Eqs. (F.24) and (F.28) into the expression for A 

given in Eq. (F.21). With B = 0, and for both possibilities given in Eq. (F.23), one can 
write the result as: 


A = — 


16 A 2 hp 

D 2 



(F.29) 


For A0 > 0 the following inequality is always valid: 


A0 . (AS 
— > sin — 
2 2 


(F.30) 
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Accordingly, from Eq. (F.29), the discriminant A must be strictly negative (A < 0), 
contradicting the necessary condition A > 0 found in Eq. (F.7). This concludes the proof of 
part (i) for B = 0. 


F.2 Second case: B = 1 


In this case, to facilitate the algebra, one can again define the two quantities: 

'SA • 


y 0 =lan 


'£ } 
Yf = tan| f 
2 J 


(F.31) 


Using Eq. (F.31), and some standard trigonometric identities 33 , one can express the six 
quantities r 0 , r/, r 0 sin^, r 0 cosl^, r^in^ f , and r,cosl^ as (recall that B = 1): 


Po(l + yo) 
2 

Pr(i+y?) 


r 0 sin^ 0 = p 0 yo ; r 0 cos^ 0 = P ° ' ~ — ( F - 32 ) 

r,sin|,=p,y, ; r,cosl;,=?£^ (F.33) 


Along such a one-segment transfer, with B = 1, Eq. (4.26) is valid. Using the definitions 
given in Eqs. (F.14) and (F.31), one can write Eq. (4.26) at the final time T f as: 

— “ — - A(y r - y,,) ~(yj - yj) CF.34) 

D. D» ' 


By defining the quantity D as: 


D = 3 - Ap 0 (y f - y 0 )(3 + y) + y f y 0 + yo) (P-35) 


p f can be written from Eq. (F.34) as: 


Pf = 


3 Po 

D 


(F.36) 


The last step consists of substituting from Eqs. (F.32), (F.33) and (F.36) into the 
expression for A given in Eq. (F.21). With B = 1, one can write the result as: 
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(F.37) 


A = — 


3A 2 h^ 

D 2 




For B 1, a full revolution around the planet during a one-segment transfer is not possible, 
because 0, or £ go through a value at which r becomes infinite, as can be seen from Eq. 

(4.21). Thus, along such a transfer 0 < § f - < 2k, and the following statement is always 

true: 


tan 




(F.38) 


Accordingly, from Eq. (F.37), the discriminant A must be strictly negative (A < 0), 
contradicting the necessary condition A > 0 found in Eq. (F.7). This concludes the proof of 
part (i) for B = 1. 


F.3 Third case: 0 < B < 1 

In this case, the generalized eccentric anomaly E, introduced in Section 4.5.2, rather 
than 0 or £, is the "natural” angular coordinate to use in the expression for A. By defining 
the quantities: 


b = (l-B 2 ) 3/2 ; 

AE = E f -E 0 

(F.39) 

and using the results of Section 4.5.2, one can 
r 0 cos£o, r fSin£ f , and r,cos£ f as: 

express the six quantities r 0 , r (1 

r o si n^ 0 ’ 

p 0 (l-BcosE 0 ) P 0 sin E n 

o b 2/3 > r o sin So b , /3 

• r _Po(cosE 0 -B) 

’ r 0 COS S 0 b 2/3 

(F.40) 

p f (l-BcosE f ) p f sinE f 

f b 2/3 ’ b l/3 

; r,cos^ r - p '( cosE >- B ) 

f ^2/3 

(F.41) 


Along such a one-segment transfer, with 0 < B < 1, Eq. (4.34) is valid. Using the 
definitions given in Eqs. (F.14) and (F.39), one can write Eq. (4.34) at the final time x f as: 
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(F.42) 


1 i 2 A[ AE - B(sin E f — sin E 0 )] 

Pr Po b 

By defining the quantity D as: 

D = b-2Ap 0 [AE-B(sinE f -sinE 0 )] (F.43) 


p f can be written from Eq. (F.42) as: 



(F.44) 


The last step consists of substituting from Eqs. (F.40), (F.41) and (F.44) into the 
expression for A given in Eq. (F.21). One can write the result as: 


A = 


16A 2 h*(l-B 2 ) 

D 2 


¥) 



(F.45) 


Since the transfer duration is strictly positive, that is, x f > 0, it must necessarily be true that 
A0 = 0 f - 0 O > 0, which in turn implies that A£ = > 0 and AE > 0. Moreover, for 

AE > 0 the following inequality is always valid: 


AE 

2 



(F.46) 


Accordingly, from Eq. (F.45), the discriminant A must be strictly negative (A < 0), 
contradicting the necessary condition A > 0 found in Eq. (F.7). This concludes the proof of 
part (i) for 0 < B < 1. 


F.4 Fourth case: B > 1 

In this case, the generalized hyperbolic anomaly H, introduced in Section 4.5.3, 
rather than 0 or is the "natural" coordinate to use in the expression for A. By defining 

the quantities: 

b = (B 2 - 1) 3 2 ; AH = H f - H 0 (F.47) 
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and using the results of Section 4.5.3, one can express the six quantities r 0 , r f , r 0 sin^ 0 , 
^cos^, r^in^, and r,cosi^ as: 


Mgcoshty-l) . rtinf _ ppSinhH,, , p 0 (B-coshH„) 

o , 2/3 ’ r o 5,111 So ,i/3 > r o cos So 7173 (r.4») 


r = P r (BcoshH, -1) r _ p,sinhH, . _ p f (B-c 0 shH f ) 

f u 2/3 ’ f sin Sf ,i/3 > r f cosq f - — ^ 


(F.49) 


Along such a one-segment transfer, with B > 1, Eq. (4.42) is valid. Using the definitions 
given in Eqs. (F.14) and (F.47), one can write Eq. (4.42) at the final time T f as: 


1 _ l 2A[AH-B(sinhH f -sinhH 0 )] 


Pf Po 


(F.50) 


By defining the quantity D as: 


D = b + 2 Ap 0 [ AH - B(sinh H f - sinh H 0 )] (F.5 1 ) 


p f can be written from Eq. (F.50) as: 


Pf = 


bp 0 

D 


(F.52) 


The last step consists of substituting from Eqs. (F.48), (F.49) and (F.52) into the 
expression for A given in Eq. (F.21). One can write the result as: 


A = — 


16A 2 h*(B 2 -l) 

D 2 


sinh 2 




(F.53) 


Since the transfer duration is strictly positive, that is, T f > 0, it must necessarily be true that 
A0 = 0 f - 0 O > 0, which in turn implies that A£ = £ f - £ 0 > 0 and AH > 0. Moreover, for 
AH > 0 the following inequality is always valid: 



(F.54) 
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Accordingly, from Eq. (F.53), the discriminant A must be strictly negative (A < 0), 
contradicting the necessary condition A > 0 found in Eq. (F.7). This concludes the proof of 
part (i) for B > 1 . 

At this point, the proof of part (i) of Theorem 4.1 is complete. Now that part (i) has 
been proven, the validity of part (ii) becomes intuitively obvious, and can be argued for in 
the following way: First, assuming that the initial and final conics are tangent to each other, 
one will have to agree that the existence of an impulsive transfer solution at the point of 
tangency is obvious. All one has to do for such an impulsive transfer is use the correct 
amount of et f at the point of tangency (with e infinitely large, x f infinitely small, and £T f 
finite) and change the speed of the vehicle (tangentially to both conics) by the desired 
amount (AV). That there can be no other impulsive transfer is also obvious, since the 
orbital speeds associated with the two conics at the point of tangency are unique, resulting 
in a unique AV requirement. Thus, the only thing that is left to complete the proof for part 
(ii) is to show that there can be no finite time (nonimpulsive) transfer in such a case. But 
this is again (almost) obvious, since, for a finite time transfer, the discriminant A found in 
Eqs. (F.29), (F.37), (F.45), and (F.53) is strictly negative, contradicting the requirement 
that it be zero (see Eq. (F.7)), coming from the assumption that the initial and final conics 
have one common point. 

Q. E. D. 
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Appendix G 


Explicit Expressions for the Power-Limited Cost 
Along a One-Segment Transfer Trajectory 

In this Appendix the goal is to obtain a workable approximation for the power- 
limited cost associated with a one-segment transfer trajectory, valid in the limit of long 
transfer duration and very low thrust (small throttling parameter A). Recall the expression 
for the (nondimensional) speed of the vehicle given by Eq. (A. 8) in Appendix A. 
Differentiating this expression and 1/r with respect to time, and using Eqs. (4.1 1) through 
(4.13) one obtains : 


dV_ AY x _d_fn _ 2 L 

dx h Vr 2 ’ dxlrj r 2 


(G.l) 


Let now S be the (nondimensional) mechanical energy (kinetic plus potential) of the 
vehicle. Explicitly, S is: 


S 



(G.2) 


Using Eq. (G.l), and substituting for the thrust acceleration program from Eq. (4. 159), the 
time derivative of S can be found as: 


dS 

dx 



(G.3) 


Equation (G.3) expresses a familiar result. It simply states that the work done per unit time 
t>y the thrust program proposed in Eq. (4.10) is equal to the time rate of change of the 
mechanical energy of the vehicle. Along a (one-segment) transfer trajectory dh = Adx, so 
that Eq. (G.3) can be rewritten as: 


dS = 


rvM 


dh 


(G.4) 
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The power-limited cost associated with a (one-segment) transfer is given by: 


'PL 


. teh . itevidt - 4 £ dh 

J 2 2J h 2 2 J h 2 


(G.5) 


and, using Eq. (b .4), it assumes the form: 


s, 


J = *1^ 

PL 2 J h 

So 


(G.6) 


Integrating by parts in Eq. (G.6) one obtains: 


= Afs L _0 A fSdh 
PL 2 (h f h 0 , 


j 


2 J h 


(G.7) 


and, substituting from Eq. (G.2), J PL can be written as: 


J = — 
Jpl 2 


'S.-Sa.'l+A.f 

‘ V 2 dh 

_±f 

dh 

^h f h 0 J 4 J 

h 2 

2 J 

rh 2 


(G.8) 


Using the definition of J PL given in Eq. (G.5), one can further simplify Eq. (G.8) as: 



(G.9) 


From Eqs. (A.7) and (G.2) S is equal to -l/(2a), where, a is the instantaneous 
semimajor axis along the transfer trajectory. Combining this with Eq. (A.9), one obtains: 


S~- 

2a 2h 2 


(G.10) 


where of course e and h are the (instantaneous) eccentricity and angular momentum along 
the transfer trajectory respectively. Using Eq. (G.10), one can write Eq. (G.9) as: 


130 



J - A K-l eg-l 
PL 2( h f 3 hj 


-aK 

J rh 2 


Consider now a (one-segment) transfer trajectory on which B is zero. On such 
transfer trajectory r is equal to (see Eq. (4. 17)), and Eq. (G.l 1) becomes: 


if B = 0, then: J P1 = — 


_ A ( e 2 - 1 e 2 -l 


21 h;* 


Of 


which after a simple integration results in: 


if B = 0, then: J p , = —\lA — I _ M 1 


6 V h f 3 


When B is nonzero things are more complicated. With £=0-C (see Eq. (4.20)), and 
d£=d0, one can now use Eqs. (4. 1 1) and (4. 14) and write Eq. (G.l 1) as: 


_ A [e 2 - 1 e 2 -l 


2 l h? h^ 


rhl'i 


(G.14) 


When B=l, one can introduce the variable y=tan(£/2), and, noting that: 

^ = ; r = — - — = -. h2 = ^ ! . ( 1 + y2 ) 

(l + y ) l + Bcos^ 1 + cos^ 2 


(G. 15) 


one can substitute in Eq. (G.14) and write the power-limited cost as: 


if B = 1, then: 


t = AfeL_!_£ozi 




(G.l 6) 


Similarly, when 0 < B < 1, or B > 1, one can introduce respectively the generalized 
eccentric and hyperbolic anomalies, use the results of Sections 4.5.2 and 4.5.3, and write 
the power-limited cost as: 
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if 0 < B < 1, then: 


JpL _ 


fe?- 1 ej-l) 

- ^ 

h ? h o J 

Vl-B 2 J 


(G.17) 


if B > 1, then: 


j = A 

Jpl 2 


ft M ej-l) 

n 

A 2 | 

hf K J 

Vb 2 -i J 

u 


dH 

h 


(G.18) 


Although the expression for the power-limited cost found in Eq. (G. 13) is exact only 
when B=0, this expression is also approximately valid when B is small: 


if B = 0, then: 


Jpl ~ 


\-h 0 
, 6x f 


v 3e 2 - 1 

h? 


3eg-f 

H J 


(G.19.a) 


For transfers from a circular orbit (see Section 4.8) e 0 is zero, and if e f is not very 
large, then B is quite close to zero (see Figs. 4.7, 4.13, 4.15, 4.17, and 4.19). In this case 
Eq. (G.19.a) assumes the form: 


if e 0 =0, B = 0, then: 


Jpl ~ 


'Vi-. v 


6x f 


1 3e f 2 - 1 

W + ~ir- 


A“o 






J 


(G.19.b) 


Transfers between two coplanar circular orbits (see Appendix D) correspond to the 
limiting case of very small e f . On such transfers B is indeed very small (see Fig. 4.7), and 
in this case Eq. (G.19.b) assumes the form: 


if e 0 = e f = 0, B = 0, then: 


Jpl = 


, 6t f 


v 


A h 0 



(G.19.C) 


The validity of the approximation expressed in Eq. (G.19) is checked in Table G.l 
that follows, by comparing the power-limited cost evaluated from Eq. (G.19.b) with the 
exact power-limited cost (see Table 4.1) for some of the transfers (the ones from a circular 
orbit). The second and third columns of Table G.l supply the throttling parameter A and 
the generalized eccentricity B for the corresponding transfer. The approximate power- 
limited cost, evaluated from Eq. (G.19.b), is given in the fourth column, and it is 
compared with the exact power-limited cost given in the fifth column, reproduced from the 
fourth column of Table 4.1. 
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Table 

G.l ADDroximate 

Dower-limited 

cost usinn Ea. 

(G.191 

Transfer 

A 

B J 

PL (Eq. (G.19)) 

J PL (Table 4.1) 

Fig. 4.7 

6.1 139e-3 

1.2228e-2 

8.9257e-4 

8.9225e-4 

Fig. 4.13 

3.6210e-2 

7.2420e-2 

5.5579e-3 

5.4860e-3 

Fig. 4.15 

1.0542e-l 

2.1084e-l 

2.0017e-2 

1.7064e-2 

Fig. 4.17 

1 .01 3 le- 1 

2.0261e-l 

2.2436e-2 

2.0457e-2 

Fig. 4.19 

2.960 le-2 

5.9202e-2 

5.0924e-3 

5.0593e-3 


It is apparent from Table G.l that the approximation expressed by Eq. (G.19.b) is 
indeed valid for low thrust levels (small A) and small B. Thus, for transfers from a circular 
to a nearly circular orbit, since both A and B are small, Eq. (G.19.b) is valid. For low- 
thrust transfers however from an arbitrary conic, B in general need not be small (see Fig. 
4.2 for example), and Eq. (G.19) needs to be corrected. For the cases in which 0 < B < 1, 
or B > 1, this correction can be obtained by noting that Eqs. (4.34) and (4.42), supplying 
the dependence of the angular momentum h on the generalized eccentric and hyperbolic 
anomalies E and H respectively, can be approximated as: 


small A, 0 < B < 1: 


1 _ 1 2 A(E — E 0 ) 

h 2 ~ h^ (1-B 2 ) 3 ' 2 


(G.20) 


small A, B > 1: 



(G.21) 


Now, substituting from Eqs. (G.20) and (G.21) into Eq. (G.17) and (G.l 8) respectively, 
one obtains for both of the above cases the approximate expression: 


if A is small enough, then: 


J PL = 


\-h 0 Y 3e 2 -2B 2 - 1 

, 6^ hf 



(G.22) 


Note that Eq. (G. 19) corresponds to a special case in Eq. (G.22), for which B=0. By 
a continuity argument it is also obvious that the approximation expressed in Eq. (G.22) 
should also be valid when B=l. Thus, Eq. (G.22) is a single, general approximation, valid 
for low-thrust, long duration, one-segment transfers. 
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One must be careful however when using Eq. (G.22), because, for the case in which 
both the initial and final conics are noncircular, low-thrust, long duration, one-segment 
transfers between the conics may not exist. It appears that, when 0 < B < 1, Eq. (G.22) is 
valid for low-thrust, many revolution, one-segment transfers, that is, transfers on which 
E f -E 0 is large. When on the other hand B > 1, Eq. (G.22) is most probably valid for low- 
thrust, one-segment transfers on which H f -H 0 is small. This condition seems to be 
necessary if the term -B(sinhH-sinhH 0 ) that was omitted when writing down Eq. (G.21) is 
to be unimportant along a transfer. Table G.2 that follows supplies all the cases of transfers 
among the examples, for which the power-limited cost evaluated from Eq. (G.22) (fourth 
column) compares well with the exact power-limited cost (fifth column) associated with the 
same one-segment transfer. 


Table 

G.2 ADDroximate 

power-limited 

cost usine Ea. 

(G.22) 

Transfer 

A 

B 

J PL (Eq. (G.22)) 

J PL (Table 4.1) 

Fig. 4.2 

3.5000e-3 

7.0000e-l 

1.5969e-4 

1.5369e-4 

Fig. 4.4 

3.5000e-2 

7.0120e-l 

3.9261e-3 

3.9537e-3 

Fig. 4.7 

6.1139e-3 

1.2228e-2 

8.9283e-4 

8.9225e-4 

Fig. 4.13 

3.6210e-2 

7.2420e-2 

5.6132e-3 

5.4860e-3 

Fig. 4.19 

2.960 le-2 

5.9202e-2 

5.1264e-3 

5.0593e-3 


For the cases of course in which B > 0 and the approximation supplied by Eq. (G.22) 
is not good one can always evaluate the power-limited cost associated with a one-segment 
transfer trajectory using the exact expressions given by Eqs. (G.16) through (G.18). 
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Appendix H 


Comparison Between Keplerian Motion and the 
Motion Uncovered in Section 4.4 

Keplerian conics corresponding to Keplerian motion (see Appendix A), and two- 
dimensional transfer trajectories, corresponding to the motion uncovered in Section 4.4, 
have in this report been unified into a single descriptive scheme, exemplified by a common 
trajectory equation. This Appendix supplies an elementary comparison of the two types of 
motion, which can be summarized as in Table H.l that follows. 

Table H.l Comparison b etween Keplerian motion and Present motion 

Type of motion: Keplerian Present (Section 4.4) 

Thrust acceleration: e = 0 e = — — (tangential) 

rcosy 

Angular momentum: h = h 0 = constant h = Ax + h 0 

Trajectory equation: r = r - !l! 

1 + e cos(0 - to) 1 + e cos(0 - to) 

In Keplerian motion the thrust acceleration is zero and the angular momentum of the 
vehicle (considered as a point mass) about the planet is constant. In the type of motion 
uncovered in Section 4.4 the thrust acceleration is tangential and equal to A/(rcosy), where 
A is a constant (the throttling parameter), and the angular momentum of the vehicle about 
the planet is a linear function of time. For both types of motion the trajectory equation has 
the same identical form! Just as Keplerian motion, the motion uncovered in Section 4.4 
corresponds to an exact analytic solution of the equations of motion. Note that the 
generalized eccentricity and orientation constants B and C of a transfer trajectory have in 
Table H.l been denoted as e and to respectively, to make clear the analogy with the 
eccentricity and orientation constants corresponding to Keplerian motion. The reader is 
advised to compare Table H.l with Table 1.1 (page 139) of Appendix I. 
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Appendix I 


The Keplerian Class of Thrust Programs 

Let Q be an arbitrary, explicit, differentiable function of the nondimensional time 
variable x. For any such function Q one can form the thrust program given explicitly by the 
components: 

xQ 2rQ . Q (U) 

6 '- h h ’ * r 

where the dot' denotes differentiation with respect to x. Q will be called the throttling 
function. Note that when Q is a constant, say equal to a throttling parameter A, one obtains 
from Eq. (1.1) the thrust program proposed in Eq. (4.10) of Section 4.3 in the main text. 
The totality of thrust programs defined by Eq. (1.1) will be given the name the Keplerian 
class. The Keplerian class is a rather privileged class of thrust programs, because for all its 
members the corresponding equations of motion afford a rather simple exact analytic 
solution. Specifically, if the right-hand-sides of Eqs. (4.6), (4.7) are forced with the 
components of the thrust acceleration vector given in Eq. (1.1), then the exact solution of 
Eqs. (4.6) through (4.9) is given by the following four quadratures: 


h = h 0 + j*Qdx 

0 

(1.2) 

2rQ + Bsin(0-C) 

(1.3) 

X “ h 


h 2 

(1.4) 

1 + Bcos(0-C) 

fdx f d0 

J h 3 J [l + Bcos(0-C)f 
0 6 0 

(1.5) 

where B, C, and 0 O are four arbitrary constants of integration. This result can be easily 
verified by differentiating Eqs. (1.2) through (1.5) with respect to x, taking into account Eq. 
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(1.1), and by checking that Eqs. (4.6) through (4.9) are identically satisfied. Eqs. (1.2) 
through (1.5) reduce to Eqs. (4.15) through (4.18) when Q=A=constant. Note that if Q is 
not explicit in T, then Eqs. (1.3), (1.4) still constitute two exact quadratures, while obtaining 
two more quadratures in such a case may or may not be possible depending on the 
functional dependence of Q. Note also that again, as in Eq. (4.17), Eq. (1.4) describing the 
trajectory is identical in form with the trajectory equation corresponding to two-body 
Keplerian motion! Thus, the constants B and C can again be called the generalized 
eccentricity and the generalized orientation respectively, since they play qualitatively the 
same role as the eccentricity and orientation constants e and go do in Keplerian motion. 

Apart from its kinship with Keplerian motion, the importance of the motion arising 
under the Keplerian class of thrust programs comes partly from the following claim 55 : If the 
throttling function Q is small, and if the time derivatives of the throttling function d n Q/dx n 
n— 1>2,3... are much smaller than the throttling function itself, then the corresponding 
motion under Keplerian thrust, described by Eqs. (1.2) through (1.5), can presumably be 
expected to be nearly-optimal when performed with a power-limited propulsion system. 
This claim can be argued upon by examining 55 the costate system of Eqs. (4.54) through 
(4.56), and the transversality conditions, Eqs. (4.106), (4.109), (4.110), under the 
assumption that the optimal power-limited thrust is approximately Keplerian. Note that for 
long duration maneuvers one can easily find members of the Keplerian class for which the 
above conditions on the throttling function Q are satisfied 55 . In particular, one can introduce 
four classes of such members each one of which is suitable for treating a special type of 
finite-time problem that may be encountered during orbital operations. These four classes 
will be given the names the Tangential , the Linear , the Quadratic, and the Cubic class, and 
correspond to the (Keplerian) thrust programs for which the throttling function Q is 
respectively a zero, first, second, and third degree polynomial of time t. 

A thrust program will be called tangential if and only if the corresponding thrust 
acceleration vector is tangent to the vehicle's flight path for all time T for which the thrust 
program is defined. 

Note now that: 

(i) The throttling function Q, and the corresponding thrust program, depend on one, 
two, three, and four arbitrary constants for the Tangential, the Linear, the Quadratic, and 
the Cubic (sub)classes of (Keplerian) thrust programs respectively. 

(ii) When dQ/dt is small compared to Q, that is, when Q is slowly varying, then the 
corresponding thrust is nearly-tangential. 
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(iii) A Kepleriart thrust program is tangential if and only if it is a member of the 
Tangential (sub)class. 

(iv) For the Tangential class of thrust programs Q is a zero degree polynomial, 
namely a constant. 

(v) The class of thrust programs proposed in Eq. (4.10) of Section 4.3 was none 
other than the Tangential (sub)class of (Keplerian) thrust programs. 

(vi) Chapter IV of this report treated in some detail the motion that can be executed 
under the Tangential class of thrust programs. 

It can be verified 55 that near-optimality (in the power-limited sense) and controllability 
appear as two opposing trends as one goes from the (simplest) Tangential class toward the 
(most sophisticated) Cubic class of thrust programs. One can satisfy the first-order 
necessary conditions for optimality extremely well using the Tangential class (see Section 
4.6.5), but the vehicle is not completely controllable under the corresponding thrust 
programs (see for example Theorem 4.1 of Section 4.11, and Appendix F). On the other 
hand, using the Cubic class one may not be able to satisfy the first-order necessary 
conditions as well, but the vehicle is much more controllable under the corresponding 
thrust programs 55 . In fact, it is straightforward to show the following theorem about 
controllability with regard to the Cubic class (the proof is given in Ref. 55): 

Theorem. If 1^, h^ r 0 , r f , and t f - x 0 are all strictly positive , then, any given initial 
state (h,),x o ,r o ,0 o ) at a given initial time x=x 0 can be driven to any given final state 
(h^x^rpGf) at a given final time x=x f by a thrust program belonging to the Cubic class. 

The four subclasses introduced above are suitable for treating the following finite- 
time problems (Note that 1^, e 0 , co 0 , e f and co f are always fixed): 

(a) Escape problems: For all such problems 0 O , 0 f , and h, are all left free. The 
time duration x f of the maneuver is fixed. The throttling function Q should depend on (at 
least) one arbitrary constant. The Tangential class of thrust programs for which the 
throttling function is an arbitrary constant is uniquely suitable for treating these problems. 
In this report this was done in Sections 4.9 and 4.10. 

(b) Transfer problems: For these problems 0 O and 0 f are left free, but h* is fixed. 
The time duration x f of the maneuver may be fixed or free. If x f is left free then Q should 
depend on (at least) one arbitrary constant, meaning that the Tangential class is the suitable 
one to choose in such cases. This was done in this report in sections 4.8 and 4. 10. If x f is 
fixed then Q should depend on (at least) two arbitrary constants and the Linear class may be 
used. 
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(c) Mixed Transfers-Rendezvous: For these problems h, is fixed, and only one 
among 0 O or 0 f is left free, while the other is fixed. The time duration x f of the maneuver 
may be fixed or left free. If x f is left free then Q should depend on (at least) two arbitrary 
constants, meaning that the Linear class is the suitable one to choose in such cases. If x f is 

fixed then Q should depend on (at least) three arbitrary constants and the Quadratic class 
may be used. 

(d) Rendezvous Problems: For these problems 0 O , 0 f , and h f are all fixed. The 
time duration x f of the maneuver may be fixed or left free (in practice it is usually desirable 
to have a fixed x f ). If x f is left free then Q should depend on (at least) three arbitrary 
constants, making the Quadratic class the suitable one to choose in such cases. If x f is fixed 

then Q should depend on (at least) four arbitrary constants and the Cubic class may be 
used. 

This Appendix concludes by supplying a comparison between Keplerian motion and 
the motion arising under Keplerian thrust, which can be summarized as in Table 1.1 that 
follows. The reader is advised to compare Table 1.1 with Table H.l (page 135) of 
Appendix H. 


Table 1.1 Compa rison between Keplerian motion 
and motion under Keplerian thrust 


Type of motion: 

Keplerian 

Present 

Thrust: 

Zero 

Keplerian 

Thrust acceleration: 

II 

o 

00 

© 

II 

o 

o _ xQ + 2rQ 

h ’ e " 

Angular momentum: 

h = h 0 = constant 

X 

h = h 0 + jQdx 

0 

Trajectory equation: 

h 2 

r — 

h 2 

l + ecos(0 -co) 

r — 

l + ecos(0-co) 


Note once more that the generalized eccentricity and orientation constants B and C in 
iq. (1.4) have in Table 1.1 been denoted as e and © respectively, to make clear the analogy 
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with the eccentricity and orientation constants corresponding to Keplerian motion. The 
detailed theory introducing the Keplerian Class of thrust programs is given in Ref. 55. For 
applications corresponding to the Keplerian Class of thrust programs the reader is referred 
to Ref. 56 and to the proceedings of the AIAA Astrodynamics Conference of 1995. 
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Figure 2. 1 Energy climb paths for an F-8 aircraft. 
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Figure 2.2 Evaluation of £(E) for an F-8 aircraft. 
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Figure 2.3 Energy climb paths for an F-15 aircraft. 



Figure 2.4 Evaluation of e(E) for an F-15 aircraft. 
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Figure 2.8 Evaluation of e(E) for a generic hypersonic vehicle. 
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vehicle path 



Figure 4.1 Polar coordinates r, 0, and nondimensional state variables 
h, x, r, and 0 (T=nondimensional time). 
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Figure 4.2 Example of transfer between two elliptic orbits. 

El: ho=l, e 0 =0.7, co 0 =0°, E2: 1.2568, epO.7097, (0,=0.195 o , 

T: A=0.0035, B=0.7000, C=O.337O°,0 o =O 0 , 9,= 1 160°, 1^73.361. 



Figure 4.3 Example of transfer between two elliptic orbits. 
El: 1^=1, e 0 =0.7, to 0 =0°, E2: h^l.4114, e^O.6468, ©,=0.731°, 
T: A=0.015, B=0.7002, C= 1.444°, 0 o =O°, 0,=28O O , t,=27.426. 
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Figure 4.4 Example of elliptic to hyperbolic transfer. 

El: ho=l, e 0 =0.7, co 0 =0°, H2: hp2.1077, e^l.0768, ©,=66.547°, 
T: A=0.035, B=0.7012, C=3.366°, 6 0 =0°, 0^196°, x,=3 1.648. 



Figure 4.5 Example of hyperbolic to hyperbolic transfer. 

HI: 1^=1, e 0 =1.01, © 0 =0°, H2: h^l.3951, e^l.5071, ©,=344.030°, 
T: A=0.03, B=1.3164, C=341.640°, 0 o =-15O°, 0,=9O°, x,= 13.169. 
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Figure 4.8 Solution search for the same circular to circular transfer 
as in Fig. 4.6, but with smaller e f . CO: 1^=1, e 0 =0, Cf: ty=2, e^O.01, co^O 



Figure 4.9 Solution search for the same circular to circular transfer 
as in Fig. 4.6, but with smaller e f . CO: 1^=1, e 0 =0, Cf: h,=2, e,=0.005, co^O 0 





Figure 4. 10 Solution search for the same circular to circular transfer 
as in Fig. 4.6, but with smaller e f . CO: 1^=1, e 0 =0, Cf: hf=2, epO.OOl, a>(=0°. 



Figure 4. 1 1 Solution search for the same circular to circular transfer 
as in Fig. 4.6, but with smaller e f . CO: 1^=1, e 0 =0, Cf: hf=2, e ,=0.0005, co,=0 0 . 
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Figure 4. 12 Solution search for a circular to elliptic transfer. 
CO: h 0 =l, e 0 =0, Ef: h^2, e^O.35, <0^=0°. 



-5 0 

Figure 4.13 Transfer corresponding to point S (with k=2) in Fig.4.12. T: A=3.6210X10 
B=7.2420X10- 2 , C=309.249°, 0 O =219.249°, 0,,,= 100.446°, t,=27.617. 






Figure 4.15 Transfer corresponding to point S (with k=l) in Fig. 4.14. 

T: A=0. 10542, B=0.21084, C=353.501°, 0 o =263.5Ol°, 0 fP =92.785°, Xf=4.143. 
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Figure 4. 16 Solution search for a circular to hyperbolic transfer. 
CO: 1^=1, e 0 =0, Hf: h,=2, e,=l.l, ©pO 0 . 



Figure 4.17 Transfer corresponding to point S (with k=l) in Fig. 4.16. 

T: A=0.10131, B= 0 . 20261 , C=312.109°, e o = 222 . 109 °, 0 fp =98.861 0 , X(=9.871 





Figure 4. 18 Solution search for an escape from a circular orbit at T^IOO. 
CO: 1^=1, e 0 =0, Pf: h^free, e,=l, copO 0 . 



Figure 4.19 Transfer corresponding to point S (with k=3, 1^3.9601) in Fig. 4.18. 
T: A=2.9601X10' 2 , B=5.9202X10 2 , C=346.123°, 0 O =256.123 O , 0fl.F9O.863 0 , tplOO. 
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Figure 4.20 Solution se 
E0: 1^=1, e 0 =0.9. 



T: A=7 


Figure 4.2 1 Transfer co 
.3577X10' 3 , B =0.87362, 


rrespon 
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Figure 4.22 Solution search for an elliptic to elliptic transfer involving a small 
change in eccentricity. E0: 1^=1, e 0 =0.9, o) 0 =O°, Ef: h,=2, epO.895, tty=0 . 



Figure 4.23 Solution search for (practically) the same elliptic to elliptic transfer 
as in Fig. 4.22. E0: 1^=1, e 0 =0.9, (D 0 =0°, Ef: 11^2, e,=0.898, CDpO . 
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Figure 4.24 Solution search for an elliptic to elliptic transfer with orientation change 
EO: h^l, e 0 =0.9, co 0 =0°, Ef: 1^2, e,=0.8, 0)^30°. 



Figure 4.25 Transfer corresponding to point S (with k=0) in Fig. 4.24. 

T: A=1.5121X10' 2 , B=0.83962, C=352.146°, 0 O =149.252°, 0*= 195. 104°, v=66.134. 
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Figure 4.27 Transfer corresponding to point S (with k=l) in Fig. 4.26. 

T: A=2.0519X10' 2 , B=0.86914, C=349.330°, 0 O =198.586°, e,p=l 53.967°, x^AZ.136. 
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Figure 4.28 Solution search for a hyperbolic to hyperbolic transfer. 
HO: h 0 =l, e 0 = 1.1, co 0 =0°, Hf: 1^2.3, e,=3.4, ©,=320°. 



Figure 4.29 Transfer corresponding to point S (with k=0) in Fig. 4.28. 

T: A=5.2151X10- 2 , B=2.8264, C=319.870°, 6 0 =-149.782°, e^SO.642 0 , x f =24.928. 
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Figure 4.30 Escape from an elliptic orbit at i,=200, first solution (© f is free). 
E0: ^=1, e 0 =0.8, © 0 =0°, Pf: hp3, e,=l, ©,=0.3562°, T: A=0.01, B=0.76944, 
C=355.827°, 6 0 = 149.787°, 0*= 104.974°, k=2, x,=200. 



Figure 4.31 Escape from the same elliptic orbit (as in Fig. 4.30) at Xp200, second 
solution (© f is free). E0: 1^=1, e 0 =0.8, © O =0°, Pf: h^3, e,=l, ©,=5 1.459°, T: A=0.01, 
B =0.78853, C=0.5593°, 6 0 =56.212°, 6 n >=192.057°, k=l, x,=200. 
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Figure 4.32 The two associated conics (dotted lines) at the initial and final times, 
for the transfer given in Fig. 4.21. ACO: h^l, e 0 =0.87362, © 0 =358.544°, 
ACf: hf=2, epO.87362, ©,=358.544°. 



Figure 4.33 The two associated conics (dotted lines) at the initial and final times, 
for the transfer given in Fig. 4.25. ACO: 1^=1, e 0 =0.83962, © 0 =352.146°, 
ACf: hf=2, epO.83962, ©,=352. 146°. 
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Figure 4.34 The geometry of departure from a circular orbit (Lemma 4.3). 



Figure 4.35 The geometry of departure from an arbitrary conic (Lemma 4.4). 
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Figure 4.36 Two hyperbolic trajectories having opposite orientation. 
A one-segment (transfer) solution cannot exist in such a case. 



Figure 4.37 A two-segment transfer between two hyperbolic orbits for a case similar 
to the one given in Fig. 4.36. HO: h 0 =1.7604, e 0 =1.2297, co 0 =180°, Hf: h^l.7604, 
ef 3 1-2297, copO 0 , Tl: A=-0.15, B=0.3, C= 1 80°, 6 0 =90°, 0 < =27O°, 1^5.0690, T2: 
Symmetric to Tl about 0=90° with A=0.15 (Note: In Table 4.1, Section 4.13, the costs are 
the total costs, while £<, and e f correspond only to T2). 
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Figure 4.38 An example of a three-segment (thrust(Tl)-coast(C2)-thrust(T3)) transfer 
for changing the orientation of an elliptic orbit by 90 degrees. E0: h^l, e 0 =0.9, to 0 =O°, Ef: 
h^l, e,=0.9, to ,=90°, Tl: A=6.3771X10 2 , B=1.4030, C=302.441°, 0 o =172.9O4°, 
0 K =2 12.441°, t ok = 36.327, T3: Symmetric to Tl about 0=225° with A=-6.3771X10 2 , C2 
(dotted line) is just a circular (Keplerian) segment with A=0 (Note: In Table 4.1, Section 
4.13, the costs are the total costs, while e 0 and e f correspond only to Tl). 





Figure 4.39 Variation of one-half times the square of the (nondimensional) thrust 
acceleration with (nondimensional) time for the transfer example of Fig. 4.2. 



Figure 4.40 Variation of one-half times the square of the (nondimensional) thrust 
acceleration with the argument of latitude for the transfer example of Fig. 4.2. 
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Figure 4.41 The Hohmann transfer (H) between two coplanar 
elliptical orbits having the same orientation. 



Figure 4.42 The Biparabolic transfer (PI), (P2) between two coplanar 
elliptical orbits having the same orientation. 
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